关注今日：9 | 主题：142451
微信扫一扫
扫一扫，下载丁香园 App
即送15丁当
【求助】如何用R语言实现 二项分布 两大样本率 U检验 啊？命令怎么写啊？
页码直达：
这个帖子发布于4年零72天前，其中的信息可能已发生改变或有所发展。
问题已关闭悬赏丁当:2
如何用R语言实现 二项分布
两大样本率 U检验 啊？命令怎么写啊？例如抽查男性120人和女性100人，发现男性36人患乙肝，女性22人患乙肝，是比较性别患病率是否不同。上面只是列出的一个例子，我做的人数是上万人的。谢谢！命令怎么写啊？急着写论文，帮忙回答一下哈，实在是感谢了！
不知道邀请谁？试试他们
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
如何用R语言实现 二项分布 两大样本率 U检验 啊？命令怎么写啊？例如抽查男性120人和女性100人，发现男性36人患乙肝，女性22人患乙肝，是比较性别患病率是否不同。上面只是列出的一个例子，我做的人数是上万人的。谢谢！命令怎么写啊？急着写论文，帮忙回答一下哈，实在是感谢了！请高手解答！！！！！！！！！！！！
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
不明白为什么 U检验？
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
两大样本率 检验都做u或Z检验的吧,例如中华流行病的一篇文章里写：两大样本率的比较采用二项分布Z检验。p≤0．05为差异有统计学意义。何用R语言实现 二项分布 两大样本率 U检验 啊？命令怎么写啊？例如抽查男性120人和女性100人，发现男性36人患乙肝，女性22人患乙肝，是比较性别患病率是否不同。上面只是列出的一个例子，我做的人数是上万人的。谢谢！命令怎么写啊？急着写论文，帮忙回答一下哈，实在是感谢了！
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
卡方检验不可以吗，U检验应该是u检验吧？
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
卡方检验不可以的，在中华流行病发的文章里都 要 两大样本率的比较采用二项分布Z检验的
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园版主
jing123fu815 卡方检验不可以的，在中华流行病发的文章里都 要 两大样本率的比较采用二项分布Z检验的你确定？不要断章取义
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
jing123fu815 卡方检验不可以的，在中华流行病发的文章里都 要 两大样本率的比较采用二项分布Z检验的1.当自由度为1时，卡方检验等同于u检验。2.两个样本率差别的统计学检验，可以根据公式手工计算。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
jing123fu815 如何用R语言实现 二项分布
两大样本率 U检验 啊？命令怎么写啊？例如抽查男性120人和女性100人，发现男性36人患乙肝，女性22人患乙肝，是比较性别患病率是否不同。上面只是列出的一个例子，我做的人数是上万人的。谢谢！命令怎么写啊？急着写论文，帮忙回答一下哈，实在是感谢了！R code如下：prop.test(x=c(36,22), n=c(120,100),alternative=&two.sided&,conf.level=.95)主要结果如下：data:
c(36, 22) out of c(120, 100) X-squared = 1.4097, df = 1, p-value = 0.2351alternative hypothesis: two.sided 95 percent confidence interval: -0.. sample estimates:prop 1 prop 2
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
以下是教材上SAS的计算方法，跟上面的计算结果不一样啊？？？？非常谢谢大家的解答！data prg11_6;n1=120;n2=110;x1=36;x2=22;p1=x1/n1;p2=x2/n2;pc=(x1+x2)/(n1+n2);sp=sqrt(pc*(1-pc)*(1/n1+1/n2));u=(p1-p2)/p=(1-probnorm(abs(u)))*2;format u p 8.4;
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
十万火急！请各位高手解答一下，结果里面的Q，I,P等代表的是什么啊？？写论文看哪些指标？如下一篇文章的表述对吗？？？R或者MetaAnalyst能检验模型吗?q值是用来检验模型有没有统计学意义的吗？？？？？2.2.1同质性检验及总体患病率情况的分析对纳入研究的10个调查资料采用随机效用模型分析，模型有统计学意义（q=0.996，P=0.000）；2.2.2老年期痴呆各个亚组患病率情况的分析对不同性别不同年龄段、文化程度和婚姻状况进行分析。在模型检验的过程中，除60～64岁人群的资料模型没有统计学意义之外（q=0.779，P=0.24﹥0.05），其余均有统计学意义；
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
首先，以上数据不是用U检验，而应该是卡方检验来完成的。除此之外，可以用OR,RR，AR等来表示性别对疾病的影响。R CODE如下：如果想获得OR,RR,AR及95%可信区间，需要安装一个PACKAGE (epiR)，接下来就方便很多了。install.packages(&epiR&)library(epiR)# male diseaseda=36# male non-diseasedb=120-36# female diseasedc=22# female not diseasedd=100-22# put 4 numbers in 2 by 2 table and saved as datamatrix(c(a,b,c,d),nrow=2,ncol=2,byrow=T)-&data## Chisq test if you need p-valuechisq.test(data)## 2by2 table, if you need OR, RR, AR and conf. intervalepi.2by2(data,method=&cohort.count&)我只会R，SAS方面不能给你意见。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
jeffylin edited on
丁香园准中级站友
jing123fu815 以下是教材上SAS的计算方法，跟上面的计算结果不一样啊？？？？非常谢谢大家的解答！data prg11_6;n1=120;n2=110;x1=36;x2=22;p1=x1/n1;p2=x2/n2;pc=(x1+x2)/(n1+n2);sp=sqrt(pc*(1-pc)*(1/n1+1/n2));u=(p1-p2)/p=(1-probnorm(abs(u)))*2;format u p 8.4;你这个公式就是大样本率的U检验公式，是正确的。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园准中级站友
zhangdog 1.当自由度为1时，卡方检验等同于u检验。一个是卡方分布，一个是正态分布，一个计算的是卡方值（见excel），一个计算的是U值(见SAS结果)，不过算出来的P值是一样的。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
mike305 edited on
丁香园准中级站友
SAS的结果和EXCEL卡方的结果不太一样。理论频数都很大，R软件的卡方搞不懂为什么要校正
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
以下摘自《中国医学百科全书－医学统计学》（可能知识比较老，仅供参考）：1.两个样本率作比较的目的在于推断两个样本各自代表的总体率是否相等 。常 用的X2检验、u检验及对数似检验，也可用平方根纸作检验。当样本含量甚小时宜用四格表的确切概率法。2.当自由度为1时，u2(alpha)=x2(alpha)，故u检验与x2检验等价。用sas教材上的u检验code的计算结果与卡方验的结果基本上一致（楼主有一个数据输错了，将100错输成110，所以差距比较大），我们以女性总数100计算，结果差不多。说明百科全书说的有一定的道理。1）SAS u检验结果（用楼主的code）,结果如下：
0.17992）用stata卡方检验（命令：tabi 36 84\22 78,all），得主要结果
Pearson chi2(1) =
Pr = 0.180 likelihood-ratio chi2(1) =
Pr = 0.178
Cram閞's V =
ASE = 0.150
Kendall's tau-b =
ASE = 0.066
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
zhangdog edited on
丁香园荣誉版主
jeffylin 首先，以上数据不是用U检验，而应该是卡方检验来完成的。除此之外，可以用OR,RR，AR等来表示性别对疾病的影响。u检验也是可以用的。R软件中的prop.test（）函数用于两个或多个样本率的比较。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
zhangdog edited on
丁香园荣誉版主
mike305 SAS的结果和EXCEL卡方的结果不太一样。理论频数都很大，R软件的卡方搞不懂为什么要校正数据不正确，所以结果不太一样
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园准中级站友
zhangdog 数据不正确，所以结果不太一样J多谢张老师，之前是数据录错了，刚更正了，用那位仁兄的R程序，算出来的是校正的卡方，而按照理论频数来看是不要校正的。我也买了徐天和写的中华医学统计百科全书，很不错。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
mike305 多谢张老师，之前是数据录错了，刚更正了，用那位仁兄的R程序，算出来的是校正的卡方，而按照理论频数来看是不要校正的。我也买了徐天和写的中华医学统计百科全书，很不错。呵呵，其实R中prop.test（）函数专门用于样本率的比较的，我看到的国外统计网站用它啊。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
zhangdog 数据不正确，所以结果不太一样不好意思教材上的例子是n2=110，不是100，开始的时候打错了，非常谢谢大家的解答！以下是教材上SAS的计算方法
。data prg11_6;n1=120;n2=110;x1=36;x2=22;p1=x1/n1;p2=x2/n2;pc=(x1+x2)/(n1+n2);sp=sqrt(pc*(1-pc)*(1/n1+1/n2));u=(p1-p2)/p=(1-probnorm(abs))*2;format u p 8.4;
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园准中级站友
zhangdog 呵呵，其实R中prop.test（）函数专门用于样本率的比较的，我看到的国外统计网站用它啊。多谢，我用PROP.TEST()函数做了一下，出来结果了，是卡方检验的结果，但P值和U检验的P值是一样的，说明关于2个率的检验卡方检验和U检验是一致的，能得到一样的P。在R中做卡方检验时，校正与不校正可通过correct = FALSE或TURE,进行指定。如上面那个R程序的结果也是可以用不校正的结果的。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
不过把例子是n2=100换成110，用R（老师提供的命令）计算的结果怎么跟sas计算的结果不一样呢。用R计算结果：& prop.test(x=c(36,22), n=c(120,110),alternative=&two.sided&,conf.level=.95)
2-sample test for equality of proportions with continuity
correctiondata:
c(36, 22) out of c(120, 110) X-squared = 2.5361, df = 1, p-value = 0.1113alternative hypothesis: two.sided 95 percent confidence interval: -0.. sample estimates:prop 1 prop 2
0.2用SAS计算结果:u=1.7745
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
jing123fu815 不过把例子是n2=100换成110，用R（老师提供的命令）计算的结果怎么跟sas计算的结果不一样呢。用R计算结果：& prop.test(x=c(36,22), n=c(120,110),alternative=&two.sided&,conf.level=.95)
2-sample test for equality of proportions with continuity
correctiondata:
c(36, 22) out of c(120, 110) X-squared = 2.5361, df = 1, p-value = 0.1113alternative hypothesis: two.sided 95 percent confidence interval: -0.. sample estimates:prop 1 prop 2
0.2用SAS计算结果:u=1.7745
p=0.0811已做解释，可从。具体可以参考：以您的数据为例，我的结果是采取连续校正的，如何采用不校正，具体code为：prop.test(x=c(36,22), n=c(120,110),alternative=&two.sided&,conf.level=.95,correct=FALSE)结果如下：data:
c(36, 22) out of c(120, 110) X-squared = 3.0433, df = 1, p-value = 0.08107alternative hypothesis: two.sided 95 percent confidence interval: -0.. sample estimates:prop 1 prop 2
0.2 虽然u值与卡方值不同，但P值相同，与SAS结果完全一致。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
zhangdog edited on
谢谢大家的解答，用R进行U检验，在论文里应该怎么表达呢，是说采用R里的prop.test进行的卡方检验？还是说进行U检验？困惑！
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
另外，可以将SAS的code稍微修改一下，可以在R上运行，具体如下：n1=120n2=110x1=36x2=22p1=x1/n1p2=x2/n2pc=(x1+x2)/(n1+n2)sp=sqrt(pc*(1-pc)*(1/n1+1/n2))u=(p1-p2)/spp=(1-pnorm(abs(u)))*2p得P值为0.。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
zhangdog edited on
还有就是有没有人会在??R软件或者SAS里计算失安全系数啊？命令如何写啊？请高手解答！！！！谢谢！！！四、失安全系数( fail-safe number. Nfs)采用失安全系数来估计发表偏倚的程度,分别计算各因素的失安全系数(详见表1)。失安全系数反映Meta分析结果的稳定性,其值越大表明结果的稳定性越好,结论被推翻的可能性越小[3]。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
jing123fu815 还有就是有没有人会在??R软件或者SAS里计算失安全系数啊？命令如何写啊？请高手解答！！！！谢谢！！！四、失安全系数( fail-safe number. Nfs)采用失安全系数来估计发表偏倚的程度,分别计算各因素的失安全系数(详见表1)。失安全系数反映Meta分析结果的稳定性,其值越大表明结果的稳定性越好,结论被推翻的可能性越小[3]。建议你先把一些基本知识弄懂了再做吧。这样做起很吃力的。建议你买《实用循证医学方法学》
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
我已经看了一些日子，断断续续大概半年了，由于是第一次做，不太清楚失安全系数而已，请高手帮忙解答！失安全系数=(s/1.645)^2-k,
s为每个纳入的研究的p值的Z值之和，率的meta分析纳入的研究提供的P值到底是哪个啊？是怎么算出来的啊？不知中华流行病《中国内地居民年被动吸烟率Meta分析》里是怎么算出来的啊？
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
丁香园荣誉版主
jing123fu815 我已经看了一些日子，断断续续大概半年了，由于是第一次做，不太清楚失安全系数而已，请高手帮忙解答！失安全系数=(s/1.645)^2-k,
s为每个纳入的研究的p值的Z值之和，率的meta分析纳入的研究提供的P值到底是哪个啊？是怎么算出来的啊？不知中华流行病《中国内地居民年被动吸烟率Meta分析》里是怎么算出来的啊？在丁香园中《实用循证医学方法学》样张试读中，有关于R实现失安全数的具体方法，可以免费下载学习。
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
谢谢！刚才下来看了，不过他是用RR值进行计算的，但率的meta不知怎么弄的啊？请高手解答一下？不知中华流行病《中国内地居民年被动吸烟率Meta分析》里是怎么算出来的啊？??
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
1.prop.test() 对于binomial experiment分2组的情形应该足够了2. fail-safe N不推荐用，这个很多文章论述过了。你要做用metafor吧library(metafor)fsn(yi, vi, data, type=&Rosenthal& )
微信扫一扫
广告宣传推广
政治敏感、违法虚假信息
恶意灌水、重复发帖
违规侵权、站友争执
附件异常、链接失效
biooo007 edited on
关于丁香园解决java 查看JDK中底层源码的实现方法_Java_ThinkSAAS
解决java 查看JDK中底层源码的实现方法
解决java 查看JDK中底层源码的实现方法
1.点 “window”-& "Preferences" -& "Java" -& "Installed JRES"2.此时"Installed JRES"右边是列表窗格，列出了系统中的 JRE 环境，选择你的JRE，然后点边上的 "Edit..."， 会出现一个窗口(Edit JRE)3.选中rt.jar文件的这一项：“c:program filesjavajre_1.5.0_06librt.jar”点 左边的“+” 号展开它，4.展开后，可以看到“Source Attachment:(none)”，点这一项，点右边的按钮“Source Attachment...”, 选择你的JDK目录下的 “src.zip”文件5.一路点"ok",结束。
PHP开发框架
开发工具/编程工具
服务器环境
ThinkSAAS商业授权：
ThinkSAAS为用户提供有偿个性定制开发服务
ThinkSAAS将为商业授权用户提供二次开发指导和技术支持
让ThinkSAAS更好，把建议拿来。
开发客服微信}