Hi r-users,
I have this code here, but I just wonder how do I use 'sapply' to make it more
efficient
lamda_cor <- eigen(winter_cor)$values
> lamda_cor
[1] 1.3459066 1.0368399 0.8958128 0.7214407
lamda_cxn <- function(dt)
{ n <- length(dt)
term <- vector(length=n, mode="numeric")
for (i in 1:n)
{ term[i] <- (dt[i]/n)*log(dt[i]/n) }
#sum(term)
cxn <- 1 + (1/log(n))*sum(term)
cxn
}
lamda_cxn(lamda_cor)
> lamda_cxn(lamda_cor)
[1] 0.01861457
Thank you so much for all helps given.
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
