首页 文章
  • 0 votes
     answers
     views

    具有分类预测因子和相互作用的二项Logistic回归(二项式族参数和p值差异)

    当我在glm模型中使用交互加上family = binomial参数时以及当我将其遗漏时,我有一个关于重要性和重要性差异的问题 . 我对逻辑回归很陌生,过去只进行了更简单的线性回归 . 我有一个 observations of tree growth rings 的数据集,其中包含 two categorical explanatory variables ( Treatment 和 Origin...
  • 28 votes
     answers
     views

    从glm中提取pvalue

    我正在运行很多回归,我只对一个特定变量的系数和p值的影响感兴趣 . 因此,在我的脚本中,我希望能够从glm摘要中提取p值(获得系数本身很容易) . 我知道查看p值的唯一方法是使用summary(myReg) . 还有其他方法吗? 例如 . : fit <- glm(y ~ x1 + x2, myData) x1Coeff <- fit$coefficients[2] # only re...
  • 2 votes
     answers
     views

    在表格中呈现Tukey HSD成对p值

    我在我的数据上运行了一个posthoc Tukey HSD,它有十个因子级别 . 该表非常庞大,我希望在成对表中将p值呈现给读者,留下附录的45行表 . 这是一个示例数据集: set.seed(42) x <- rnorm(100,1,2) category <- letters[1:10] data <- cbind.data.frame(x, category) summar...
  • 0 votes
     answers
     views

    微阵列数据的多个配对t检验

    我正在研究来自同一患者的转移和主要肿瘤中数千个基因的mRNA表达数据 . 从多发转移的患者中,所有转移灶都获得单独的数据 . 该表看起来像下面的那个 . Subjects Type GeneA GeneB GeneC A Met 1 2 5 A Main 3 6 43 B Met 5 ...
  • 0 votes
     answers
     views

    如何获得ks测试和t检验的p值矩阵?

    我是R的新人 . 如果你能帮助我那将是伟大的 . 我的问题如下: 假设我有5个组,Group1,Group2,Group3,Group4和Group5,每个组包含100个数据点 . 现在我想使用t-test或ks-test将这些组相互比较,并希望生成p值矩阵 . 基本上,存在5x5的p值矩阵 . 我使用 corr.mat 函数完成了类似的相关工作 . 这里,5组只是为了说明的目的,在一天结束时,我...
  • 1 votes
     answers
     views

    相关矩阵与R中的显着性检验

    我想知道在R中是否已经有一个可以用SPSS或Statistica进行显着性测试的相关矩阵?目前,我不得不生成一个相关矩阵(使用 ggplot2 geom_tile )并使用 cor.test 分离显着性测试 . 将两者合并为一个图形会很棒 . 这是我正在使用的真实数据的子集 . X Var2 CorrValue pValue 1 AAAA ...
  • -1 votes
     answers
     views

    相关系数的Bootstrap p值(重采样方法)

    我有这个大的数据集(N = 300.000)并且通过功率分析我得出的结论是,如果它存在,我只需要250个观察来找到相关性 . 因此,我想使用自举来挑选1000个大小为n = 250的样本,以找到这1000个样本中的p值范围 . 我对bootstrap方法很不熟悉,但是在这里我给出了一个关于我在启动包中有多远的例子 . 我用Iris数据集来说明 . 我想要的输出是直方图,显示1000个获得的p值的频...
  • 0 votes
     answers
     views

    通过p值和皮尔逊相关项来评估预测

    我正在使用R对两种不同的预测模型进行一些评估 . 评估的基本思想是使用cor . ()函数比较Pearson相关性和相应的p值 . 下图显示了相关系数的最终结果及其p值 . 我们建议具有较低相关系数和相应较低p值(小于0.05)的模型更好(或者,相关系数更高但具有相当高的相应p值) . 所以,在这种情况下,总的来说,我们会说model1比model2更好 . 但问题是,是否还有其他具体的统计方法...
  • 0 votes
     answers
     views

    从相关矩阵中提取特定的成对相关性,并为相应的显着性水平添加适当数量的星号

    我正在从一个更大的矩阵创建一个具有特定相关性的(apa)表,如果相关性很大,我还想添加一个星号 . 我使用Hmisc包来创建所有可能的相关性和相应的p值 . 然后我使用MOTE包来舍入相关性并摆脱前导零 . 然后我将p值更改为星号 . 我把那些感兴趣的相关性拉出来并把它们放在一个新的矩阵中 . 假设我只想创建一个新的相关矩阵(3乘4),其中'am','gear'和'carb'定义3行,'mpg',...
  • 0 votes
     answers
     views

    Manova Strange P值

    我在这方面工作了一天,并希望你们能对这个奇怪的事情做出暗示 . y0,y1和y2通过相同的方法独立生成 . 它们通过相同的方法分成20组 . 然而,manova说它们有很大的不同?为什么? Manova测试(存储在变量s中)的摘要说: Pr(> F)值小于2.2e-16 . y0 <- runif(100, 0, 1) y1 <- runif(100, 0, 1) y2 <...
  • 0 votes
     answers
     views

    R:使用模拟计算p值

    我编写了这段代码来对两个随机分布的观察x和y运行一个测试统计 mean.test <- function(x, y, B=10000, alternative=c("two.sided","less","greater")) { p.value <- 0 alternative <- match.arg(alternat...

热门问题