CalculateMatrixofMeansAndStandardErrors: CalculateMatrixofMeansAndStandardErrors

In lsmeans: Calculates the lsmeans and their variances for lm, gls,lme,mer,lmer and glmer models

Description

Returns the means and their standard errors for each factor or factors combination included in formula argument of the call to lm, gls, mle,lmer or glmer functions.

Usage

CalculateMatrixofMeansAndStandardErrors(myModel, myData,myModelTerm, Mmatrix,MissingCells)

Arguments

myModel	An object of class lm, gls, mle, glm, mer, lmerMod or glmerMod
myData	A data frame included in the data argument of the call to the lm, gls, mle, glm, lmer and glmer functions.
myModelTerm	myModelTerm is a string identifying the term of the model for which the means and estandar errors are requested. These terms must not contain references to covariables. Only the factors involved in the model formula can be included. Let suppose that the model formula is Y~FactorA+FactorB+FactorA:FactoB +CovariableX, then admissible calls to CalculateMatrixofMeansAndStandardErrors are: CalculateMatrixofMeansAndStandardErrors(myModel, myData,"FactorA",Mmatrix,MissingCells) CalculateMatrixofMeansAndStandardErrors(myModel, myData,"FactorB",Mmatrix,MissingCells) CalculateMatrixofMeansAndStandardErrors(myModel, myData,"FactorA:FactorB",Mmatrix,MissingCells)
Mmatrix	A matrix returned by CalculateMmatrix(myModel,myData) function.
MissingCells	MissingCells=FALSE. If TRUE it will return, when possible, estimations for means for levels combinations not present in the data set.

Value

A two columns matrix containing means and their estandard error for the levels of the model term identified by myModelTerm.

Author(s)

Julio A. Di Rienzo <email: [e-mail]>

See Also

MeanAndVariance

Examples

library(nlme)
library(lme4)
library(lsmeans)
Oats[,"nitro"]=as.factor(Oats[,"nitro"])

myModel<-lme(yield~1+Variety+nitro+Variety:nitro,random =~1|Block/Variety,,method="REML",data=Oats)
Mmatrix=CalculateMmatrix(myModel,Oats)
CalculateMatrixofMeansAndStandardErrors(myModel,Oats,"Variety",Mmatrix,MissingCells=FALSE)
CalculateMatrixofMeansAndStandardErrors(myModel,Oats,"nitro",Mmatrix,MissingCells=FALSE)
CalculateMatrixofMeansAndStandardErrors(myModel,Oats,"Variety:nitro",Mmatrix,MissingCells=FALSE)

lsmeans documentation built on May 31, 2017, 4:45 a.m.