search for: jameslenihan

Dear Colleagues   I have the following code that generates a boxplot for one specific labtest:   boxplot.n(LBSTRESN~COHORT, main="Boxplot of laboratory data for XLXXX-XXX test=Creatinine", subset = LBTEST=="Creatinine", xlab = "Cohort Number", ylab = "Units = umol/L", varwidth=TRUE   I would like to know if there is a way to loop through the dataset and

I am running a BUGS function with following   schools.sim <-bugs(data,inits,                              parameters,                              model.file="schools.txt",                              n.chains=3,                              n.iter=1000,                              bugs.directory="E:/Rprograms")   My model.file IS in the directory

I have downloaded the Umacs (Universal Markov chain sampler) and submitted the following sample code from Kerman and Gelman.   s <-Sampler( J=8, sigma.y  =c(15,10,16,11,9,11,10,18),           y  =c(28, 8,-3,7,-1,1,18,12),      theta =Gibbs(theta.update,theta.init),           V =Gibbs(V.update,mu.init),         mu =Gibbs(mu.update,mu.init),         tau =Gibbs(tau.update,tau.init),       

I am running the sm.density.compare function amd I am getting the following error:   my code is  > sm.density.compare(LBSTRESN,COHORT,xlab="Units = umol/L"subset = LBTEST=="Creatinine")   Error in if (from == to) rep.int(from, length.out) else as.vector(c(from,  :   missing value where TRUE/FALSE needed   I do not understand the error and I have had no help when searching

I found this example in an Introductory R book in the chapter on Matrices and Arrays The array is > m [,1] [,2] [,3] [,4] [,5] [1,] 0 12 13 8 20 [2,] 12 0 15 28 88 [3,] 13 15 0 6 9 [4,] 8 28 6 0 33 [5,] 20 88 9 33 0 The code is #returns the minimum value of d[i,j], i !=j and the row attaining #the minimum, for square