t1waysub <-
function(tm,sqse,hval){
#
# Used by t1waybt to compute Welch test statistic based on trimmed means
# and squared standard errors stored in tm and sqse
#
w<-1/sqse
uval<-sum(w)
xtil<-sum(w*tm)/uval
A<-sum(w*(tm-xtil)^2)/(length(tm)-1)
B<-sum((1-w/uval)^2/(hval-1))
t1waysub<-A/(B+1)
t1waysub
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.