[R] Problem with lapply
arun
smartpink111 at yahoo.com
Tue Oct 15 20:32:48 CEST 2013
Hi,
Try:
Testcor <- read.table("Testcor.txt",sep="",header=TRUE,stringsAsFactors=FALSE)
Testcor$DATE <- as.Date(Testcor$DATE)
firm.names <- colnames(Testcor)[!grepl("DATE|Market",colnames(Testcor))]
res <- sapply(firm.names, function(x) {cor.results <- cor.test(Testcor[,x],Testcor[,"Market"],na.action=na.exclude); cor.results$estimate})
# 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
A.K.
On Tuesday, October 15, 2013 1:57 PM, rissa <rissa_r at gmx.ch> wrote:
Hi together
I'm pretty new to R, so excuse me if it is a basic question.
I have a big dataset (extract of it found in the attachment) of returns from
firms. I'd like to compute the Pearson correlation of each firm with the
"Market" and the corresponding p-Value. So I thought of making a list of
'cor.test's and then extract the needed values with a for loop. What I did
so far...
dataset=Testcor
dataset$DATE<-as.Date(dataset$DATE,"%Y-%m-%d")
firm.names=colnames(dataset)[2:length(colnames(dataset))]
cor.results1=lapply(firm.names,function(x)cor.test(as.formula(paste(firm.names[firm.names==x],Market)),data=dataset,na.action=na.exclude))
But I always get error messages like
Error in paste(firm.names[firm.names == x], Market) :
object 'Market' not found
or (when I attach dataset)
Error in parse(text = x) : <text>:1:6: unexpected numeric constant
1: PAXN NA
(Next step would have been something like
for(i in 1:ncol(cor.results1)){
newmatrix[1,i]= cor.results1 [[i]]$estimate}
etc)
Could you please help me?
Thank you so much!
Rissa
Testcor.txt <http://r.789695.n4.nabble.com/file/n4678290/Testcor.txt>
--
View this message in context: http://r.789695.n4.nabble.com/Problem-with-lapply-tp4678290.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