similar to: Anova Help?

Hello dear R members. I have been learning the Anova syntax in order to perform an SS type III Anova with repeated measures designs (thank you Prof. John Fox!) And another question came up: where/what are the (between/within) residuals for my model? ############ Play code: phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),

(Have searched for this already) Hi, How do you find the strength of correlation between two variables using an ANOVA table? "Pr(>F)" gives the statistical significance of the association, but not the strength of the correlation. See data (from R) below Readable: "Df" "Sum Sq" "Mean Sq" "F value"

Hello everyone, I am doing an analysis on some data for a project, and do not have much background in either R or statistics, but have been pointed in the direction of doing ANOVA and then multiple comparisons on the data to give me the information I need. I was wondering if some of you could help me out. What I want to do is to do ANOVA on a data set like the following: ID G1 G2 G3 NG1 NG2

I created a 2 D plan: > plot(range(-2.5,0.95),range(0.00,1.00),type="n",axes=TRUE) I made a projection of points with their coordonates (X,Y) in that plan > fullpoints=read.csv2("fullpoints.csv",h=T) > plot(fullpoints) The points are listed in that .csv file it is organized this way: 0,48875 0,142857143 0,409 0,142857143 0,45611 0,25 0,49833 0,222222222 0,61158

Hello dear R-help group. My task looks simple, but I can't seem to find a "smart" (e.g: non loop) solution to it. Task: I wish to randomize a data.frame by one column, while keeping the inner-order in the second column as is. So for example, let's say I have the following data.frame: xx <-data.frame(a= c(1,2,2,3,3,3,4,4,4,4) , b =

Hello R help members. I came across today with an article on Barnard's exact test (http://www.cytel.com/Papers/twobinomials.pdf), that is supposed to give a more powerful fisher.test - Because it doesn't assume that we know the row and column totals are in advance. Any pointers to such a function ? Thanks, Tal -- ---------------------------------------------- My contact information:

When I perform a two-way anova on my dataframe "pin", I can't get any indication about the interaction among the two factors "gen" and "con" while the statistics about the significance relative to the two factors alone is correct. I wrote this: *pvalues_genotipo<-sapply(pin[,1:(length(pin)-2)],FUN=function(x) summary(aov(x~gen*con)))* Where is the mistake?

Hi All, I am new comer here and for Statistics R. I want to know how to use R to calculate one way repeated measure ANOVA, for example, on three groups of data. I appreciate if someone can help me solve this problem. John

I have the attached data set (csv) and I want to run an analysis of variance on the wingsize data (comparing infected vs non-infected) within and between experiments. Can anyone help me with the command I should use? Many thanks Georgina PhD Student Division of Infection and Immunity B5-29, GBRC 120 University Place Glasgow G12 8TA Tel: 0141 330 5650

Hello!!! I'm trying to do multiple ANOVA tests with R (testing the affect off different factors on the same response). As a result I get many ANOVA tables, and I want to extract a list of the "Pr(>F)" from all the tables. Maybe someone have an idea how to do this? Thanks Imri -- View this message in context:

*(note*: This is an R community question, not a statistical nor coding question. Since this is my first time writing such a post, I hope no one will take offence of it.) Hello all, I will be attending useR 2009 next week, and was wondering if there are any of you who are *bloggers *intending to participate and report on useR 2009? If so - I would love to know your blogs URL so as to follow you.

Hi, I'm trying to use the function "error.bars", but the program don't find it, and I dont't found any package with this function. Is there some another functin to draw barplots with error bars? Sueli Rodrigues Eng. Agr?noma - UNESP Mestranda - USP/ESALQ PPG-Solos e Nutri??o de Plantas Fones (19)93442981 (19)33719762

Hello Martin Maechler and All, A simple question (I hope): How can I compute the "sum of the dissimilarities" that appears in the pam command (from the cluster package) ? Is it the "manhattan" distance (such as the one implemented by "dist") ? I am asking since I am running clustering on a dataset. I found 7 medoids with the pam command, and from it I have the

Hi all, I am trying to run a two factor anova, but one of the factors is a random factor, now I am also running in SPSS and it seems its dividing by the wrong term to get the appropriate F term. here is my data. In SPSS the F scores about double the ones in R, how can I specify one of my factors as a random factor or change it to where it does the right model fitting? I am using the lm command

Hello dear R mailing list members. I have recently became curious of the possibility applying model selection algorithms (even as simple as AIC) to regressions of large datasets. I searched as best as I could, but couldn't find any reference or wrapper for using step or stepAIC to packages such as biglm. Any ideas or directions of how to implement such a concept ? Best, Tal --

Hello all. I wish to run k-means with "manhattan" distance. Since this is not supported by the function "kmeans", I turned to the "pam" function in the "fpc" package. Yet, when I tried to have the algorithm run with different starting points, I found that pam ignores and keep on starting the algorithm from the same starting-points (medoids). For my

(I use R 2.9.1 with win XP) If I run this code: seq(-0.1,.9, by = .05)[seq(-0.1,.9, by = .05) <= 0.5] I get this output: [1] -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 Why is 0.50 not in the results ? (It seems that it gives a slightly bigger number then 0.5 but I don't understand why it does that) Where as if I try: seq(-0.1,.9, by = .05)[seq(-0.1,.9, by = .05) <=

Hello Prof. Harrell and dear R-help mailing list, I wish to perform a non-parametric repeated measures anova. If what I read online is true, this could be achieved using a mixed Ordinal Regression model (a.k.a: Proportional Odds Model). I found two packages that seems relevant, but couldn't find any vignette on the subject: http://cran.r-project.org/web/packages/repolr/

Hello dear R help group, I am trying to read a .txt file, with Hebrew column names, while keeping the column names looking well in R - but without success. I uploaded an example file to: http://www.talgalili.com/files/aa.txt And am trying the command: read.table("http://www.talgalili.com/files/aa.txt", header = T, sep = "\t") This returns me with: X.....ª X...ª......

Hello dear R help group, I would like to download the tRNA data on: http://gtrnadb.ucsc.edu/download.html And then import it into R. Can anyone direct me as to how to do so? Thanks, Tal ---------------------------------------------- My contact information: Tal Galili E-mail: Tal.Galili@gmail.com Phone number: 972-52-7275845 FaceBook: Tal Galili My Blogs: http://www.talgalili.com (Web and