Displaying 2 results from an estimated 2 matches for "randomformula".
2012 May 11
0
contrasts with an imbalance in a factor
...NA,
-200L))
## this is what I really want to be able to do, but without contrast
## complaining about the imbalance in the number of rows
cat("### With scanner\n")
fixedFormula=as.formula("fmri ~ group * task + scanner")
randomFormula = as.formula("random = ~ 1 | subject")
mylme = lme(fixed=fixedFormula, random=randomFormula, data=my.model)
## now look at the risky choices (40outcome and 80outcome) versus the
## safe choices (20outcome)
con=contrast(
mylme,
a=list(group="Group1", task=c("40outcome&q...
2012 Apr 18
0
Error in eval when using contrast and nlme
...;,
"group", "brik", "scanner", "fmri"), row.names = c(NA, -210L), class = "data.frame")
runLme = function (inData, inZ, inNumberOfOutputBriks, inContrasts, inAnovaIndices, inModel, inModelFormula, inRandomFormula) {
cat("There will be " , inNumberOfOutputBriks, " stats briks\n")
outStats <- vector(mode="numeric", length=inNumberOfOutputBriks)
##cat (inData, "\n")
if ( ! all(inData == 0 ) ) {
##try(
##if( inherits(
print(inModelFormula)
print...