Displaying 2 results from an estimated 2 matches for "treatvar".
Did you mean:
tratar
2007 Jan 23
1
How to evaluate an lm() object for generating warning statement in a function
...,2], length) #number
of cases per treatment level
nisq<-ni^2
no<-(1/(length(ni)-1))*(sum(ni)-(sum(nisq)/sum(ni)))
#required to calculate variance components
s2a<-((MStreat-MSErr)/no)
stot<-s2a + MSErr
treatvar<-s2a/stot*100 #calculate variance components as
a percentage of the total
errorvar<-MSErr/stot*100
list(treat.var.comp=s2a,
err.var.comp=MSErr,
p.var.treat=treatvar, p.var.err=errorvar)
}
comp.var.estimates(object1)
I'...
2007 Jan 22
0
Advice on an appropriate condition for warning in a function
...t;-anovmod[1,3]; MSErr<-anovmod[2,3] #extract Mean Squares
dataframe<- as.data.frame(object.lm[13])
ni <- tapply(dataframe[,1], dataframe[,2], length)
nisq<-ni^2
no<-(1/(length(ni)-1))*(sum(ni)-(sum(nisq)/sum(ni)))
s2a<-((MSTreat-MSErr)/no)
stot<-s2a + MSErr
treatvar<-s2a/stot*100
errorvar<-MSErr/stot*100
list(treat.var.comp=s2a, err.var.comp=MSErr,
p.var.treat=treatvar, p.var.err=errorvar)
}
comp.var.estimates(object1)
I'd like to include a "warning" statement in the function that
returns 'function requires arguments that...