search for: anamika

...et[,3]=1))} Error: subscript out of bounds In addition: Warning messages: 1: numerical expression has 11627 elements: only the first used in: mincls:maxcls 2: numerical expression has 11627 elements: only the first used in: mincls:maxcls Any help will be greatly appreciated. Thanks, Anamika --------------------------------- [[alternative HTML version deleted]]

...quation. I would then like to use the simulated data to run the Beta-Binomial model with covariates on it. Appreciate any help. set.seed(111) k<-20 n<-60 x<-NULL p<-rbeta(k,3,3)# so that the mean nausea rate is alpha/(alpha+beta) for(i in 1:k) x<-cbind(x,rbinom(300,n,p[i])) Thanks Anamika [[alternative HTML version deleted]]

...se rates pi1, . . . , piC for site i. The individual variation in category response rates can be modeled by a Dirichlet distribution. Just wondering if I am thinking through this correctly. If so, could someone share some thoughts on how this could be done in R? Really Appreciate any help. Thanks Anamika [[alternative HTML version deleted]]

...12 [5,] 8 14 12 > n [1] 10 36 31 > l_cl_exact [,1] [,2] [,3] [1,] 0.44390454 0.2184996 0.2314244 [2,] 0.04667766 NaN 0.2454760 [3,] 0.05452433 0.1855618 NaN [4,] 0.44390454 0.4862702 0.1855618 [5,] 0.10115053 NaN 0.2184996 Thanks for your help. Anamika [[alternative HTML version deleted]]

...e 2 functions: library(car) library(mvtnorm) mu = c(0,0) sigma = matrix(c(20,0,0,45),nrow=2) z = rmvnorm(10000,mu,sqrt(sigma)) dataEllipse(z,levels=.95) car::ellipse(mu, sigma*qchisq(.05,2), col="blue", radius=sqrt(2 * qf(.975, 2, 9998)) ) Any help is appreciated. Thanks Anamika [[alternative HTML version deleted]]

...terval using Clopper > Pearson method. I am using the following code but it seems to get into a > loop and not get out of it, it goes on forever although I am looping it > only 10 times across 63 sites with 10 observations per site. I was hoping > to get some help. > > Thanks > Anamika > set.seed(111) #so that you can regenerate the same values > # another time if you need to > k<-63 > n<-10 > x<-NULL > eta<-rnorm(63,-0.085,.990) # generating data using Binomial Logit Normal > p<-exp(eta)/(1+exp(eta)) > for(i in 1:k) > x<-cbind(x,rb...

I am using a random intercept model with SITEID as random and NAUSEA as outcome. Thanks. Anamika I tried using a dataset without missing values and changed my model statement accordingly but still get the same error. Follwoing in an excerpt. > anal.data <- read.table("nausea.txt", header=T, sep="\t") > list(names(anal.data)) [[1]] [1] "SITEID" "NAU...

Hi: Is there a way we can set up WINBUGS to run 100 simulated datasets on the same model and output results? Or do we have to call in each dataset at a time and repeat the process 100 times manually? Thanks Anamika [[alternative HTML version deleted]]

Hi All: I was looking for some help to code the equation of an ellipse in R, given I have a Bivariate Normal distribution (X,Y) with mean (m1,m2) and var-covariance matrix (var1,cov12,cov12,var2). Thanks Anamika [[alternative HTML version deleted]]

Hi, I cant understand where I am going wrong.Below is my code.I would really appreciate your help. Thanks. > genfile<-read.table("c:/tina/phd/bs871/hw/genfile.txt",skip=1) > > #read in SNP data > snp.dat <- as.matrix(genfile) > snp.name <- scan("c:/tina/phd/bs871/hw/genfile.txt",nline=1,what="character") Read 100 items

