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