search for: 3373767

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

...,] new.xx <- rbind(new.xx , xx.subset) } return(new.xx) } randomize.by.column.a(xx) # END of - R example I would love for a better, faster, way of doing it. Thanks, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...80 0.3256 0.992814 ## gender:phase:hour 6.636 8 96.167 80 0.6900 0.699124 ## treatment:gender:phase:hour 14.155 16 96.167 80 0.7359 0.749562 -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: www.talgalili.com www.biostatistics.co.il [[alternative HTML version deleted]]

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

...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: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...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 -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: www.talgalili.com www.biostatistics.co.il

...URL so as to follow you. I am planning on coming with a laptop and see what I'll be able to report on: www.r-statistics.com (which is very empty for now) See you there, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...rm(15,5,0.5))) pam(x, 2, medoids = c(1,16)) output: Medoids: ID [1,] 3 -0.1406026 0.1131493 [2,] 17 4.9564839 4.6480520 ... So the initial medeoids where 3 and 17, not 1 and 16 as I asked. Thanks, Tal -- ---------------------------------------------- Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: www.talgalili.com www.biostatistics.co.il [[alternative HTML version deleted]]

...the minimum manhattan distance to the medoids that pam predicted. If this is the manhattan distance that is used, I will create some toy data to see if I can reproduce this. Thanks, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...as if I try: seq(-0.1,.9, by = .05)[seq(-0.1,.9, by = .05) <= 0.4] and get: [1] -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 Then 0.40 WILL be in the results. Thanks,Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

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

...never ran R in "batch", and any "step by step" instructions (either by referring to links or explaining here) would be of great help. Thanks in advance, Tal ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...split(temp.file[i], "\t"))) } new.df # or maybe apply(as.data.frame(temp.file), 1, function(b) unlist(strsplit(b, "\t"))) # but this transposes the matrix Thanks, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...d ( Stuart-Maxwell<http://ourworld.compuserve.com/homepages/jsuebersax/mcnemar.htm>) ? 3) Do you have a recommended link to an explanation of the algorithm employed? Thanks, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

...n(x.vector) } x1.m <- put.missings(x1, .3) x2.m <- put.missings(x2, .6) tX.m <- rbind(table(na.omit(x1.m)), table(na.omit(x2.m))) fisher.test(tX.m) With regards, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

(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"

...re and diastolic blood pressure do they substitute for each other or do they tend to occur in the same model. Thanks all (and again - thanks Thomas for all your patient answers so far) Tal ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: www.talgalili.com www.biostatistics.co.il [[alternative HTML version deleted]]

...(mpg , list(round(wt), am) , data = mtcars, las = 2) Hope this will benefit others, and also hope that others will improve on this code and give it back to the community, Tal galili -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]

..."red", cex = 2) text( .4,.1, "and how to find coordinates ?",col = "blue", cex = 1.5) text( .89,.3, "and rotation ?!",col = "white", srt = -90, cex = 2) Thanks, Tal. ---------------------------------------------- Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blog: www.talgalili.com [[alternative HTML version deleted]]

...ar.test(aa) # p-value = NA install.packages("outliers") require(outliers) cochran.test(aa) # Error in x$terms : $ operator is invalid for atomic vectors Thanks, Tal -- ---------------------------------------------- My contact information: Tal Galili Phone number: 972-50-3373767 FaceBook: Tal Galili My Blogs: http://www.r-statistics.com/ http://www.talgalili.com http://www.biostatistics.co.il [[alternative HTML version deleted]]