...{alpha[j]~dnorm(0.0,1.0E-10)} > > tau ~ dgamma(0.001,0.001) > } > *Dataset:* > SITEID NAUSEA > 1 0 > 1 1 > 1 1 > 1 0 > 1 1 > 1 1 > 1 0 > 1 1 > 1 1 > 1 1 > 1 0 > 1 1 > 1 0 > 1 1 > 1 1 > 1 1 > 1 0 > 1 0 > 1 0 > 1 1 > > -Anamika > [[alternative HTML version deleted]]

...0.6 0.7 0.8 1.0 1 630 alpha0 -2.0 0.2 -2.3 -2.1 -2.0 -1.9 -1.7 1 330 p[1] 0.2 0.0 0.1 0.1 0.2 0.2 0.3 1 1000 p[2] 0.1 0.0 0.0 0.1 0.1 0.1 0.2 1 1000 p[3] 0.0 0.0 0.0 0.0 0.0 0.1 0.1 1 1000 etc > Thanks > Anamika [[alternative HTML version deleted]]

...30.13 1 3080101562 01 25.8 1 3080101564 01 31.07 1 3080101565 01 26.53 0 3080101566 01 29.3 0 30801 01694 01 28.4 1 3080101695 01 26.03 1 3080101696 01 27.43 1 3080101697 01 31.96 0 3080101698 01 28.76 0 3080101699 01 28.4 1 3080101700 01 29.33 0 3080101701 01 36.02 0 3080101768 01 31.05 1   Thanks Anamika   [[alternative HTML version deleted]]

...rnin.cycles, mcmc=mcmc.cycles, + thin=5, tune=0.5, beta.start=vec, user.prior.density=logpriorfun, logfun=TRUE, + mu=mu0, gshape=a0, grate=b0) + return(posterior) + }) Error in eval(expr, envir, enclos) : object "x1" not found > Any suggestions will be greatly appreciated. Thanks, Anamika --------------------------------- We won't tell. Get more on shows you hate to love [[alternative HTML version deleted]]

Trevor, Thanks for your reply. That doesnot help Any other suggestions? Anamika On Tue, Aug 10, 2010 at 11:58 AM, Trevor Davies <tdavies@mathstat.dal.ca>wrote: > I think you need to load is R2WinBUGS again. > > require(R2WinBUGS) > > Trevor > > > Hi All: > > > > I am getting an error while trying to plot in R2Winbugs > > > &...

Hey Guys, I just joined this list a couple days ago...what''s up with all the spam? I maintain about 8 or so googlegroups lists and have never received any <cross fingers>. Is this something to (unfortunately) look forward to? -- ______________________________________________________________ Clayton Dukes ______________________________________________________________

...("p","rho","sigma2"), n.chains=1, n.iter=12000, n.burnin=5000, bugs.directory="H://AChaudhuri/winbugs14/WinBUGS14", working.directory=NULL) } Error in FUN(X[[1L]], ...) : .C(..): 'type' must be "real" for this format Thanks for your help. Anamika [[alternative HTML version deleted]]

Hi, I am trying to runn a age-period-cohort model, but here is what I am having problem with, hope you can help me! This is what I am trying to do: sumzero_a<-((A-min(A))/5+1) - mean((A-min(A))/5+1) where A is my age variable (numeric, the mid-point of a five-year age group), but I got the following error: Error in min(..., na.rm = na.rm) : invalid 'mode' of argument I am pretty

Hi, Anybody know what kind of tools are available in windows NT/2K for tracing win32 api calls ? Something on the lines of truss/strace. Apimon.exe, apispy just don't cut it. Not as robust as truss/strace. What kind of tracing tools does the wine development team use I wonder ? Thanks, -A

*How early sex affects teenagers ? Girls who have sex at an early age are at slightly greater risk than their peers for feeling <http://tomorrowsmedicinetoday.blogspot.com/2007/06/early-sex-affects-teenagers.html> * *depressed, a new study has found<http://tomorrowsmedicinetoday.blogspot.com/2007/06/early-sex-affects-teenagers.html> * *FOR MORE INFORMATION LOG