Dear experts,
I try to run a F1-LD-F3 experimental design with one between-subject factor and three within-subject repeated factors, using the nonparametric ANOVA analyses on SAS. Yet despite the absence of errors or warnings, no outputs except for the Model Info, Class Level Info, Dimensions, and Number of Observations can be shown. I even tried to add ods html file or ods tests3 in the code, but the ods html codes received the warning of a lack of authorization, while the ods tests3 shows the error that tests3 cannot be created.
I thereby wonder how can I retrieve the parameter estimation statistics. Or where is the problem in my code causing such an issue? Your advice will be greatly appreciated. I attach a sample of my data here, along with my codes for your information. Thank you in advance for your valuable guidance~
The code:
ABMDATA.RANK WITH SUBJECT and TIME CLASS
Obs Runs Region CT RSUW RSUP RSSP Time Visits rank_Visits SUBJECT TimeCls12345678910
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 1 | 0 | 9438.0 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 2 | 0 | 9438.0 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 3 | 0.0833333333 | 19643.5 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 4 | 0.1666666667 | 23718.0 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 5 | 0.3333333333 | 37300.5 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 6 | 0.3333333333 | 37300.5 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 7 | 0.5 | 50259.5 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 8 | 0.5833333333 | 58853.0 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 9 | 0.75 | 71033.5 | 1 | 1 |
1 | SouthEu | HpHd | LrwLsw | LruLsu | LrsLss | 10 | 0.8333333333 | 77468.0 | 1 | 1 |
Looks like there's some kind of issue with the data and the model you're trying to fit. I'm not familiar enough with PROC MIXED to tell you what that exact problem is, but the SAS error message tells you where to start looking.
NOTE: An infinite likelihood is assumed in iteration 0 because of a nonpositive definite estimated R matrix for Subject 1.
Since this is statistical related, I'll also move your question to the stats forum.
@luhawkyesag wrote:
Dear experts,
I try to run a F1-LD-F3 experimental design with one between-subject factor and three within-subject repeated factors, using the nonparametric ANOVA analyses on SAS. Yet despite the absence of errors or warnings, no outputs except for the Model Info, Class Level Info, Dimensions, and Number of Observations can be shown. I even tried to add ods html file or ods tests3 in the code, but the ods html codes received the warning of a lack of authorization, while the ods tests3 shows the error that tests3 cannot be created.
I thereby wonder how can I retrieve the parameter estimation statistics. Or where is the problem in my code causing such an issue? Your advice will be greatly appreciated. I attach a sample of my data here, along with my codes for your information. Thank you in advance for your valuable guidance~
The code:
OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;libname ABMDATA v9 '/home/yezhang3/myfolders/';NOTE: Libref ABMDATA was successfully assigned as follows:Engine: V9Physical Name: /home/yezhang3/myfoldersoptions user=ABMDATA;proc import out=ORIGIN datafile='/home/yezhang3/myfolders/All_Interactions.xlsx' dbms=XLSX replace;run;NOTE: One or more variables were converted because the data type is not supported by the V9 engine. For more details, run withoptions MSGLEVEL=I.NOTE: The import data set has 896282 observations and 8 variables.NOTE: USER.ORIGIN data set was successfully created.NOTE: PROCEDURE IMPORT used (Total process time):real time 1:13.69user cpu time 1:13.56system cpu time 0.05 secondsmemory 3923.71kOS Memory 31144.00kTimestamp 01/21/2019 01:55:00 AMStep Count 18 Switch Count 4Page Faults 0Page Reclaims 1232Page Swaps 0Voluntary Context Switches 493Involuntary Context Switches 84Block Input Operations 0Block Output Operations 112664proc rank data=ABMDATA.ORIGIN out=ABMDATA.RANK;var Visits;ranks rank_Visits;run;NOTE: The data set ABMDATA.RANK has 896282 observations and 9 variables.NOTE: PROCEDURE RANK used (Total process time):real time 0.49 secondsuser cpu time 0.25 secondssystem cpu time 0.10 secondsmemory 25069.18kOS Memory 52932.00kTimestamp 01/21/2019 01:55:01 AMStep Count 19 Switch Count 4Page Faults 0Page Reclaims 5895Page Swaps 0Voluntary Context Switches 2585Involuntary Context Switches 3Block Input Operations 112416Block Output Operations 126736proc print data=ABMDATA.RANK(obs=10);title "Subset of ABMDATA.RANK";run;NOTE: There were 10 observations read from the data set ABMDATA.RANK.NOTE: PROCEDURE PRINT used (Total process time):real time 0.04 secondsuser cpu time 0.04 secondssystem cpu time 0.01 secondsmemory 3453.84kOS Memory 30888.00kTimestamp 01/21/2019 01:55:01 AMStep Count 20 Switch Count 1Page Faults 0Page Reclaims 896Page Swaps 0Voluntary Context Switches 15Involuntary Context Switches 1Block Input Operations 288Block Output Operations 8NOTE: There were 896282 observations read from the data set ABMDATA.RANK.NOTE: The data set ABMDATA.SUB has 896282 observations and 11 variables.NOTE: DATA statement used (Total process time):real time 0.69 secondsuser cpu time 0.42 secondssystem cpu time 0.10 secondsmemory 3670.59kOS Memory 33456.00kTimestamp 01/21/2019 01:55:01 AMStep Count 21 Switch Count 3Page Faults 0Page Reclaims 655Page Swaps 0Voluntary Context Switches 2420Involuntary Context Switches 55Block Input Operations 126208Block Output Operations 154632proc print data=ABMDATA.SUB(obs=10);title "ABMDATA.RANK WITH SUBJECT and TIME CLASS";run;NOTE: There were 10 observations read from the data set ABMDATA.SUB.NOTE: PROCEDURE PRINT used (Total process time):real time 0.04 secondsuser cpu time 0.04 secondssystem cpu time 0.00 secondsmemory 2073.31kOS Memory 31400.00kTimestamp 01/21/2019 01:55:01 AMStep Count 22 Switch Count 1Page Faults 0Page Reclaims 302Page Swaps 0Voluntary Context Switches 15Involuntary Context Switches 1Block Input Operations 288Block Output Operations 32proc sort data=ABMDATA.SUB out=ABMDATA.SUB;by SUBJECT CT RSUW RSUP TimeCls;run;NOTE: There were 896282 observations read from the data set ABMDATA.SUB.NOTE: The data set ABMDATA.SUB has 896282 observations and 11 variables.NOTE: PROCEDURE SORT used (Total process time):real time 0.48 secondsuser cpu time 0.53 secondssystem cpu time 0.18 secondsmemory 139025.35kOS Memory 169060.00kTimestamp 01/21/2019 01:55:02 AMStep Count 23 Switch Count 7Page Faults 0Page Reclaims 34214Page Swaps 0Voluntary Context Switches 6541Involuntary Context Switches 7Block Input Operations 154368Block Output Operations 154632proc mixed data=ABMDATA.SUB anovaf method=mivque0 noprofile plots(maxpoints=none);CLASS CT RSUW RSUP TimeCls;MODEL rank_Visits = CT|RSUW|RSUP|TimeCls / CHISQ;Repeated RSUW*RSUP*TimeCls / SUB=SUBJECT TYPE=CS GRP=CT;LSMEANS CT*RSUW*RSUP*TimeCls / pdiff;run;NOTE: 282 observations are not included because of missing values.NOTE: An infinite likelihood is assumed in iteration 0 because of a nonpositive definite estimated R matrix for Subject 1.NOTE: PROCEDURE MIXED used (Total process time):real time 45.85 secondsuser cpu time 40.90 secondssystem cpu time 4.90 secondsmemory 2253651.15kOS Memory 2284604.00kTimestamp 01/21/2019 01:55:48 AMStep Count 24 Switch Count 30Page Faults 0Page Reclaims 57464Page Swaps 0Voluntary Context Switches 1083Involuntary Context Switches 51Block Input Operations 154656Block Output Operations 9141280OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
ABMDATA.RANK WITH SUBJECT and TIME CLASS
Obs Runs Region CT RSUW RSUP RSSP Time Visits rank_Visits SUBJECT TimeCls12345678910
1 SouthEu HpHd LrwLsw LruLsu LrsLss 1 0 9438.0 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 2 0 9438.0 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 3 0.0833333333 19643.5 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 4 0.1666666667 23718.0 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 5 0.3333333333 37300.5 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 6 0.3333333333 37300.5 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 7 0.5 50259.5 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 8 0.5833333333 58853.0 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 9 0.75 71033.5 1 1 1 SouthEu HpHd LrwLsw LruLsu LrsLss 10 0.8333333333 77468.0 1 1
Registration is now open for SAS Innovate 2025 , our biggest and most exciting global event of the year! Join us in Orlando, FL, May 6-9.
Sign up by Dec. 31 to get the 2024 rate of just $495.
Register now!
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.