[R] Converting SAS Code

[Michael Dewey]

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
> 
> 

-- 
Michael
http://www.dewey.myzen.co.uk/home.html