Dear List,

Im new to R - making a transition from SAS. I have a space delimited file 
with the following structure. Each line in the datafile is identified by 
the first letter. 

A = Inventory (Inventory)
X = Stratum (Stratum_no Total Ye=year established)
P = Plot (Plot_no age slope= species)
T = Tree (tree_no frequency)
L = Leader (leader diameter height)
F = Feature (start_height finish_height feature)

On each of these lines there are some 'line specific' variables (in 
brackets). The data is hierarchical in nature - A feature belongs to a 
leader, a leader belongs to a tree, a tree belongs to a plot, a plot 
belongs to a stratum, a stratum belongs to inventory. There are many 
features in a tree. Many trees in a plot etc. 

In SAS I would read in the data in a procedural way using first. and last. 
variables to work out where inventories/stratums/plots/trees  finished and 
started so I could create summary statistics for each of them. For 
example, how many plots in a stratum? How many trees in a plot? An example 
of the sas code I would (not checked for errors!!!). If anybody could give 
me some idea on what the right approach in R would be for a similar 
analysis it would be greatly appreciated.

regards Andrew


Data datafile;
infile 'test.txt';
input @1 tag $1. @@;
retain inventory stratum plot tree leader;
if tag = 'A' then input @3 inventory $.;
if tag = 'X' then input @3 stratum_no $. total $. yearest $. ;
if tag = 'P' then input @3 plot_no $. age $. slope $. species $;
if tag = 'T' then input @3 tree_no $. frequency  ;
if tag = 'L' then input @3 leader_no $ diameter  height  ;
if tag = 'F' then input @3 start $ finish $ feature $;
if tag = 'F' then output;
run;
proc sort data = datafile;
by inventory stratum_no  plot_no  tree_no  leader_no;

* calculate mean dbh in each plot
data dbh
set datafile;
by inventory stratum_no  plot_no  tree_no leader_no
if first.leader_no then output;

proc summary data = diameter;
by inventory stratum plot tree;
var diameter;
output out = mean mean=;
run;

A BENALLA_1
X 1 10 YE=1985
P 1 20.25 slope=14 SPP:P.RAD
T 1 25
L 0 28.5 21.3528
F 0 21.3528 SFNSW_DIC:P
F 21.3528 100 SFNSW_DIC:P
T 2 25
L 0 32 23.1
F 0 6.5 SFNSW_DIC:A
F 6.5 23.1 SFNSW_DIC:C
F 23.1 100 SFNSW_DIC:C
T 3 25
L 0 39.5 22.2407
F 0 4.7 SFNSW_DIC:A
F 4.7 6.7 SFNSW_DIC:C
P 2 20.25 slope=13 SPP:P.RAD
T 1 25
L 0 38 22.1474
F 0 1 SFNSW_DIC:G
F 1 2.3 SFNSW_DIC:A
T 1001 25
L 0 38 22.1474
F 0 1 SFNSW_DIC:G
F 1 2.3 SFNSW_DIC:A
T 2 25
L 0 32.5 21.7386
F 0 2 SFNSW_DIC:A
F 2 3.3 SFNSW_DIC:G
F 3.3 10.4 SFNSW_DIC:C
X 2 10 YE=1985
P 1 20.25 slope=14 SPP:P.RAD
T 1 25
L 0 28.5 21.3528
F 0 21.3528 SFNSW_DIC:P
F 21.3528 100 SFNSW_DIC:P
T 2 25
L 0 32 23.1
F 0 6.5 SFNSW_DIC:A
F 6.5 23.1 SFNSW_DIC:C
F 23.1 100 SFNSW_DIC:C
T 3 25
L 0 39.5 22.2407
F 0 4.7 SFNSW_DIC:A
F 4.7 6.7 SFNSW_DIC:C
P 2 20.25 slope=13 SPP:P.RAD
T 1 25
L 0 38 22.1474
F 0 1 SFNSW_DIC:G
F 1 2.3 SFNSW_DIC:A
T 1001 25
L 0 38 22.1474
F 0 1 SFNSW_DIC:G
F 1 2.3 SFNSW_DIC:A
T 2 25
L 0 32.5 21.7386
F 0 2 SFNSW_DIC:A
F 2 3.3 SFNSW_DIC:G
F 3.3 10.4 SFNSW_DIC:C




	[[alternative HTML version deleted]]


