similar to: package GLLMGibbs

I am trying to write a function that will simulate observed counts X and Y as below, I want the function to be able give many replicates, can somebody advise? thanks Oarabile #alpha,n sdx and sdy are constant, and N is a vector of length n. unstructured<-function(n,N,alpha,sdx,sdy){ Vx<-rnorm(n,mean=0,sd=sdx) Vy<-rnorm(n,mean=0,sd=sdy) thetax<-exp(alpha+Vx) thetay<-exp(alpha+Vy)

I'm trying to write a function that takes a vector of length n and then takes the first value of the vector i.e j=1 and forms a new vector of length n (i.e replicate the first value n times). This function will then calculate the absoulte difference of the original vector and the new vector and store the results omitting the difference between the value and itself. This function should

I am trying to determine whether two samples are identical or not. I'm aware that somebody can use the Kolmogorov-Smirnoff test to compare empirical distributions, but since my samples have ties I'm not sure if I'm getting the right p-values for the comparison. Can the Kolmogorov-Smirnoff test be adjusted for the case when ties exists and are there any functions that already

I'm trying to produce empirical bayes estimates based on the lognormal model in disease mapping Is there a way this can be done in R? thanks Oarabile

I want to do kernel smoothing on the standardised mortality rates at specified location(x-y coordinates) and be able to produce contour map. I've tried to use the package called fields in R (function smooth.2d) but the problem is I get negative results after smoothing and I do not understand why? Your contributions will be highly appreciated or if there is any package that can do this I

-------------- next part -------------- An embedded message was scrubbed... From: Martin Maechler <maechler at stat.math.ethz.ch> Subject: Re: producing maps Date: Fri, 25 Feb 2005 00:18:59 +0100 Size: 2511 Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050225/75bb1c86/producingmaps.mht

Hello, can anyone explain the following error I get when trying to use the CHECK command to check a new version of my pakcage under 1.4.0? ****** ./R CMD check ~/GLMMGibbs.0.5.1/GLMMGibbs * checking for working latex ... OK * using log directory `/homef/jonm/R-1.4.0/bin/GLMMGibbs.Rcheck' ... <Installs library, documentation, and then performs various tests, including the example,

Dear R-users, I wondered if anyone has extended GLMMGibbs to include correlated random effects, and if so, whether they would be willing to let me use their code? Jonathan Myles has no plans to extend glmm in this manner within the foreseeable future. With thanks, Patty -- -------------------------------------------------------------------------------- Assoc Prof Patty Solomon

Hi all I know GLMMGibbs is still in beta but has anyone experienced (and solved ;-) this problem? I decided to look at the seeds example but I get a core dump on two intel linux boxes and also a sun workstation. All are running R1.3.0 but different hardware/OS's so I think I've done something wrong > library(GLMMGibbs) > data(seeds) > seeds$plate <- as.factor(1:21) >

Dear All, I wonder if there is a way to fit a generalized linear mixed models (for repeated binomial data) via a direct Maximum Likelihood Approach. The "glmm" in the "repeated" package (Lindsey), the "glmmPQL" in the "MASS" package (Ripley) and "glmmGIBBS" (Myle and Calyton) are not using the full maximum likelihood as I understand. The

Dear all; Installing the GLMMGibbs package to my Solaris Unix box, I got an compiling error: ars.c:497:10: missing terminating " character ars.c: In function `dump_arse': ars.c:498: error: parse error before "mylesj" ..... The compiling error was reported to the list on Jul 3, 2003. According to Prof. Brain Ripley this is a known problem with the package and gcc 3.3,

I getting compilation errors when trying to install GLMMGibbs (see below). I'm running R v 1.7.1 on SuSE Linux 8.2. Has anyone else had this problem? I tried it on a Win2000/R 1.5.1 combination and it worked fine. Any hints are greatly appreciated. Thank you in advance, Damien >install.packages("GLMMGibbs") trying URL `http://cran.r-project.org/src/contrib/PACKAGES'

Greetings- I'm working with data that are multilevel in nature and have a dichotomous outcome variable (presence or absence of an attribute). As far as I can tell from reading archives of the R and S lists, as well as Pinheiro and Bates and Venables and Ripley, - nlme does not have the facility to do what amounts to a mixed-effects logistic regression. - The canonical alternative is

I''m attempting to use caches_page to speed up the loading of certain pages on my site, and have run into a couple questions: 1) The site has ''states'' (as in U.S. states), so I have a State model and a States controller. Going to http://sitename/states/ routes you to the States controller, ''index'' action. The States controller looks like so: class

Happy new month, everyone! I am planning to write a NIH grant proposal to study ways to speed Monte Carlo based maximum likelihood algorithm for hierarchical models with a focus on generalized linear mixed models (GLM with random effects). I thought it would be nice and also increase the chance of funding if I could produce an R package in the process. I understand that Prof. Pinheiro ang Bates

I am having trouble understanding how the 'by' function works. Using this bit of code: i <- data.frame(x=c(1,2,3), y=c(0,0,0), B=c("red","blue","blue")) j <- data.frame(x=c(1,2,3), y=c(1,1,1), B=c('red','blue','green')) plot(0, 0, type="n", xlim=c(0,4), ylim=c(0,1)) by(i, i$B, function(s){ points(s$x, s$y, col=s$B) })

I have a working dialplan for our phone system with Mon-Fri, business hours identification, etc. But what I'm lacking right now is support for company holiday dates. What I'd like to do is to create a database of these dates and just update them as new years rollover. I suspect others have done this sort of thing with Asterisk before, but I've not found any resources so far.

I''m following the Depot tutorial in the Rails book, and have just hit a wall. As instructed by the tutorial, I added the following code to the Product model: protected def validate errors.add(:price, "should be positive") unless price > 0.0 end When I attempt to add a new product, I''m getting an error: undefined method `>'' for nil:NilClass

The answer was (thanks to Mark Leeds) to do with the use of a factor instead of a vector. on [2017-08-05] at 08:57 Myles English writes: > I am having trouble understanding how the 'by' function works. Using > this bit of code: > > i <- data.frame(x=c(1,2,3), y=c(0,0,0), B=c("red","blue","blue")) > j <- data.frame(x=c(1,2,3),

Hello, I am trying to use glmm() in library GLMMGibbs, but I don't seem to have the package and it is not listed on CRAN. Is it no longer supported? (I am using R 2.1.1 on Windows.) Many thanks in advance for your help. Regards, Mark The information contained in this e-mail is confidential and...{{dropped}}