擅长:python、mysql、java
<p>我知道我不想在这里提出后续问题,如果有人能告诉我如何在不单击“回答问题”的情况下继续讨论我的原始问题(带代码,不带字符数),我很乐意将此文本移到该部分;我道歉。在</p>
<p>我的最后一个问题是如何“允许独立于体重调查与死亡率的关系”。我的下一个问题是出于好奇,如果要扩展这个问题,比如说,我是否想独立于体重和身高来研究死亡率?在</p>
<p>我写了这段代码,对于我的数据,这些残差加起来等于0,但我只想和专家们核实一下,我会这样做,以备将来参考:</p>
<pre><code>import numpy as np
import statsmodels.formula.api as smf
import sys
dat = np.loadtxt(sys.argv[1],dtype={"names":("SpeciesName","Mortality","Height","Weight"),"formats":("S40","f4","f4","f4")})
mymodel = smf.ols("Mortality~Height+Weight",data=dat).fit()
Residues = list(mymodel.resid_pearson)
SpeciesList = list(dat["SpeciesName"])
for species,residue in zip(SpeciesList,Residues):
print species + "\t" + str(residue)
</code></pre>
<p>再次道歉,如果我写错了这一节;我不觉得这是一个新的问题,作为一个评论,我不能添加代码;如果这是更合适的话,我很乐意把这个问题变成一个新的问题。在</p>