如何解决是否有一些R软件包用于生成偏移的卡方分布?
我想从左方平方f(Xa)生成随机变量。我试过rchisq(10,1)-a,但这是错误的,它等于f(X)-a而不是f(Xa )。是否存在一些从f(Xa)生成随机数的解决方案,f =是卡方分布?
desired_length<-1000
empty_list <- vector(mode = "list",length = desired_length)
empty_list1 <- vector(mode = "list",length = desired_length)
empty_list2<-vector(mode="list",length=desired_length)
empty_list3<-vector(mode="list",length=desired_length)
empty_list4<-vector(mode="list",length=desired_length)
for (i in 1:1000) {
set.seed(i)
h<-rchisq(10,3)
set.seed(i)
g<-rchisq(10,3)
set.seed(i)
g1<-rchisq(10,3)-0.5
set.seed(i)
g2<-rchisq(10,3)-1
set.seed(i)
g3<-rchisq(10,3)-1.5
set.seed(i)
g4<-rchisq(10,3)-2
a<-cbind(h,g)
b<-cbind(h,g1)
c<-cbind(h,g2)
d<-cbind(h,g3)
e<-cbind(h,g4)
empty_list[[i]]<-a
empty_list1[[i]]<-b
empty_list2[[i]]<-c
empty_list3[[i]]<-d
empty_list4[[i]]<-e
}
pvalue<-numeric(1000)
pvalue1<-numeric(1000)
pvalue2<-numeric(1000)
pvalue3<-numeric(1000)
pvalue4<-numeric(1000)
x<-numeric(5)
for (i in 1:1000){
pvalue[i]<-t.test(empty_list[[i]][,1],empty_list[[i]][,2])$p.value
pvalue1[i]<-t.test(empty_list1[[i]][,empty_list1[[i]][,2])$p.value
pvalue2[i]<-t.test(empty_list2[[i]][,empty_list2[[i]][,2])$p.value
pvalue3[i]<-t.test(empty_list3[[i]][,empty_list3[[i]][,2])$p.value
pvalue4[i]<-t.test(empty_list4[[i]][,empty_list4[[i]][,2])$p.value
}
x[1]<-sum(pvalue<0.05)/1000
x[2]<-sum(pvalue1<0.05)/1000
x[3]<-sum(pvalue2<0.05)/1000
x[4]<-sum(pvalue3<0.05)/1000
x[5]<-sum(pvalue4<0.05)/1000
location<-seq(0,2,by =0.5)
plot(location,x,ylab="Power for chi3^2 distribution",xlab="location difference",type = "l",col="blue")
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。