similar to: method of rpart when response variable is binary?

Dear R users, I’m performing some GLMMs analysis with a negative binomial link. I already performed such analysis some months ago with the lmer() function but when I tried it today I encountered this problem: *Erreur dans famType(glmFit$family) : unknown GLM family: ‘Negative Binomial’* Does anyone know if the negative binomial family has been removed from this function? I really appreciate any

i read a book from WOOD, there's an example which is talking about the pollutant. library(gamair) library(mgcv) y<-gam(death~s(time,bs="cr",k=200)+s(pm10median,bs="cr")+s(so2median,bs="cr")+s(o3median,bs="cr")+s(tmpd,bs="cr"),data=chicago,family=Possion) lag.sum<-function(a,10,11) {n<-length(a) b<-rep(0,n-11) for(i in 0:(11-10))

Make a plot in R where you compare the probability distributions for the Binomial distributions with N=1000 trials and succes probability 0.005 and the Possion Distribution with rate lambda=5. My answer is b<-binom(x, 1000, 0.005, log = FALSE) plot(b) p<-dpois(x, lambda, log = FALSE) plot(p) *Is this the right way to solve the problem?* -- View this message in context:

I think I should have posted this question here as well. I am posting my question here since it is R related. Please see below. I originally posted this to sci.stat.math "Nasser Abbasi" <nma at 12000.org> wrote in message news:Mm4gj.28843$R92.4987 at newsfe16.phx... > > I think R documentation is a bit hard for me to sort out at this time. > > I was wondering if

Hi I'm doing a glm analysis and I have two doubts (at least :) 1) When I run the function it gives a lot of warnings (see below) what they mean ? (may be I'm ignorant about this analysis ...) glm.poisson<-glm(log(Jkij+1)~fac.ano+fac.tri+fac.icesr+fac.mat+fac.ano:fac.icesr+fac.ano:fac.tri,family=poisson()) warnings() 40: non-integer x = 1.252763 41: non-integer x = 1.864785 42:

Hello, from a search of the archives and functions, I am looking for information on creating random correlated counts from a multivariate Poisson distribution.  I can not seem to find a function that does this. Perhaps, it has not yet  been created. Has anyone created an R package that does this.   thanks,   Jourdan Gold     [[alternative HTML version deleted]]

Hello, Anyone who knows a fast and accurate algorithm for generating draws from an underdispersed Poisson distribution. Or even better, if there is a package containing such an implementation. Thanks Rene

I am trying to find a definition for the two-tailed exact binomial test but have been unsuccessful. Can you help? --------------------------------- [[alternative HTML version deleted]]

Hi this is my first post. I am trying to run a simulation for a computer playing Von Neumann poker and adjusting it's expectation of an opponent's behavior according to how the opponent plays. This program involves random generation of "hands" and shifting of parameters. However, when I run the code, no errors come up, but the program doesn't do anything. Could someone

Dear all, I stumbled over a problem recently when trying to use srt with text() on a windows device. What I intended to do was to plot a simple regression line, and to rotate a piece of text such that the text has the same angle as the regression line. However, the text is always plotted in a slightly wrong angle: #### x=1:10 #create arbitrary x and y values y=x*2-rnorm(1:10)

Hello All, Am I misreading the documentation? The text.rpart documentation says: "label a column name of x$frame; values of this will label the nodes. For the "class" method, label="yval" results in the factor levels being used, "yprob" results in the probability of the winning factor level being used, and ?specific yval level? results in the probability of

I tried using rpart, as below, and got this error message "rpart Error in yval[, 1] : incorrect number of dimensions". Thinking it might somehow be related to the large number of missing values, I tried using complete data, but with the same result. Does anyone know what may be going on, and how to fix it? I have traced two similar error messages in the Archive, but following the

Hello, I am in the process of writing a Perl program to carry out and analyze large numbers of regressions using R as the regression engine, and I am using Statistics::R to create the communication bridge between the two programs. Statistics::R creates a pipe between R and Perl and uses Rterm.exe (hidden during run-time) to carry out the R commands. This communication bridge works well for

Hello: I am using rpart and would like more background on how the splits are made and how to interpret results - also how to properly use text(.rpart). I have looked through Venables and Ripley and through the rpart help and still have some questions. If there is a source (say, Breiman et al) on decision trees that would clear this all up, please let me know. The questions below pertain to a

Hi R-users, I have a collection of dataframes and know how to build a string that refers to it, in this example, name_infra_alg_inc. Then, I have a character string yval, which the user can select from a drop down list. It contains the column names of the dataframes. assign(paste(name_infra_alg_inc, "ci", sep="."), ci(get(name_infra_alg_inc)[[yval]], confidence=0.95)) My

Hi, I am using the rpart package to build a classification tree. I did manage to build a tree with data on a previous project. However, when attampting to build a tree on a project I am working on, I seem to be getting the error shown below: > nhg3.rp <- rpart(profitresp ~., nhg3, method="class") > plot(nhg3.rp, branch=0.4, uniform=T); text(nhg3.rp, digits=3) Error in

Thanks for the tip using a smoothing technique before plotiing in order to get a curve instead of a line connecting the observations. But that's not the solution for my main problem with the "unclean" line plot. In order to show my problem let's take this simple example: > xval <- c(1, 2, 3, 4, 5, 6, 7, 8) > yval <- c(10, 30, 40, 50, 70, 90, 100, 110) > plot

Hi. How do I plot the straight line and r? in a scatter plot using a simple file x~y? Sueli Rodrigues Eng. Agr?noma - UNESP Mestranda - USP/ESALQ PPG-Solos e Nutri??o de Plantas Fones (19)93442981 (19)33719762

I am trying to make a decision tree using rpart. The function runs very quickly considering the size of the data (1742, 163). When I call the summary command I get this: > summary(bookings.cart) Call: rpart(formula = totalRev ~ ., data = bookings, method = "class") n=1741 (1 observation deleted due to missingness) CP nsplit rel error 1 0 0 1 Error in yval[, 1] :

Hi all, I am using R1.5.0 under Unix, I have a couple of questions here. 1. My program is running out of memory. I am writing a program to grow a list of trees using rpart() on a subset of a large dataset(5807x693) with a different response for every tree. I saw that after each tree was constucted, 116 MB of data was being added to the Vcells. I have no idea what this data is. My dataset is