search for: exmp

Hi! Is there a way to get biplot and identify to work togheter. Having the output of prcomp I would like to draw a biplot that. Instead of plotting the sample (row-names) names plots some "pch" symbols. (thats easy with xlabs) But now I would like to add using identify the names to only some of the points. I have noticed that both biplot.prcomp and biplot.default does a lot of

...f the code Thanks Dimitris ####################################################### quant<-c(0.25, 0.5, 0.75) qrnox<-c("qrnox1", "qrnox2", "qrnox3") for (i in 1:3){ qrnox[i]<-rq(nox~factor(year)+factor(state)+pcinc+I(pcinc^2)+I(pcinc^3), tau=quant[i], data=exmp) } summary(qrnox1) -- View this message in context: http://www.nabble.com/Problem-to-male-an-Index-in-looping-tp19706138p19706138.html Sent from the R help mailing list archive at Nabble.com.

....num) for(i in 1:sample.size){ genocase[i] = rbinom(1, 1, prob[i]) } case.idx = sample(which(genocase==1),1000) control.idx = sample(which(genocase==0),1000) return(rbind(gt[case.idx,],gt[control.idx,])) } How I can modify code in a way that it will simulate: 50 % of causal SNPs** ( exmp. 24 causal variants and 24 non causal SNPs) 10 % of causal SNP (exmpl. 5 causal and 43 non causal SNPs) Thanks a lot for any suggestion. [[alternative HTML version deleted]]

...se[i] = rbinom(1, 1, prob[i]) > } > case.idx = sample(which(genocase==1),1000) > control.idx = sample(which(genocase==0),1000) > > return(rbind(gt[case.idx,],gt[control.idx,])) > } > > How I can modify code in a way that it will simulate: > 50 % of causal SNPs** ( exmp. 24 causal variants and 24 non causal SNPs) > 10 % of causal SNP (exmpl. 5 causal and 43 non causal SNPs) > > Thanks a lot for any suggestion. > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailin...