An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110714/2614adb0/attachment.pl>
Repating a loop of lm function with different columns of database
2 messages · Jon Toledo, Weidong Gu
If I understood your question x<-data.frame(matrix(rnorm(2000,10,10),ncol=50)) sapply(1:5,function(i) summary(lm(x[,i]~x[,i+10]+x[,50]))) Weidong Gu
On Thu, Jul 14, 2011 at 2:27 PM, Jon Toledo <tintin_jb at hotmail.com> wrote:
Hi,
First let me thank you for the incredible help and resource that this forum is.
I am trying to compare the ?repeated measurement of more than 100 analytes that have been take in 70 subjects at 2 time points adjusted for the time difference of sample times(TimeDifferenceDays), therefore I wanted to do it with a function that allows me to do all at once. (131 is the column difference that separates the two different measurements of the same anlyte)
I have this one:
for(i in 1:125){return(summary(lm(Data[,i] ~ Data[,(i+131)] + Data$TimeDifferenceDays)))}
But it only gives me one result
I also wanted to get the p-value in a dataframe with three columns:
Fitst column: analyte name (that? the name of the column)Second column: pvalue of the first measure of the analte predicting the second measureThird column: The effect of time
I copy a samller example:
txt <- "V1a ?V2a ?V3a ?V1b ?V2b ?V3b ? TimeDifferenceDays2.42 ?72.4 3.75 2.46 55.4 4.44 6081.66 89.7 2.54 2.17 94.0 2.15 4192.45 112. 0.46 2.40 129.0 .42 7142.58 55.6 5.05 2.44 135.0 5.39 7212.61 332.0 22.6 3.55 238.0 16.4 729"
Data <- read.table(textConnection(txt), header = TRUE)
Thank you
Jon Toledo, MD
Postdoctoral fellow
University of Pennsylvania School of Medicine
Center for Neurodegenerative Disease Research
? ? ? ?[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org 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.