search for: ideffect

Displaying 2 results from an estimated 2 matches for "ideffect".

2012 Sep 06
0
lme( y ~ ns(x, df=splineDF)) error
...ed from http://jacobwegelin.net/tmp/r-help/. wrapper <-function(WhichApproach=1, splineDF=4 ){ # Create toy dataset nID<-25 IDnames<-LETTERS[1:nID] longdat<-data.frame( x=rep(-10:10, nID) / 3 , ID= rep( IDnames , each=21) ) set.seed(5) longdat$x<-longdat$x + rnorm(nrow(longdat))/10 IDeffect<-rnorm(nID) * 20 names(IDeffect) <- IDnames longdat$IDeffect<-IDeffect[as.character(longdat$ID)] longdat$y<- (longdat$x + longdat$x^3 + (longdat$x-1)^4 / 5 + 1/(abs(longdat$x/50) + 0.02) + longdat$IDeffect + rnorm(1:nrow(longdat)) * 2 ) longdat<-longdat[order...
2012 Sep 26
0
lme(y ~ ns(x, df=splineDF)) error
...2.1 stringr_0.6 [16] tools_2.15.1 wrapper <-function(WhichApproach=1, splineDF=4 ){ # Create toy dataset nID<-25 IDnames<-LETTERS[1:nID] longdat<-data.frame( x=rep(-10:10, nID) / 3 , ID= rep( IDnames , each=21) ) set.seed(5) longdat$x<-longdat$x + rnorm(nrow(longdat))/10 IDeffect<-rnorm(nID) * 20 names(IDeffect) <- IDnames longdat$IDeffect<-IDeffect[as.character(longdat$ID)] longdat$y<- (longdat$x + longdat$x^3 + (longdat$x-1)^4 / 5 + 1/(abs(longdat$x/50) + 0.02) + longdat$IDeffect + rnorm(1:nrow(longdat)) * 2 ) longdat<-longdat...