我希望能够查看由 coxme 包生成的 lmekin 对象的 p 值。
例如。
model= lmekin(formula = height ~ score + sex + age + (1 | IID), data = phenotype_df, varlist = kinship_matrix)
我试过：
summary(model)
coef(summary(model))
summary(model$coefficient$fixed)
fixef(model)/ sqrt(diag(vcov(model)) #(Calculates Z-scores but not p-values)
但这些都不起作用。那么如何查看这个线性混合模型的 p 值呢？
我正在将一个短语读入R脚本作为参数。如果该短语包含星号（*），则我不希望脚本运行。
但是，使用grepl时，我在识别星号时遇到了问题。例如：
> asterisk="*"
> phrase1="hello"
> phrase2="h*llo"
> grepl(asterisk,phrase1)
[1] TRUE
> grepl(asterisk,phrase2)
[1] TRUE
的结果grepl(asterisk,phrase1)应为FALSE。有谁知道我如何grepl识别短语中是否有星号？
r ×2