R语言脚本函数

加权方差函数

myfunction<-function(x,w){
  xbar<-weighted.mean(xi)
  v<-sum((x-xbar)^2*w)/sum(w)
  v}

t检验函数

myfunction<-function(x,y){
  xbar<-mean(x)
  ybar<-mean(y)
  n1<-lengths(x)
  n2<-lengths(y)
  s1<var(x)
  s2<-var(y)
  s=((n1-1)*s1+(n2-1)*s2)/(n1+n2-2)
  t=(xbar-ybar)/sqrt(s*(1/n1)+s*(1/n2))
  return(t)
}