多重假设检验<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />大家好,
我很快就需要进行大量的假设检验。因为众所周知,在这种情况下依赖p值是错误的,由于I类错误,我想使用错误发现率,更准确地说是积极的错误发现率(Storey),并且计算q值作为我的“调整后的p值”。
我的问题分为两个:
1。是否有任何代码来计算JMP 9中的pFDR和相应的q值?我在互联网上看到了一些东西,我认为这些功能已经存在于JMP Clinical中。我只有JMP 9的许可证,我不知道JMP的临床费用是多少,或者我是否有任何预算。因此,如果某处有一个简单的JSL代码,它对我来说非常有用。
2。如果存在这样的代码,则很可能需要p值向量作为输出。考虑到我有一个带有1个分类变量和500个连续变量的JMP文件,我想为500中的每一个运行一个t检验,并在一个新文件中保存一个p值的变量,包含其中500个。我如何在JSL中做到这一点?有可能吗?JMP将自动运行500次测试并仅保存p值?
先感谢您....
多重假设检验感谢大家的帮助。
周末过后,我将尝试使用代码并阅读幻灯片。
感谢Dave的链接,它看起来就像我需要的那样,谢谢大家的帮助,我很感激!
多重假设检验这些幻灯片中有很多好的信息。谢谢你分享戴夫。
多重假设检验以下是一些解释如何操作的幻灯片使用JSL连接到R
多重假设检验对于矩阵||是一个水平连接运算符。Concat(m,owr ......)是等价的。
我没有FDR的经验,但我在JMP-R情感方面做了一些尝试。以下是一个例子。前半部分与上面的例子相同。输出显示在jmp日志中。要将输出作为字符串输入,请将R Submit放入Log Capture()函数中。
DT = 打开 ( “$ ENGLISH_SAMPLE_DATA / Big Class.jmp” ) ;
米 = [] ;
对于 ( 一世 = 4 ,一世 <= N Col ( DT ) ,一世 ++,
流 = 单程 ( ÿ ( 柱 ( 一世 ) ) , X ( : 性别 ) , t测试 ( 1 ) ,无形 ) ;
OWR = 报告 ( 流 ) ;
米 = 米 || OWR [ “t测试” ] [ 数字Col Box ( 4 )] [ 1 ] ;
) ;
DTP = 新表 ( “t检验” ,新栏目 ( “P” ,数字 ,设定值 ( 米 ) ) ) ;
VARNAME = “P” ;
R初始 () ; //初始化R连接
R发送 ( DTP ) ; //将数据表发送给R.
//使用R命令提交字符串
R提交 ( “数据< - dtp
名称(数据)
连接(数据)
意思(” || VARNAME || “)” ) ;
R项 () ;
多重假设检验谢谢,这是一个好主意......
还有一件事我需要完成它(接受我的道歉,因为我太忙了,我没有时间真正学习JSL ......)。从你的代码中我得到一个包含1列p值的JMP表,对吧?我现在应该编写一个连接到R的JSL代码,如果是这样,我需要一点指导,或者,如果我导出表并直接从R工作,它应该可以工作,不是吗?
如果有人能给我写一个快速简单的例子,说明如何在重新启动新表之后,连接到R并运行像mean(p)这样简单的东西,我可以稍后为pFDR修改它。
和最后一个问题(承诺),在行:
米 = 米 || OWR [ “t测试” ] [ 数字Col Box ( 4 )] [ 1 ] ;
为什么要么 owr [...]?我的意思是,为什么||需要吗?
万分感谢,你节省了我很多时间将大数据集导出到R并在那里编码....
多重假设检验1。我不知道可以自由使用的JMP代码,但R中有几个选项。 http://strimmerlab.org/notes/fdr.html 列出了一些与FDR相关的包。JMP 9可以连接到R并从中检索结果。这可能比从头开始编写JSL中的算法要简单得多。
2。那应该是可能的。下面是一个可能会给你一些想法的例子。它通过以单向拟合循环遍历固定x列右侧的所有数字y列并从报告中检索p值(双侧t检验)并将其插入矩阵来工作。矩阵最终被添加到新数据表中的列中。
DT = 打开 ( “$ ENGLISH_SAMPLE_DATA / Big Class.jmp” ) ;
米 = [] ;
对于 ( 一世 = 4 ,一世 <= N Col ( DT ) ,一世 ++,
流 = 单程 ( ÿ ( 柱 ( 一世 ) ) , X ( : 性别 ) , t测试 ( 1 ) ,无形 ) ;
OWR = 报告 ( 流 ) ;
米 = 米 || OWR [ “t测试” ] [ 数字Col Box ( 4 )] [ 1 ] ;
) ;
新表 ( “t检验” ,新栏目 ( “P” ,数字 ,设定值 ( 米 ) ) ) ;
多重假设检验我认为pFDR存在于JMP Genomics中,但不存在于标准(或Pro)JMP中