Displaying 1 result from an estimated 1 matches for "limconf".
Did you mean:
libconf
2004 Oct 26
1
indexing within the function "aggregate"
...ues showed in the graph
above.
Since "aggregate" can manage functions that
return a single value (V&R, pag.35),
I wrote the following
My.fun.conf <-
function(y,x){
dieffe <- anova(lm(y~x))$Df[1]
enne <- dieffe+1
esse <- sqrt(anova(lm(y~x))$Mean[1])
limconf.95 <- qt(.975,df=dieffe)*esse/sqrt(enne)
limconf.95}
that can be used like
tag <- My.df$SEQUENCE == "1" & My.df$STD == "INT"
My.fun.conf(My.df$PPM[tag],My.df$VIAL[tag])
but it do not works in "aggregate"
since the indexing of values is no longer done
o...