K*******5
发帖数: 98 | 1 我有2个方程，有两个未知parameters.其他variables已知。 google之后我用以下sas
code：
已知得是yi，ni， 未知的是pi, rhoeqls.
equations one and two are set to zero.
sas 一直run,没有给解。有没有用过的给些指点。谢谢
proc model data=outt;
eq.one = sum((yi-ni*pi)/(pi*(1-pi)*(1+(ni-1)*rhoeqls)));
eq.two = sum(((ni-1)/(1+(ni-1)*rhoeqls)**2)*
((2*((yi*log(yi/(ni*pi))+(ni-yi)*log((ni-yi)/(ni-ni*pi))))*(yi/ni-
pi)-1+(ni- 1) *rhoeqls)));
solve pi rhoeqls/solveprint;
id ni yi;
where x1=0;
run; | K*******5
发帖数: 98 | 2 有没有知道的。不用SAS，用R也行。具体怎么解呢？
还有我之前用的sum does not make sense. I want to sum vertically from 1 to k.
please ignore sum().
xiexie a.
sas
ni-
【在 K*******5 的大作中提到】
: 我有2个方程，有两个未知parameters.其他variables已知。 google之后我用以下sas
: code：
: 已知得是yi，ni， 未知的是pi, rhoeqls.
: equations one and two are set to zero.
: sas 一直run,没有给解。有没有用过的给些指点。谢谢
: proc model data=outt;
: eq.one = sum((yi-ni*pi)/(pi*(1-pi)*(1+(ni-1)*rhoeqls)));
: eq.two = sum(((ni-1)/(1+(ni-1)*rhoeqls)**2)*
: ((2*((yi*log(yi/(ni*pi))+(ni-yi)*log((ni-yi)/(ni-ni*pi))))*(yi/ni-
: pi)-1+(ni- 1) *rhoeqls)));
|
|
