小编DKa*_*yan的帖子

Python 中 GLM 的方差分析测试

我正在尝试获取 GLM 中每个协变量的 F 统计量和 p 值。在 Python 中,我使用 stats mode.formula.api 来执行 GLM。

formula = 'PropNo_Pred ~ Geography + log10BMI + Cat_OpCavity + CatLes_neles + CatRural_urban + \
        CatPred_Control + CatNative_Intro + Midpoint_of_study'

mod1 = smf.glm(formula=formula, data=A2, family=sm.families.Binomial()).fit()
mod1.summary()
Run Code Online (Sandbox Code Playgroud)

之后,我尝试使用 statsmodels.stats 中的方差分析对该模型进行方差分析测试

table1 = anova_lm(mod3)
print table1
Run Code Online (Sandbox Code Playgroud)

但是我收到一条错误消息:“GLMResults”对象没有属性“ssr”

看起来这个 anova_lm 函数只适用于线性模型,python 中有一个模块可以对 GLM 进行 anova 测试吗?

python statistics glm statsmodels

7
推荐指数
2
解决办法
5054
查看次数

检查说明元信息...注意

我正在R中开发一个程序包,运行devtools::check()时会收到以下说明。

checking DESCRIPTION meta-information ... NOTE Malformed Description field: should contain one or more complete sentences.

我没有在说明中使用包装名称或“包装”一词。我也使用完整的句子进行描述,但是我不断得到此注释。所以我想知道在这种情况下完整的句子是什么意思。

r package

7
推荐指数
2
解决办法
1329
查看次数

R - 如何从距离矩阵中获取匹配元素的行和列下标

我有一个整数向量vec1,我正在使用dist函数生成一个远程矩阵。我想获取距离矩阵中某个值的元素的坐标(行和列)。本质上,我想获得相距 d 距离的一对元素。例如:

vec1 <- c(2,3,6,12,17)
distMatrix <- dist(vec1)

#   1  2  3  4
#2  1         
#3  4  3      
#4 10  9  6   
#5 15 14 11  5
Run Code Online (Sandbox Code Playgroud)

说,我对向量中相距 5 个单位的一对元素感兴趣。我想得到坐标 1 是行和坐标 2 是距离矩阵的列。在这个玩具示例中,我希望

coord1  
# [1] 5
coord2
# [1] 4
Run Code Online (Sandbox Code Playgroud)

我想知道是否有一种有效的方法来获取这些值而不涉及将dist对象转换为矩阵或遍历矩阵?

r distance matrix euclidean-distance

6
推荐指数
1
解决办法
3895
查看次数