[R] Problem with lapply

arun smartpink111 at yahoo.com
Tue Oct 15 21:56:35 CEST 2013


Hi,

No problem.
Regarding your code:

#changed name from 'dataset' to "Testcor"
#Your code

  firm.names=colnames(Testcor)[2:length(colnames(Testcor))]
 firm.names
# [1] "PAXN"   "PED"    "PEDP"   "PM"     "PRFN"   "PRWN"   "ZG"     "ZUBN"  
 #[9] "ZURN"   "ZWM"    "Market" #####includes "Market"

#If you look my code:
firm.names <- colnames(Testcor)[!grepl("DATE|Market",colnames(Testcor))]
 firm.names
# [1] "PAXN" "PED"  "PEDP" "PM"   "PRFN" "PRWN" "ZG"   "ZUBN" "ZURN" "ZWM" 

If you wanted to use the "formula" method:
res1 <- sapply(firm.names,function(x) { cor.results <- cor.test(as.formula(paste("~", paste(x,"Market", sep="+")) ),data=Testcor,na.action=na.exclude); cor.results$estimate})
res1
#   PAXN.cor     PED.cor    PEDP.cor      PM.cor    PRFN.cor    PRWN.cor 
#-0.31113122 -0.09359550 -0.17056943  0.40025112  0.34385888  0.17935045 
 #    ZG.cor    ZUBN.cor    ZURN.cor     ZWM.cor 
 #0.20762797  0.39238270  0.74336028 -0.09166795 
 identical(res,res1)
#[1] TRUE

#or just
res2 <- sapply(firm.names,function(x) { cor.test(~get(x) + Market,data=Testcor,na.action=na.exclude)$estimate })
 identical(res,res2)
#[1] TRUE



A.K.






On Tuesday, October 15, 2013 3:15 PM, rissa <rissa_r at gmx.ch> wrote:
Thank you! Worked perfectly!

By the way: What was wrong with my original code?



--
View this message in context: http://r.789695.n4.nabble.com/Problem-with-lapply-tp4678290p4678305.html

Sent from the R help mailing list archive at Nabble.com.

______________________________________________
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.




More information about the R-help mailing list