Tag: chi squared

从JAVA调用R以获得卡方统计和p值

我在JAVA中有两个4 * 4矩阵，其中一个矩阵保存观察计数和其他预期计数。 我需要一种自动的方法来计算这两个矩阵之间的卡方统计量的p值; 但是，就我所知，JAVA没有这样的function。 我可以通过将两个矩阵读成R作为.csv文件格式来计算卡方和它的p值，然后使用chisq.test函数，如下所示： obs<-read.csv("obs.csv") exp<-read.csv("exp.csv") chisq.test(obs,exp) 其中.csv文件的格式如下： A, C, G, T A, 197.136, 124.32, 63.492, 59.052 C, 124.32, 78.4, 40.04, 37.24 G, 63.492, 40.04, 20.449, 19.019 T, 59.052, 37.24, 19.019, 17.689 给定这些命令，R将给出格式的输出： X-squared = 20.6236, df = 9, p-value = 0.01443 其中包括我正在寻找的p值。 有谁知道自动化过程的有效方法： 1）将我的矩阵从JAVA输出到.csv文件中2）将.csv文件上传到R 3）将.csv文件上的chisq.test调用到R中4）将输出的p值返回到JAVA？ 谢谢你的帮助….