similar to: About zero-inflation poisson model

Dear Mr. or Ms.,   I used the R-software to run the zero-inflatoin negative binomial model (zeroinfl()) .   Firstly, I introduced one dummy variable to the model as an independent variable, and I got the estimators of parameters. But the results are not satisfied to me. So I introduced three dummy variables to the model. but I could not get the results. And the error message is

Hi All, My dependent variable is a ratio that takes a value of 0 (zero) for 95% of the observations and positive non-integer values for the other 5%. What model would be appropriate? I'm thinking of fitting a GLM with a Poisson ~. Now, becuase it takes non-integer values, using the glm function with Poisson family issues warning messages. Warning messages: 1: In dpois(y, mu, log = TRUE) :

Dear Mr. or Ms.,   I used the R-software to run the zero-inflatoin negative binomial model, but I could not get the results. And the error message is "solve.default(as.matrix(fit$gaussian)) ". In the model, I introduced 3 dummy variables. I do not know the reasons fully.   I will be very appreciate if you can help me. Thank you.   Best regards,   Sincerely, Xiongqing Zhang

Dear Mr. or Ms.,   I would like to use the R-software to run the zero inflated mixed model and hurdle mixed model. But I do not know how to do? Would you please tell me the code and data format?   I will be very appreciated if you can help me. Thank you very much.   Best regards,   Sincerely, Xiongqing Zhang [[alternative HTML version deleted]]

Hi, I try to use the visual studio 2008(VC.net) to call R code, but there is en error:  Exception from HRESULT: 0x80040013 How can I fix the problem? Thanks.   My PC systern is XP, R version is "R-2.14.1" and "R_Scilab_DCOM3.0-1B5.exe".   best regards, Xiongqing [[alternative HTML version deleted]]

Hi all, First post for me here, but I have been reading on the forum for almost two years now. Thanks to everyone who contributed btw! I have a dataset of 4000 observations of count of a mammal and I am trying to predict abundance from a inflated-zero model as there is quite a bit of zeros in the response variable. I have tried multiple options, but I might do something wrong as every

Dear all, I am analyzing counts of seabirds made from line transects at sea. I have been fitting Poisson and negative binomial distributions to the data using the goodfit function from the vcd library. I would also like to evaluate how well a quasi-poisson distribution fits the data. However, none of the potentially suitable functions I have identified (goodfit(vcd), fitdistr(MASS),

Hi! I am a psychologist who suspects that the only sensible way to analyse a particular data set is to use generalised linear mixed models. I am hoping that someone might be able to point me in the right direction to find some very practical hands on documentation that might be able to talk me through actually doing such an analysis? So far in my searches the most useful document I have turned

For example, subtracting 1:2 from the rows of a two-column matrix: > t(apply(matrix(1:6,ncol=2),MARGIN=1,function(y) y - 1:2)) [,1] [,2] [1,] 0 2 [2,] 1 3 [3,] 2 4 > sweep(matrix(1:6,ncol=2),MARGIN=2,1:2,FUN="-") [,1] [,2] [1,] 0 2 [2,] 1 3 [3,] 2 4 Is there a logic to this difference, or is it just a quirk of the history of these

Hi. As a follow-up on [1], I'd like to raise the question of whether it's practically possible to compile R to the Google Android mobile platform? The best way would probably be to use the Native Developer Kit [2], NDK, and in that way get a library. This could then be interfaced to by a Java-program. I know that several questions have been raised against the idea, e.g. with the input

Hi all: Here's a question about result of loglinear analysis. There're 2 factors:area and nation.The raw data is in the attachment. I fit the saturated model of loglinear with the command: glm_sat<-glm(fre~area*nation, family=poisson, data=data_Analysis) After that,I extract the coefficients: result_sat<-summary(glm_sat) result_coe<-result_sat$coefficients I find that all the

Hello, I have a question regarding estimation of the dispersion parameter (theta) for generalized linear models with the negative binomial error structure. As I understand, there are two main methods to fit glm's using the nb error structure in R: glm.nb() or glm() with the negative.binomial(theta) family. Both functions are implemented through the MASS library. Fitting the model using these

Hi, > update.packages(ask='graphics') gives me multiple warning (one per updated package?) similar to ... Warning: unable to move temporary installation '\\Server02\stats\R\library\2.10\file3de56e0d\locfit' to '\\Server02\stats\R\library\2.10\locfit' The final, updated, folders do not end up where they should be. I can move them 'by hand', but it is an

While putting my R code into functions, I've encountered a ddply function nesting issue and need a bit of advice on the proper way to fix it.? I've tried several approahces, but neither worked and I need to have the ability to include the "cut", "range", and "fullseq" methods within ddply.? (For a bit of that explanation refer to

Hi all, hope someone can help me out with this. Background Introduction I have a data set consisting of data collected from a questionnaire that I wish to validate. I have chosen to use confirmatory factor analysis to analyse this data set. Instrument The instrument consists of 11 subscales. There is a total of 68 items in the 11 subscales. Each item is scored on an integer scale between 1 to 4.

Hi, I've been following the list for a few months but this just my first post to the list. I've just started to use compiz again a few days ago and I have a few things (probably bugs) here that bugged me a little bit. I am not sure whether the bugs is in compiz or somewhere else so it would be nice if someone could confirm them. Here they are: - on default window decoration the maximize

I'm trying to apply a function defined in the VW R docs, that attemps to convert a data.table object to Vowpal Wabbit format. In the process i'm getting the error in printf mentioned in the subject. The original function is here: https://github.com/JohnLangford/vowpal_wabbit/blob/master/R/dt2vw.R Below there is a small example that reproduces the error. The function works great with

Hi all, I am getting the following error message: > mymodel = glm(response ~ . , family=binomial, data=C); Error: cannot allocate vector of size 734.2 Mb In addition: Warning messages: 1: In array(0, c(n, n), list(levs, levs)) : Reached total allocation of 1535Mb: see help(memory.size) 2: In array(0, c(n, n), list(levs, levs)) : Reached total allocation of 1535Mb: see help(memory.size) 3:

Has a test for bimodality been implemented in R? Thanks, John NIWA is the trading name of the National Institute of Water & Atmospheric Research Ltd.

Hi, I'm having a problem using sample() within a function. Basically I get an error reading: Error in sample(v, 1, prob = h) : non-positive probability Can anyone advise me as to the possible origin of this error? Here is my code #Discretised Gillespie algorithm function (From SMfSB, D.J. Wilkinson) gillespied=function (N, T=100, dt=1, ...) { tt=0 n=T%/%dt x=N$M