BookmarkSubscribeRSS Feed
deleted_user
Not applicable
Hi SAS Procedures Group,

I am needing some help in analyzing an experiment in a vineyard in a generalized randon block design. I have 3 blocks, 6 light treatments and I have separate results from both sides of the grapevines (east and west). Each light treatment is replicated two times within each block, so I have a total of 36 plots. I have as response variable the sugar content of the fruit. I would like to know if someone could help me with the sas statements to analyze these experiment.
In advance thank you very much
4 REPLIES 4
data_null__
Jade | Level 19
This is the question that you should have the answer to before you apply any treatments, collect any data or spend any money. 🙂

I'm not a statistician but I know when I need one.

I'm not sure the following is correct, but it seems reasonable. Perhaps a statistician will offer some advice.

[pre]
/* Make sample data */
proc plan seed=1150759636;
factors
block = 3 ordered
light = 6 random
side = 2 ordered
plot = 1 of 36 perm
;
treatments sugar = 1 of 500;
output out=plan side cvals=('E' 'W');
run;
proc print;
run;
proc glm;
class block light side;
model sugar = block light block*light side(block);
test h=light e=block*light;
run;

proc mixed;
class block light side;
model sugar = light;
random block block*light side(block);
lsmeans light / pdiff;
run;
[/pre]
deleted_user
Not applicable
Dear Data_null,
Thank you very much. I am missing in your porgramathe replication within the block. I have really 72 data, because I have 6 treatments, 3 blocks, 2 replicate of the same light treatment, and both sides of the canopy, East and West. The other question I have, but maybe as you said I have to ask a statistician, is how worry I have to be here about normality and homegeinity of variance.
Maybe you have a program in order to consider the 72 data.

Thanks again
data_null__
Jade | Level 19
> 2 replicate of the same light treatment,

You say you have 36 plots. Does this imply 72 plots? Or 2 samples on the sample of the same plot.

You will need to "draw a picture" I don't understand
data_null__
Jade | Level 19
This seems reasonable but you really need to get a statistician to help you...

[pre]
proc plan seed=1150759636;
factors
block = 3 ordered
rep = 2 ordered
light = 6 random
plot = 1 of 36 perm
side = 2 ordered
i = 1 of 72 perm
/ noprint
;
treatments sugar = 1 of 500;
output out=plan side cvals=('E' 'W');
run;
proc print;
run;
proc glm;
class block rep light side;
model sugar = block light block*light block*light*rep side side*block;
test h=light block*light e=block*light*rep;
run;

proc mixed method=type3;
class block rep light side;
model sugar = light block*light side;
random block block*light*rep side*block;
lsmeans light / pdiff;
run;
[/pre]

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 4 replies
  • 743 views
  • 0 likes
  • 2 in conversation