BookmarkSubscribeRSS Feed
luhawkyesag
Fluorite | Level 6

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: V9
Physical Name: /home/yezhang3/myfolders
 
options 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 with
options 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.69
user cpu time 1:13.56
system cpu time 0.05 seconds
memory 3923.71k
OS Memory 31144.00k
Timestamp 01/21/2019 01:55:00 AM
Step Count 18 Switch Count 4
Page Faults 0
Page Reclaims 1232
Page Swaps 0
Voluntary Context Switches 493
Involuntary Context Switches 84
Block Input Operations 0
Block Output Operations 112664
 
 
proc 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 seconds
user cpu time 0.25 seconds
system cpu time 0.10 seconds
memory 25069.18k
OS Memory 52932.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 19 Switch Count 4
Page Faults 0
Page Reclaims 5895
Page Swaps 0
Voluntary Context Switches 2585
Involuntary Context Switches 3
Block Input Operations 112416
Block Output Operations 126736
 
 
proc 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 seconds
user cpu time 0.04 seconds
system cpu time 0.01 seconds
memory 3453.84k
OS Memory 30888.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 20 Switch Count 1
Page Faults 0
Page Reclaims 896
Page Swaps 0
Voluntary Context Switches 15
Involuntary Context Switches 1
Block Input Operations 288
Block Output Operations 8
 
 
NOTE: 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 seconds
user cpu time 0.42 seconds
system cpu time 0.10 seconds
memory 3670.59k
OS Memory 33456.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 21 Switch Count 3
Page Faults 0
Page Reclaims 655
Page Swaps 0
Voluntary Context Switches 2420
Involuntary Context Switches 55
Block Input Operations 126208
Block Output Operations 154632
 
 
proc 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 seconds
user cpu time 0.04 seconds
system cpu time 0.00 seconds
memory 2073.31k
OS Memory 31400.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 22 Switch Count 1
Page Faults 0
Page Reclaims 302
Page Swaps 0
Voluntary Context Switches 15
Involuntary Context Switches 1
Block Input Operations 288
Block Output Operations 32
 
 
proc 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 seconds
user cpu time 0.53 seconds
system cpu time 0.18 seconds
memory 139025.35k
OS Memory 169060.00k
Timestamp 01/21/2019 01:55:02 AM
Step Count 23 Switch Count 7
Page Faults 0
Page Reclaims 34214
Page Swaps 0
Voluntary Context Switches 6541
Involuntary Context Switches 7
Block Input Operations 154368
Block Output Operations 154632
 
 
proc 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 seconds
user cpu time 40.90 seconds
system cpu time 4.90 seconds
memory 2253651.15k
OS Memory 2284604.00k
Timestamp 01/21/2019 01:55:48 AM
Step Count 24 Switch Count 30
Page Faults 0
Page Reclaims 57464
Page Swaps 0
Voluntary Context Switches 1083
Involuntary Context Switches 51
Block Input Operations 154656
Block Output Operations 9141280
 
 
 
OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
 

 

 

 

ABMDATA.RANK WITH SUBJECT and TIME CLASS

 Obs Runs Region CT RSUW RSUP RSSP Time Visits rank_Visits SUBJECT TimeCls12345678910

1SouthEuHpHdLrwLswLruLsuLrsLss109438.011
1SouthEuHpHdLrwLswLruLsuLrsLss209438.011
1SouthEuHpHdLrwLswLruLsuLrsLss30.083333333319643.511
1SouthEuHpHdLrwLswLruLsuLrsLss40.166666666723718.011
1SouthEuHpHdLrwLswLruLsuLrsLss50.333333333337300.511
1SouthEuHpHdLrwLswLruLsuLrsLss60.333333333337300.511
1SouthEuHpHdLrwLswLruLsuLrsLss70.550259.511
1SouthEuHpHdLrwLswLruLsuLrsLss80.583333333358853.011
1SouthEuHpHdLrwLswLruLsuLrsLss90.7571033.511
1SouthEuHpHdLrwLswLruLsuLrsLss100.833333333377468.011
2 REPLIES 2
Reeza
Super User

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: V9
Physical Name: /home/yezhang3/myfolders
 
options 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 with
options 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.69
user cpu time 1:13.56
system cpu time 0.05 seconds
memory 3923.71k
OS Memory 31144.00k
Timestamp 01/21/2019 01:55:00 AM
Step Count 18 Switch Count 4
Page Faults 0
Page Reclaims 1232
Page Swaps 0
Voluntary Context Switches 493
Involuntary Context Switches 84
Block Input Operations 0
Block Output Operations 112664
 
 
proc 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 seconds
user cpu time 0.25 seconds
system cpu time 0.10 seconds
memory 25069.18k
OS Memory 52932.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 19 Switch Count 4
Page Faults 0
Page Reclaims 5895
Page Swaps 0
Voluntary Context Switches 2585
Involuntary Context Switches 3
Block Input Operations 112416
Block Output Operations 126736
 
 
proc 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 seconds
user cpu time 0.04 seconds
system cpu time 0.01 seconds
memory 3453.84k
OS Memory 30888.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 20 Switch Count 1
Page Faults 0
Page Reclaims 896
Page Swaps 0
Voluntary Context Switches 15
Involuntary Context Switches 1
Block Input Operations 288
Block Output Operations 8
 
 
NOTE: 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 seconds
user cpu time 0.42 seconds
system cpu time 0.10 seconds
memory 3670.59k
OS Memory 33456.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 21 Switch Count 3
Page Faults 0
Page Reclaims 655
Page Swaps 0
Voluntary Context Switches 2420
Involuntary Context Switches 55
Block Input Operations 126208
Block Output Operations 154632
 
 
proc 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 seconds
user cpu time 0.04 seconds
system cpu time 0.00 seconds
memory 2073.31k
OS Memory 31400.00k
Timestamp 01/21/2019 01:55:01 AM
Step Count 22 Switch Count 1
Page Faults 0
Page Reclaims 302
Page Swaps 0
Voluntary Context Switches 15
Involuntary Context Switches 1
Block Input Operations 288
Block Output Operations 32
 
 
proc 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 seconds
user cpu time 0.53 seconds
system cpu time 0.18 seconds
memory 139025.35k
OS Memory 169060.00k
Timestamp 01/21/2019 01:55:02 AM
Step Count 23 Switch Count 7
Page Faults 0
Page Reclaims 34214
Page Swaps 0
Voluntary Context Switches 6541
Involuntary Context Switches 7
Block Input Operations 154368
Block Output Operations 154632
 
 
proc 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 seconds
user cpu time 40.90 seconds
system cpu time 4.90 seconds
memory 2253651.15k
OS Memory 2284604.00k
Timestamp 01/21/2019 01:55:48 AM
Step Count 24 Switch Count 30
Page Faults 0
Page Reclaims 57464
Page Swaps 0
Voluntary Context Switches 1083
Involuntary Context Switches 51
Block Input Operations 154656
Block Output Operations 9141280
 
 
 
OPTIONS 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

 

luhawkyesag
Fluorite | Level 6
Thank you, Reeza, for your effortful assistance and valuable guidance. I’ll try to figure out a solution following that direction. Also, I’ll wait for any feedback from the stats forum 🙂

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

Register now!

What is ANOVA?

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.

Discussion stats
  • 2 replies
  • 471 views
  • 0 likes
  • 2 in conversation