Running 'invariance.R' [5s]
  Running 'lavaan.R' [34s]
  Running 'tree.R' [4s]
  Running 'vim.R' [8s]
Running the tests in 'tests/invariance.R' failed.
Complete output:
  > #
  > # testing invariance 
  > #
  > 
  > # simulation:
  > # factor structure holds over SES but not over age
  > # both SES and age predict a mean difference in cognitive outcome
  > 
  > set.seed(123)
  > 
  > require("semtree")
  Loading required package: semtree
  Loading required package: OpenMx
  > require("lavaan")
  Loading required package: lavaan
  This is lavaan 0.6-3
  lavaan is BETA software! Please report any bugs.
  
  Attaching package: 'lavaan'
  
  The following object is masked from 'package:OpenMx':
  
      vech
  
  > #
  > lambda <- list()
  > age <- c(0,0,1,1) # 0=young, 1 =old
  > ses <- c(0,1,0,1) # 0=low, 1=high
  > lambda[[1]] <- c(1,0.9,0.8,0.8)
  > lambda[[2]] <- c(1,0.9,0.8,0.8)
  > lambda[[3]] <- c(1,0.4,0.2,0.9)
  > lambda[[4]] <- c(1,0.4,0.2,0.9)
  > cogmean <- 80-age*30 + ses*20
  > cogsd <- 1
  > errsd <- 1
  > 
  > Nsub <- 200 # persons per agexSES group
  > 
  > cbind(age,ses,cogmean)
       age ses cogmean
  [1,]   0   0      80
  [2,]   0   1     100
  [3,]   1   0      50
  [4,]   1   1      70
  > 
  > # simulate data from 4 groups
  > data <- c()
  > for (i in 1:4) {
  +   cogsim <- rnorm(n = Nsub,mean = cogmean[i],cogsd)
  +   scores <- as.matrix(t(outer(lambda[[i]],cogsim))) + rnorm(Nsub*4,0,errsd)
  +   data <- rbind(data,scores)
  + }
  > data <- data.frame(data)
  > names(data) <- paste0("x",1:4)
  > 
  > fulldata <- cbind(data, age=rep(age,each=Nsub),ses=rep(ses,each=Nsub))
  > #
  > model<-"
  + ! regressions 
  + F=~1.0*x1
  + F=~l2*x2
  + F=~l3*x3
  + F=~l4*x4
  + ! residuals, variances and covariances
  + x1 ~~ VAR_x1*x1
  + x2 ~~ VAR_x2*x2
  + x3 ~~ VAR_x3*x3
  + x4 ~~ VAR_x4*x4
  + F ~~ 1.0*F
  + ! means
  + F~1
  + x1~0*1;
  + x2~0*1;
  + x3~0*1;
  + x4~0*1;
  + ";
  > result<-lavaan(model, data=data, fixed.x=FALSE, missing="FIML");
  Warning message:
  In lav_object_post_check(object) :
    lavaan WARNING: some estimated ov variances are negative
  > 
  > manifests<-c("x1","x2","x3","x4")
  > latents<-c("F")
  > model <- mxModel("Unnamed_Model", 
  +                  type="RAM",
  +                  manifestVars = manifests,
  +                  latentVars = latents,
  +                  mxPath(from="F",to=c("x1","x2","x3","x4"), free=c(FALSE,TRUE,TRUE,TRUE),
  +                         value=c(1.0,1.0,1.0,1.0) , arrows=1, label=c("F__x1","l2","l3","l4") ),
  +                  mxPath(from="one",to=c("F"), free=c(TRUE), value=c(1.0) , arrows=1, label=c("const__F") ),
  +                  mxPath(from="one",to=c("x2","x3","x4"), free=c(TRUE,TRUE,TRUE), value=c(1.0,1.0,1.0) , arrows=1, label=c("const__x2","const__x3","const__x4") ),
  +                  mxPath(from="x1",to=c("x1"), free=c(TRUE), value=c(1.0) , arrows=2, label=c("VAR_x1") ),
  +                  mxPath(from="x2",to=c("x2"), free=c(TRUE), value=c(1.0) , arrows=2, label=c("VAR_x2") ),
  +                  mxPath(from="x3",to=c("x3"), free=c(TRUE), value=c(1.0) , arrows=2, label=c("VAR_x3") ),
  +                  mxPath(from="x4",to=c("x4"), free=c(TRUE), value=c(1.0) , arrows=2, label=c("VAR_x4") ),
  +                  mxPath(from="F",to=c("F"), free=c(FALSE), value=c(1.0) , arrows=2, label=c("VAR_F") ),
  +                  mxPath(from="one",to=c("x1"), free=F, value=0, arrows=1),
  +                  mxData(data[1:50,], type = "raw")
  + );
  > 
  > result <- mxTryHard(model)
  
  Beginning initial fit attempt
  
   Lowest minimum so far:  611.424050364123
  
   OpenMx status code 6 not in list of acceptable status codes, (0,0)
  
  Beginning fit attempt 1 of at maximum 10 extra tries
  
   Lowest minimum so far:  611.424045697128
  
  Solution found
  
  
  Solution found!  Final fit=611.42405 (started at 196710.97)  (2 attempt(s): 2 valid, 0 errors)
  
  Error in `[[<-.data.frame`(`*tmp*`, i, value = c("l2", NA)) : 
    replacement has 2 rows, data has 1
  Calls: mxTryHard ... [<-.data.frame -> xpdrows.data.frame -> [[<- -> [[<-.data.frame
  Execution halted
Running the tests in 'tests/tree.R' failed.
Complete output:
  > set.seed(789)
  > require("semtree")
  Loading required package: semtree
  Loading required package: OpenMx
  > data(lgcm)
  > 
  > lgcm$agegroup <- as.ordered(lgcm$agegroup)
  > lgcm$training <- as.factor(lgcm$training)
  > lgcm$noise <- as.numeric(lgcm$noise)
  > 
  > # LOAD IN OPENMX MODEL.
  > # A SIMPLE LINEAR GROWTH MODEL WITH 5 TIME POINTS FROM SIMULATED DATA
  > 
  > manifests <- names(lgcm)[1:5]
  > lgcModel <- mxModel("Linear Growth Curve Model Path Specification",
  +                             type="RAM",
  +                             manifestVars=manifests,
  +                             latentVars=c("intercept","slope"),
  +                             # residual variances
  +                             mxPath(
  +                               from=manifests,
  +                               arrows=2,
  +                               free=TRUE,
  +                               values = c(1, 1, 1, 1, 1),
  +                               labels=c("residual1","residual2","residual3","residual4","residual5")
  +                             ),
  +                             # latent variances and covariance
  +                             mxPath(
  +                               from=c("intercept","slope"),
  +                               connect="unique.pairs",
  +                               arrows=2,
  +                               free=TRUE,
  +                               values=c(1, 1, 1),
  +                               labels=c("vari", "cov", "vars")
  +                             ),
  +                             # intercept loadings
  +                             mxPath(
  +                               from="intercept",
  +                               to=manifests,
  +                               arrows=1,
  +                               free=FALSE,
  +                               values=c(1, 1, 1, 1, 1)
  +                             ),
  +                             # slope loadings
  +                             mxPath(
  +                               from="slope",
  +                               to=manifests,
  +                               arrows=1,
  +                               free=FALSE,
  +                               values=c(0, 1, 2, 3, 4)
  +                             ),
  +                             # manifest means
  +                             mxPath(
  +                               from="one",
  +                               to=manifests,
  +                               arrows=1,
  +                               free=FALSE,
  +                               values=c(0, 0, 0, 0, 0)
  +                             ),
  +                             # latent means
  +                             mxPath(
  +                               from="one",
  +                               to=c("intercept", "slope"),
  +                               arrows=1,
  +                               free=TRUE,
  +                               values=c(1, 1),
  +                               labels=c("meani", "means")
  +                             ),
  +                     mxData(lgcm,type="raw")
  +                     )
  > 
  > 
  > # TREE CONTROL OPTIONS.
  > # TO OBTAIN BASIC/DEFAULT SMETREE OPTIONS, SIMPLY TPYE THE FOLLOWING:
  > 
  > controlOptions <- semtree.control(method = "naive")
  > controlOptions$alpha <- 0.05
  > 
  > # RUN TREE.
  > 
  > tree <- semtree(model=lgcModel, data=lgcm, control = controlOptions)
  Error in `[[<-.data.frame`(`*tmp*`, i, value = c("residual1", NA)) : 
    replacement has 2 rows, data has 1
  Calls: semtree ... [<-.data.frame -> xpdrows.data.frame -> [[<- -> [[<-.data.frame
  Execution halted
Running the tests in 'tests/vim.R' failed.
Complete output:
  > set.seed(789)
  > require("semtree")
  Loading required package: semtree
  Loading required package: OpenMx
  > data(lgcm)
  > 
  > lgcm$agegroup <- as.ordered(lgcm$agegroup)
  > lgcm$training <- as.factor(lgcm$training)
  > lgcm$noise <- as.numeric(lgcm$noise)
  > 
  > # LOAD IN OPENMX MODEL.
  > # A SIMPLE LINEAR GROWTH MODEL WITH 5 TIME POINTS FROM SIMULATED DATA
  > 
  > manifests <- names(lgcm)[1:5]
  > lgcModel <- mxModel("Linear Growth Curve Model Path Specification",
  +                             type="RAM",
  +                             manifestVars=manifests,
  +                             latentVars=c("intercept","slope"),
  +                             # residual variances
  +                             mxPath(
  +                               from=manifests,
  +                               arrows=2,
  +                               free=TRUE,
  +                               values = c(1, 1, 1, 1, 1),
  +                               labels=c("residual1","residual2","residual3","residual4","residual5")
  +                             ),
  +                             # latent variances and covariance
  +                             mxPath(
  +                               from=c("intercept","slope"),
  +                               connect="unique.pairs",
  +                               arrows=2,
  +                               free=TRUE,
  +                               values=c(1, 1, 1),
  +                               labels=c("vari", "cov", "vars")
  +                             ),
  +                             # intercept loadings
  +                             mxPath(
  +                               from="intercept",
  +                               to=manifests,
  +                               arrows=1,
  +                               free=FALSE,
  +                               values=c(1, 1, 1, 1, 1)
  +                             ),
  +                             # slope loadings
  +                             mxPath(
  +                               from="slope",
  +                               to=manifests,
  +                               arrows=1,
  +                               free=FALSE,
  +                               values=c(0, 1, 2, 3, 4)
  +                             ),
  +                             # manifest means
  +                             mxPath(
  +                               from="one",
  +                               to=manifests,
  +                               arrows=1,
  +                               free=FALSE,
  +                               values=c(0, 0, 0, 0, 0)
  +                             ),
  +                             # latent means
  +                             mxPath(
  +                               from="one",
  +                               to=c("intercept", "slope"),
  +                               arrows=1,
  +                               free=TRUE,
  +                               values=c(1, 1),
  +                               labels=c("meani", "means")
  +                             ),
  +                     mxData(lgcm,type="raw")
  +                     )
  > 
  > 
  > fr <- semforest(lgcModel, lgcm,control = semforest.control(num.trees = 25))
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  No traceback available 
  > 
  >   
  > vimp <- varimp(fr)
  Error in t.default(tree$params) : argument is not a matrix
  Calls: varimp ... eval.fun -> getIdParameterMapping -> cbind -> t -> t.default
  Execution halted