Hello programmers,   I am running a PLS analysis for a diet data and i want to know how i can retain the factors scores in PLS. I am very familiar with retaining factor scores of the PCA analysis for use in my regression model but how to retain factors that are derived from PLS is a little problematic. Attached is my pls code. proc pls data = zthree method = pls nfac=4 varss details; model iron calc = T_mean l_mean vc_mean vd_mean vr_mean cn_mean e_mean mp_mean mr_mean /solution; output out = pattern xscore = scorex yscore = scorey; run; ... View more

Hello statisticians,   I currently have two datasets.   In the first dataset, i carried out a PCA on the dataset and i want to retain  3 principal components . This dataset has 750 people (Subject IDs) The second dataset has 200 people (subject IDs) and contain my relevant outcomes of interest.   I would normally want to merge the two datasets by subject IDs. But the problem would be that the final dataset will contain only 200 samples (since the second dataset has 200 people) and iam afraid  that it will change the initial principal components that i got which is looking pretty good.   My question therefore is: Is there a way i can retain the values of the 3 principal components in the first dataset, so that when i merge the two data, i can still use those values in my analysis?   I would appreciate any response.  ... View more

Hello programmers,   I am trying to change the numbers in my variables to another set of numbers. I am using the if and then statement but it seems to be longer as i would be changing the values of about 50 variables. I am not an advanced sas user. Is there are way i can do this task faster with if --the statement or any other easier method? Any help will be appreciated. My codes are something like this. data one; set fooddata; if xa= 0 then xa1= 0; if xa= 1 then xa1= 0.033; if xa=2 then xa1=0.067; if xa=3 then xa1= 0.143; if xa=4 then xa1=0.429; if xa=5 then xa1=0.857; if xa=6 then xa1= 1; if xa=7 then xa1=2.5; if xa=8 then xa1=4.5; if xa=9 then xa1=7; run; ... View more

Hello programmers,   Please i would be glad to know what could be wrong with my code.   I am trying to plot my data to be the line format as opposed to being in the ugly looking format as in the pic. I'm i missing an option that needs to be specified? Do i need to specify line? data one; set two;run; goptions reset=global ; axis1 label =(f= 'arial/bo' h=1.9 "Age of patients" justify=c f= 'arial/bo' h= 1.3 "(in years)") order = (60 to 120 by 10) value=(f= 'arial/bo'h=1.3 "60.0" "70.0" "80.0" "90.0" "100.0" "110.0" "120.0"); axis2 label =(a =90 f= 'arial/bo' h=1.9 "Systolic Blood Pressure mmHg") order =(80 to 200 by 10) value=(a= 90 f= 'arial/bo'h=1.3 "80.0" "90.0" "100.0" "110.0" "120.0" "130.0" "140.0" "150.0" "160.0" "170.0" "180.0" "190.0" "200.0"); symbol1 value=dot color=black interpol=none h=1.2; symbol2 value=triangle color=black interpol=none h=1.5; Legend1 across=1 down=2 noframe position= (bottom right inside) mode=protect label =(f= 'arial/bo' h=1.4 "Obesity") value=(f= 'arial/bo' h=1.4 "0" "1"); proc gplot data= one; plot systolic*aage=obese/haxis=axis1 vaxis=axis2 legend=legend1; run; ... View more