search for: zaihra

hi, a) i have something like: ecdfgrp1<-ecdf(subset(mydata,TMT_GRP==1)$Y); ecdfgrp2<-ecdf(subset(mydata,TMT_GRP==2)$Y); how can i plot the difference between these 2 step functions? i could begin with ecdfrefl<-function(x){ecdfgrp2(x)-ecdfgrp1(x);} ... what next? b) if i have a vector with repeated numeric values how can i get the subset without repeated values .e.g (0,4,0,2,2)

Dear R community, I am using a dataset and would like to define new variables using a R variable: e.g. for (i in 1:10){ dataset$i<-something } But this is not the right way, I get only one variable in "dataset"... How can I change this? Thank you! Georg. ************************ Georg Ehret Johns Hopkins Baltimore, MD [[alternative HTML version deleted]]

Hi, can anyone help me fit betabinomial model to the following dataset where each iD is a cluster in itself , if i use package aod 's betabinom model it gives an estimate of zero to phi(the correlation coeficient ) and if i fix it to the anova type estimate obtained from icc( in package aod) then it says system is exactly singular. And when i try to fit my loglikelihood by

I have a problem, would really appreciate if anyone could help me. I am trying to use bootstrap to resample a dataframe n then use that sample of dataframe to calcuate parameter estimates using betabinomial model. I am kind of lost coz the bootstrap function requires us to give statistic (A function which when applied to data returns a vector containing the statistic(s) of