DATA Step, Macro, Functions and more

model for GxE biplot

Reply
N/A
Posts: 1

model for GxE biplot

I'm a new sas user. I have a dataset from 2 environments with about 50 genotypes, I want to do a GxE biplot analysis and have done the proc gml for the AMMI model and things are not working well with the proc iml which is where i'm currently stuck. see the code:

proc iml;

use outres2;

read all into resid;

ngen=nrow(resid);

nenv=ncol(resid);

use stats2;

read var{mse} into msem;

read var{dfe} into dfem;

use stats3;

read var{nrep} into nrep;

call SVD (U,L,V,resid);    [error message- statement invalid]

minimo=min(ngen,nenv);

L=L[1:minimo,];

ss=(L##2)*nrep;

suma=sum(ss);

percent=((1/suma)#ss)*100;

minimo=min(ngen,nenv);

percenta=0;

do I=1 to minimo

df=(ngen-1)+(nenv-1)-(2*I-1);

dfa=dfa//df;

porceacu=percent [I, ];

percenta=percenta+porceacu;

percenac=percenac//percenta;

run;

any tips?

Ask a Question
Discussion stats
  • 0 replies
  • 268 views
  • 0 likes
  • 1 in conversation