如何解决从SAS中的Proc Mixed复制到python的变异系数C.V.?
我不是统计学家,我需要在python中转换一些SAS功能。我需要一个公式或python代码,可以提供此语句创建的变异系数(C.V。):
proc mixed covtest CL;
by loc;
class line rep;
model &var=line ;
random rep;
estimate 'testmean' intercept 1;
lsmean line/diff=all;
ods output diffs=db30;
ods output lsmeans;
ods listing exclude all;
我目前正在使用此python代码,但CV已关闭。数据集仅是3列-品种的Entrycode,Rep列,Plot列和yield列:
#Calculate the mean and CV
mean = round(new_data['YLD'].mean(),1)
x = new_data[['ENTRYCODE','REP']]
cv = gs.calccv(new_data,x,''YLD ~ C(ENTRYCODE) + C(REP)','YLD',mean)
def calccv(dataset,formula,var,varmean):
#########################################
# Purpose: #
# Calculate the Coefficient of Variance #
# Inputs: #
# dataSet - test dataset #
# X = factors for the X axis #
# formula - statistical model used #
# var - variable to be checked #
#########################################
model = ols(formula,data=dataset).fit()
y_predict = model.predict(x)
rsme_calc = rmse(dataset[var],y_predict)
cv = round(100 * (rsme_calc / varmean),1)
std = dataset[var].std()
cv = round(100 * (std/varmean))
return cv
任何帮助将不胜感激。
谢谢!
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。