Re: [R] sas to r

From: Uwe Ligges <ligges_at_statistik.uni-dortmund.de>
Date: Sat 17 Jul 2004 - 20:43:43 EST

Adaikalavan Ramasamy wrote:

> On Fri, 2004-07-16 at 23:18, Greg Adkison wrote:
>

>>I would be incredibly grateful to anyone who'll help me translate some 
>>SAS code into R code.

>
>
> Searching for "SAS code OR script OR translate" on
> http://maths.newcastle.edu.au/~rking/R/ gives a few results, one of
> which looks promising is

> http://tolstoy.newcastle.edu.au/R/help/04/04/0009.html and

I anticipated the above citation from 01 Apr 2004 (hint: 01 April !!!!!) when I looked at the thread and realized URLs were cited ...

Uwe

> http://tolstoy.newcastle.edu.au/R/help/04/02/0660.html
>
>

>>Say for example that I have a dataset named "dat1" that includes five 
>>variables:  wshed, site, species, bda, and sla.  I can calculate with the 
>>following SAS code the mean, CV, se, and number of observations of 
>>"bda" and "sla" for each combination of "wshed," "species," and "site," 
>>restricting the species considered to only three of several species in 
>>dat1 (b, c, and p).  Moreover, I can output these calculations and 
>>grouping variables to a dataset named "dat2" that will reside in RAM 
>>and include the variables  wshed, site, species, mBdA, msla, cBda, 
>>sBdA, ssla, nBda, and nsla.

>
>
> data(iris)
> attach(iris)
> iris[c(1,2,51,52,101,102), ]
> Sepal.Length Sepal.Width Petal.Length Petal.Width Species
> 1 5.1 3.5 1.4 0.2 setosa
> 2 4.9 3.0 1.4 0.2 setosa
> ...
> 51 7.0 3.2 4.7 1.4 versicolor
> 52 6.4 3.2 4.5 1.5 versicolor
> ...
> 101 6.3 3.3 6.0 2.5 virginica
> 102 5.8 2.7 5.1 1.9 virginica
>
>
>>tapply(Sepal.Length, Species, function(x) c( mean(x), sd(x)/mean(x),

>
> length(x) ))
> $setosa
> [1] 5.00600000 0.07041344 50.00000000
>
> $versicolor
> [1] 5.93600000 0.08695606 50.00000000
>
> $virginica
> [1] 6.58800000 0.09652089 50.00000000
>
>
>
>>proc sort data=dat1;
>>  by wshed site species;
>>proc means data=dat1 noprint mean cv stderr n;
>>  by wshed site species;
>>  where species in ('b', 'c', 'p');
>>  var BdA sla;
>>  output out=dat2
>>    mean=mBdA msla
>>    cv=cBdA csla
>>    stderr=sBdA ssla
>>    n=nBdA nsla;
>>
>>Thanks,
>>Greg
>>
>>______________________________________________
>>R-help@stat.math.ethz.ch mailing list
>>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>>PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>>

>
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


R-help@stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Jul 17 21:05:44 2004

This archive was generated by hypermail 2.1.8 : Fri 18 Mar 2005 - 02:36:41 EST