我正在使用Python尝试负二项式回归。我使用R和一个数据集找到了这个例子:
import pandas as pdimportstatsmodels.api as sm
df = pd.read_stata("http://www.karlin.mff.cuni.cz/~pesta/prednasky/NMFM404
我正在尝试使用Statsmodels模块将以下代码从R转换为Python:我使用pandas得到了一个类似的dataframe (名为df),目前在Python中有以下代码行(版本3.4,如果有区别的话):
model = sm.Poisson.from_formula("goals ~ att + def + hom
我对python是新手,我尝试做一个伽马回归,我希望得到类似的估计R,但是我不能理解python的语法,它会产生一个错误,一些如何解决它的想法。我的R代码:y = rgamma(18,10,.1)[1] 76.67251 140.40808 138.26660 108.20993 53.46417Residual deviance: 1.1761 on 17 degrees of freedom
Number of Fisher Scoring
我正在使用twang包创建倾向性分数,这些分数在使用survey::svyglm的二项式glm中用作权重。- get.weights(pscore, stop.method="es.mean")
这将生成以下警告: In eval(expr, envir, enclos) : non-integer #successes in a binomial <em