我是否应该告诉SAS在每个时间进行非线性回归使用NLIN的重复次数？

问题描述

我有一个关于植物物种的伪茎高（代码中称为 pseudostem ）与热时间（代码中名为 gdd ）的数据集。对于每个事件（对于每个gdd），我有几次重复的测量（在代码中称为 rep ）。为了执行折断回归，我应该告知SAS重复，还是重复（因为多个度量具有相同的gdd）。我使用下面的代码：

数据

data pseudostemsum;
  input 
gdd rep pseudostem; datalines;
33.46 1 .
33.46 2 .
33.46 3 .
33.46 4 .
33.46 5 .
33.46 6 .
33.46 7 .
33.46 8 .
33.46 9 .
33.46 10 .
33.46 11 .
33.46 12 .
33.46 13 .
33.46 14 .
33.46 15 .
33.46 16 .
33.46 17 .
33.46 18 .
33.46 19 .
33.46 20 .
93.23 1 .
93.23 2 . 
93.23 3 .
93.23 4 12.8
93.23 5 19.9
93.23 6 9.8
93.23 7 16.6
93.23 8 20.8
93.23 9 26.5
93.23 10 22.5
158.61 1 12.4
158.61 2 15.7
158.61 3 13.1
158.61 4 9.9
158.61 5 16.6
158.61 6 12.5
158.61 7 10.3
158.61 8 14.1
158.61 9 12.3
158.61 10 9.4
158.61 11 16.9
158.61 12 13.9
158.61 13 10
158.61 14 7.8
158.61 15 11.5
158.61 16 13.4
226.98 1 17.3
226.98 2 11.8
226.98 3 14.5
226.98 4 16
226.98 5 11
226.98 6 18.8
226.98 7 18.6
226.98 8 17
226.98 9 18.5
226.98 10 9.1
226.98 11 17.5
226.98 12 18.8
226.98 13 21
226.98 14 14
226.98 15 17.7
226.98 16 16.5
325.84 1 19.7
325.84 2 16.1
325.84 3 29.4
325.84 4 22.6
325.84 5 17.6
325.84 6 29
325.84 7 27.8
325.84 8 20.6
325.84 9 24
325.84 10 18.2
325.84 11 18
325.84 12 25.8
325.84 13 14.4
325.84 14 20.9
325.84 15 28.5
410.01 1 22.5
410.01 2 20.9
410.01 3 25.3
410.01 4 25
410.01 5 13.9
410.01 6 27
410.01 7 19.1
410.01 8 28.9
410.01 9 27.2
410.01 10 24
410.01 11 26.1
410.01 12 30.4
410.01 13 24
410.01 14 22
410.01 15 30.2
509.95 1 26.3
509.95 2 19.8
509.95 3 18.9
509.95 4 26.3
509.95 5 30
509.95 6 22.5
509.95 7 23
509.95 8 29.8
509.95 9 27.5
509.95 10 29.5
509.95 11 20.7
509.95 12 26
509.95 13 27.7
509.95 14 25.1
509.95 15 23.5
622.94 1 24
622.94 2 20.7
622.94 3 21.5
622.94 4 28.5
622.94 5 29.9
622.94 6 22
622.94 7 24.8
622.94 8 30.2
622.94 9 27
622.94 10 28.5
622.94 11 29
622.94 12 24
622.94 13 26.5
622.94 14 29.5
622.94 15 28.5
699.53 1 27.6
699.53 2 25
699.53 3 25.9
699.53 4 29
699.53 5 26.8
699.53 6 17.7
699.53 7 18.7
699.53 8 27.7
699.53 9 26.5
699.53 10 28.8
699.53 11 14.1
699.53 12 25.3
699.53 13 23.7
699.53 14 23.5
699.53 15 32.4
774.97 1 25.1
774.97 2 29
774.97 3 26.1
774.97 4 25
774.97 5 22.4
774.97 6 26.3
774.97 7 19.4
774.97 8 28.2
774.97 9 26.7
774.97 10 34.9
774.97 11 28.1
774.97 12 30
774.97 13 26.7
774.97 14 22.5
774.97 15 18.5
872.68 1 19
872.68 2 18.8
872.68 3 12.2
872.68 4 18.7
872.68 5 23.9
872.68 6 19.8
872.68 7 16.2
872.68 8 23.4
872.68 9 25.4
872.68 10 22.8
872.68 11 22
872.68 12 27.8
872.68 13 17.8
872.68 14 24.2
872.68 15 27.5
998.51 1 22.8
998.51 2 24
998.51 3 22.7
998.51 4 15.2
998.51 5 26.8
998.51 6 20
998.51 7 20
998.51 8 22.5
998.51 9 22.5
998.51 10 20
998.51 11 25.5
998.51 12 22
998.51 13 27.5
998.51 14 27
998.51 15 16
;

NLIN

proc nlin data=PSEUDOSTEMSUM plots (stats=none);

  parms a1 = 0.029 b1 = 13 gdd0 = 415;

  if (gdd<gdd0) then
    mean = a1 * gdd + b1;
  else
    mean = a1 * gdd0 + b1;

  model pseudostem = mean;
run;

即使我在代码中写了重复（rep 1,2,3 ....），我也没有说明任何重复。输出似乎还可以，但是我对此不太确定（我是使用SAS的新手。）

感谢您的宝贵时间！

亲切的问候，

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）