[R] Converting SAS Code

[Kevin E. Thorpe]

Regarding point 3, as a moderator I have been helping Andrew get this 
post out to the list over the past week. His previous attempts were 
encoded in some way that the listserv rejected. He sent me the post via 
his gmail account and viewing the source I saw it had at least both 
plain test and HTML an I said it was worth a try to post it. Certainly 
on my mail client his post displays acceptably with the notice that the 
HTML alternative was removed.

Kevin

On 09/29/2017 09:51 AM, Michael Dewey wrote:
> You might get better answers if you
>
> 1 - break this down into separate issues
> 2 - tell us what you want to achieve in words rather than SAS, we all 
> read English but few of us speak SAS
> 3 - post in plain text not HTML as HTML mangles your post
>
> On 29/09/2017 13:47, Andrew Harmon wrote:
>> Hello all,
>>
>> My statistical analysis training up until this point has been 
>> entirely done
>> in SAS. The code I frequently used was:
>>
>> *Yield Champagin;
>>
>> data yield;
>>
>> set stress;
>>
>> if field='YV' then delete;
>>
>> if field='HB' then delete;
>>
>> if barcode='16187DD4015' then delete;
>>
>> if barcode='16187DD6002' then delete;
>>
>> if barcode='16187DD2007' then delete;
>>
>> if barcode='16187DD5016' then delete;
>>
>> if barcode='16187DD8007' then delete;
>>
>> if barcode='16187DD7010' then delete;
>>
>> if barcode='16187DD7007' then delete;
>>
>> if barcode='16187DD8005' then delete;
>>
>> if barcode='16187DD6004' then delete;
>>
>> if barcode='16187DD5008' then delete;
>>
>> if barcode='16187DD7012' then delete;
>>
>> if barcode='16187DD6010' then delete;
>>
>> run; quit;
>>
>>
>>
>> Title'2016 Asilomar Stress Relief champagin yield';
>>
>> proc mixed method=reml data=yield;
>>
>> class rep Management Foliar_Fungicide Chemical_Treatment;
>>
>> model Grain_Yield__Mg_h_ =Management|Foliar_Fungicide|Chemical_Treatment
>> Final_Stand__Plants_A_ / outpred=resids residual ddfm=kr;
>>
>> random rep rep*Management rep*Management*Foliar_Fungicide;
>>
>> lsmeans Management|Foliar_Fungicide|Chemical_Treatment / pdiff;
>>
>> ods output diffs=ppp lsmeans=means;
>>
>> ods listing exclude diffs lsmeans;
>>
>> run; quit;
>>
>> %include'C:\Users\harmon12\Desktop\pdmix800.sas';
>>
>> %pdmix800(ppp,means,alpha=0.10,sort=yes);
>>
>> ods graphics off;
>>
>> run; quit;
>>
>> proc univariate data=resids normal plot; id Barcode Grain_Yield__Mg_h_
>> pearsonresid; var resid;
>> proc print data=resids (obs=3);run;
>>
>> Can someone please help me convert my code to R? Any help would be much
>> appreciated.
>>
>>
>> Thanks,
>>
>>
>> Andrew Harmon
>>
>>     [[alternative HTML version deleted]]
>>
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide 
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>> ---
>> This email has been checked for viruses by AVG.
>> http://www.avg.com
>>
>>
>

-- 
Kevin E. Thorpe
Head of Biostatistics,  Applied Health Research Centre (AHRC)
Li Ka Shing Knowledge Institute of St. Michael's Hospital
Assistant Professor, Dalla Lana School of Public Health
University of Toronto
email: kevin.thorpe at utoronto.ca  Tel: 416.864.5776  Fax: 416.864.3016