是否有一些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 (将#修改为@）