similar to: poisson mean hypothesis

Dear useRs, I very much like the effect display of the proportional odds model on page 29 (Figure 8) of the following paper by John Fox: http://socserv.mcmaster.ca/jfox/Papers/logit-effect-displays.pdf It really gives a very concise overview of the model. I would like to use it to illustrate the proportional odds mixed models we fit here for a project on Diabetes but I can't seem to reproduce

Dear R-users, Does anyone of you have binary files for the packages rpart.permutation, snow and rsprng. I would like to use them in my classification tree. I know they are still at the 0.x development stage, though. Where can I get information on how to compile tar.gz files? I'm using windows XP with R 2.1.0. Thanks, Jan --------------------------------------------------------------------

Dear R-users, I tried to install the R(D)COM server version 1.35, to connect MS Excel and R. I am using R-version 2.1.0 on a Windows environment. And I am using MS Excel 2000 in a dutch version. I get the following error when I open Excel: "Compile error, can not find the project or the library" (translated from dutch). And it indicates LEFT in blue bold. Does anyone of you know what

Dear R-users, I would like to use the package flexmix to fit latent classes to a regression model. My data are repeated measurements of bernouilli variables so I can use the binomial family link to the glm function. The design is not balanced, meaning that for some individuals in my data set I have 10 measurements or more, for others I only have 5 or even less. My question is the following. Can

Craig, Thanks for your follow-up note on using the asypow package. My problem was not only constructing the "constraints" vector but, for my particular situation (Poisson regression, two groups, sample sizes of (1081,3180), I get very different results using asypow package compared to my other (home grown) approaches. library(asypow) pois.mean<-c(0.0065,0.0003) info.pois <-

Hi I'm quite new to optimization algorithms and I could use some advice or pointers. I'm using ?optim (method L-BFGS-B) to optimize a function over a 60-dimensional parameter space. The function itself takes about 1 to 6 minutes to compute. It finds an optimum after 6 tot 24 hours, depending on the problem I need to solve. For me, speed is an issue and I would like to reduce this to 1

Hi, We recently switched from R2.7.0 to R2.8.1 but having problems tracking down this 'invalid multibyte character' encoding issue. Can someone point us how to solve this? > sessionInfo() R version 2.8.1 (2008-12-22) x86_64-unknown-linux-gnu locale:

I am attempting to run a Poisson EWMA model using Patrick Brandt's source code. I get the following error when I run the code: Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent Dimnames(x) looks like this: [[1]] NULL [[2]] [1] "mip" "div" "nom" "unity" "mood"

Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary = TRUE,border=0.1) map.axes() filled.contour(mslp, zlim=c(1000,1020),color.palette = colorRampPalette(c("blue", "white",

Dear helpeRs, I would like to specify a newline command in R and pass it to latex via Sweave such that it corresponds to latex' \\ command. But that doesn't seem to be possible. If I Sweave the \n character, it just makes a new line in latex but not the \\ command. Is there a way such that the following code would result in latex in blablabla \\ blablabla on different line

Dear list, I am using glm's to predict count data for a fish species inside and outside a marine reserve for three different methods of monitoring. I run glms and figured out the best model using step function for each methods used. I predicted two values for my fish counts inside and outside the reserve using means of each of the covariates (using predict() ) therefore I have only one value

Hi R-helpers I'm looking for a vectorised function which does missing value replacement as in last observation carried forward in the zoo package but instead of a locf, I would like the locf function to add +1 to each time a missing value occurred. See below for an example. > require(zoo) > x <- 5:15 > x[4:7] <- NA > coredata(na.locf(zoo(x))) [1] 5 6 7 7 7 7 7 12 13

Hello, I just managed to compile openssh-3.6.1p1 on Ultrix/MIPS. One of the fixes needed for this is the following. It's because of: --- According to the Changelog - markus at cvs.openbsd.org 2002/12/10 19:26:50 [packet.c] move tos handling to packet_set_tos; ok provos/henning/deraadt --- This IP_TOS is now nicely stuffed in a function, but that function should be #defined

Dear Sir/Madam, I apologize for any cross-posting. I got a simple question, which I thought the R list may help me to find an answer. Suppose we have Y_1, Y_2, ., Y_n ~ Poisson (Lambda_i) and Lambda_i ~Gamma(alpha_i, beta_i). Empirical Bayes Estimator for hyper-parameters of the gamma distr, i.e. (alpha_t, beta_t) are needed. y=c(12,5,17,14) n=4 What about a Hierarchal B ayes

Hi everyone I am stuck on specifying my own maximum likelihood function for a special poisson model. My poisson model is as follow: O ~ Pois(b*N + b*RR*E) With O = observed cases b = constant (known) N = number of unexposed persons (known) E = number exposed persons (known) RR = relative risk (value is assumed under a scenario, e.g. RR=2.0) I used rpois to simulate the values of O for several

Dear R Users, I'm working on a problem where I have a multivariate response vector of counts and a continuous predictor. I've thought about doing this the same way you would do a Multvariate regression model with normally distributed data, but since these data are counts, they are probably better modeled with a Poisson distribution. For example y1<-rpois(100,3.5) y2<-rpois(100,1.5)

Hi, I'm totally new to R so I apologise for the basic request. I am looking at the incidence of a disease over two time periods 1990-1995 and 2003-2008. I have counts for each year, subdivided into three disease categories and by males/females. I understand that I need to analyse the data using poisson regression and have managed to use the pois.daly function to get age-sex adjusted rates and

Dear, I tried to install the DBI package in R-1.7.1, but this gave the following error: /volume1/scratch/jallemee/R/lib/R/bin/INSTALL: line 1: 1856 Done ( echo "options(save.image.defaults=${save_image_defaults})"; if test -s R_PROFILE.R; then cat R_PROFILE.R; fi; echo

Dear, I have problems with installing R-1.7.1 for Redhat 9. When I applied `./configure' and `make', I get the following error: make check make[1]: Entering directory `/volume1/scratch/jallemee/R-1.7.1/tests' make[2]: Entering directory `/volume1/scratch/jallemee/R-1.7.1/tests' make[3]: Entering directory `/volume1/scratch/jallemee/R-1.7.1/tests/Examples' make[4]: Entering

Hi, I'm sorry for this not being related to R but I think this is a good place to ask. I'm looking for DOE examples(experiments) that can be done at home or in class, such as Paper Helicopter, Paper Towel etc.. I'm thankful for any comment. Thomas [[alternative HTML version deleted]]