search for: yik

...for two variables, I want to convert it to work for any number of variables. Any ideas? yk <- numeric(nrow(xk)) for (j in 1:nrow(xk)) { w1 <- abs(x[,1] - xk$x1[j]) w2 <- abs(x[,2] - xk$x2[j]) Ik1 <- which(w1 <= h[1]) Ik2 <- which(w2 <= h[2]) Ik <- intersect(Ik1, Ik2) YIk <- yvals[Ik] yk[j] <- mean(YIk) } The result is the one-dimensional vector yk. Does anyone know how to extend this code to cover any number of variables? Thanks in advance. dxc13 -- View this message in context: http://www.nabble.com/extending-code-to-handle-more-variables-tp1560423...

Hi, I'm using samba 3.0.6. I see computer_name folders created in homedir. This never happened before I switch from Domain to ADS. Could someone please give me a few pointers? Regards, Norman Zhang [global] workgroup = ARKONDOMAIN realm = HQ.ARKONNETWORKS.COM server string = Samba Server %v security = ADS obey pam restrictions = Yes

Hi folks I am using the lmer function (in the lme4 library) to analyse some data where individuals are clustered into sets (using the SetID variable) with a single fixed effect (cc - 0 or 1). The lmer model and output is shown below. Whilst the fixed effects are consistent with stata (using xtmixed, see below), the std dev of the random effect for SetID is very very small (3.5803e-05)compared to

Hi, I have a genotype data for both case and controls and would like to calculate the HW p-value. However, since the number of one genotype is 0, I got wired result. Would someone help me to figure it out? Or confirm it's right? Thanks a lot. ============ > library( "genetics" ) NOTE: THIS PACKAGE IS NOW OBSOLETE. The R-Genetics project has developed an set of enhanced