@mahdiamah wrote:
Dear ballardw and Astounding
Thank you for your attentions.
I opened and copied all variables from my data to the model. The data set is correct too.I even tried for other responses. But the error is still the same.
{NOTE: The parameter repeat is assigned the default starting value of 1.0, because it is not
assigned an initial value with the PARMS statement.
WARNING: Parameter group_a is in the input data set.
WARNING: Parameter Time_num is in the input data set.
WARNING: Parameter group_b is in the input data set.
WARNING: Parameter Time_num is in the input data set.
ERROR: Dependent variable is not in the input data set.
NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE NLMIXED used (Total process time):
real time 0.05 seconds
cpu time 0.04 seconds
}
Here is an example of what I mean by "copy from the log" and paste into code box opened with the {I}
908 proc nlmixed data=work.theoph;
909 parms beta1=-3.22 beta2=0.47 beta3=-2.45
910 s2b1 =0.03 cb12 =0 s2b2 =0.4 s2=0.5;
911 cl = exp(beta1 + b1);
912 ka = exp(beta2 + b2);
913 ke = exp(beta3);
914 pred = dose*ke*ka*(exp(-ke*time)-exp(-ka*time))/cl/(ka-ke);
915 model conc ~ normal(pred,s2);
916 random b1 b2 ~ normal([0,0],[s2b1,cb12,s2b2]) subject=subject;
917 run;
ERROR: QUANEW Optimization cannot be completed.
NOTE: Optimization routine cannot improve the function value.
NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE NLMIXED used (Total process time):
real time 2.38 seconds
cpu time 2.00 seconds
This Error is likely because the example data set used is limited in content. Put please notice that the entire procedure is present from the Proc to the Run statement. That way it is easy to see every single option used.
If I change the spelling of the variable on the Model statement:
918 proc nlmixed data=work.theoph;
919 parms beta1=-3.22 beta2=0.47 beta3=-2.45
920 s2b1 =0.03 cb12 =0 s2b2 =0.4 s2=0.5;
921 cl = exp(beta1 + b1);
922 ka = exp(beta2 + b2);
923 ke = exp(beta3);
924 pred = dose*ke*ka*(exp(-ke*time)-exp(-ka*time))/cl/(ka-ke);
925 model conco ~ normal(pred,s2);
926 random b1 b2 ~ normal([0,0],[s2b1,cb12,s2b2]) subject=subject;
927 run;
ERROR: Dependent variable is not in the input data set.
NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE NLMIXED used (Total process time):
real time 0.01 seconds
cpu time 0.00 seconds
Which shows the exact same error you claim to be getting. I know the set used in my example Nlmixed code does not have a variable named Conco so the error is not unexpected. This is why I strongly suspect either you have a typo for the variable on the model statement or possibly are using an incorrect data set.
If I run Proc Contents on the dataset used in the Proc Nlmixed data= such as
proc contents data=work.theoph;
run;
Part of the output looks like:
Alphabetic List of Variables and Attributes
# Variable Type Len
3 conc Num 8
4 dose Num 8
1 subject Num 8
2 time Num 8
5 wt Num 8
Where in my case it is very easy to see that the variable name is not conco.
Which is why @PaigeMiller suggesting showing us the result of Proc contents on the data set used in Proc Nlmixed along with the code with messages from the log.
