Re: [R] OT - use of R
A quick search of google scholar found 2556 cites to the classic paper by Ihaka et al. http://scholar.google.com/scholar?hl=enlr=cites=15992947024900415641 On 8/15/07, Andy Bunn [EMAIL PROTECTED] wrote: Hello, is there an up-to-date reference for how many people use R? I'm giving an R demo and want to cite wonderful R usage stats. How many people use it (or download it)? How often is R used in peer-reviewed pubs, etc. Is there any whiz-bang citation that says something like R is great and developed by the best minds in statistical computing. Any thoughts? -Andy __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Meta-Analysis of proportions
OpenBUGS should be something related to Bayesian statistics. You may refer to Chapter 12 of Handbook http://cran.r-project.org/doc/vignettes/HSAUR/Ch_meta_analysis.pdf It talks about meta-regression. On 6/28/07, Monica Malta [EMAIL PROTECTED] wrote: Dear colleagues, I'm conducting a meta-analysis of studies evaluating adherence of HIV-positive drug users into AIDS treatment, therefore I'm looking for some advice and syntax suggestion for running the meta-regression using proportions, not the usual OR/RR frequently used on RCT studies. Have already searched already several handbooks, R-manuals, mailing lists, professors, but... not clue at all... Does anyone have already tried this? A colleague of mine recently published a similar study on JAMA, but he used OpenBUGS - a software I'm not familiar with... If there is any tip/suggestion for a possible syntax, could someone send me? I need to finish this paper before my PhD qualify, but I'm completely stuck... So, any tip will be more than welcome...I will really appreciate it!!! Thanks in advance and congrats on the amazing mailing-list. Bests from Rio de Janeiro, Brazil. Monica Monica Malta Researcher Oswaldo Cruz Foundation - FIOCRUZ Social Science Department - DCS/ENSP Rua Leopoldo Bulhoes, 1480 - room 905 Manguinhos Rio de Janeiro - RJ 21041-210 Brazil phone +55.21.2598-2715 fax +55.21.2598-2779 [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Date and selection
Dear R experts, Suppose I have a data.frame recording the date and test results of some subjects like this: Name Date results John 01/01/1991 2 John 02/01/1991 3 John 09/0101991 4 Micheal 02/01/1991 4 Micheal 04/01/1991 5 How to select the earliest (or latest) test result from all subjects? Thank you. Regards, CH -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] How to load a big txt file
Easy solution will be split your big txt files by text editor. e.g. 5000 rows each. and then combine the dataframes together into one. On 6/7/07, ssls sddd [EMAIL PROTECTED] wrote: Dear list, I need to read a big txt file (around 130Mb; 23800 rows and 49 columns) for downstream clustering analysis. I first used Tumor - read.table(Tumor.txt,header = TRUE,sep = \t) but it took a long time and failed. However, it had no problem if I just put data of 3 columns. Is there any way which can load this big file? Thanks for any suggestions! Sincerely, Alex [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Proxy Under Mac OS X
Dear R programmers, I can only config. proxy under Mac OS X terminal and launch R under Terminal by Terminal: export http_proxy=http://un:[EMAIL PROTECTED]:port Under R: chooseCRANmirror(graphics=FALSE) update.packages() I don't know how to config this in R for Mac OS X Aqua GUI. I checked the relevant document but no specific steps for Mac OS X. Regards, Ch Chan -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] R's Spearman
Hi, You can try with cor.test(rank(y[1]),rank(y[2])) On 5/29/07, Raymond Wan [EMAIL PROTECTED] wrote: Hi all, I am trying to figure out the formula used by R's Spearman rho (using cor(method=spearman)) because I can't seem to get the same value as by calculating by hand. Perhaps I'm using cor wrong, but I don't know where. Basically, I am running these commands: y=read.table(file=tmp,header=TRUE,sep=\t) y IQ Hours 1 106 7 2 86 0 3 9720 4 11312 5 12012 6 11017 cor(y[1],y[2],method=spearman) Hours IQ 0.2319084 [it's an abbreviated example of one I took from Wikipedia]. I calculated by hand (apologies if the table looks strange when pasted into e-mail): IQHoursrank(IQ) rank(hours)diffdiff^2 110673 2 11 2 8601 1 00 3 9720 2 6-416 411312 5 3.5 1.52.25 512012 6 3.5 2.56.25 611017 4 5-11 26.5 rho=0.242857 where rho = (1 - ((6 * 26.5) / 6 * (6^2 - 1))). I kept modifying the table and realized that the difference in result comes from ties. i.e., if I remove the tie in rows 4 and 5, I get the same result from both cor and calculating by hand. Perhaps I'm handling ties wrong...does anyone know how R does it or perhaps I need to change how I'm using it? Thank you! Ray __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Simple question about function with glm
Thank you akk. I know it is not statistically sounded to check the distribution of response before glm. I will check the distribution of xmodel$residuals later on. About the program problem. It can print summary(xmodel) but not confint(xmodel) by amending my code as suggested by Bill Venables. Regards, CH On 5/7/07, Ben Bolker [EMAIL PROTECTED] wrote: Bill.Venables at csiro.au writes: Finally, I'm a bit puzzled why you use glm() when the simpler lm() would have done the job. You are fitting a linear model and do not need the extra paraphernaila that generalized linear models require. Bill Venables. Perhaps the original poster is confused about the difference between general (a la PROC GLM) and generalized (glm) linear models? The code is also a little puzzling because the same tests seem to be run whether p0.05 or not. Perhaps the code will eventually be written to log-transform the data if it fails the normality test? [ hint: ?boxcox in the MASS package might be a better way to go ] Ben Bolker __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Simple question about function with glm
Dear all, I coded a function called u.glm u.glm - function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt - lillie.test(x) if (lil.rslt$p.value 0.05) { cat(Logtrans=0, lillie=,lil.rslt$p.value,\n) xmodel-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat(Logtrans=1, lillie=,lil.rslt$p.value,\n) xmodel-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } } Basically I just want to test the response variable for normality before modeling. When I try to use this function, it can do the lillie's test but failed to do the glm. What's wrong with my code? Regards, CH -- The scientists of today think deeply instead of clearly. One must be sane to think clearly, but one can think deeply and be quite insane. Nikola Tesla http://www.macgrass.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.