similar to: clusters in zero-inflated negative binomial models

Hello! I am doing an analysis on a questionnaire of hunters taken in 4 different districts of some mysterious foreign country. The aim of the study was to gather info on the factors that determine the hunting success of a peculiarly beautiful bird in that area. All variables are factors, i.e. they are variables such as "Use of Guns - yes / no", "Use of Dogs - yes / no" and

For a sanitation project in Bangladesh I need to design a three stage sample survey to be representative of around 40 million people. I find myself suddenly with several challenges with which I struggle and would be gratefully for any help. As the questions are linked I kept them together rather than creating multiple posts 1) SURVEY DESIGN So far I designed mainly two stage cluster surveys but

Hello, I apologize in advance as this example seems really elementary. Below I have created a simple scatterplot with lines. I would like to label each line with the name of the village, instead of using a legend. Can someone please let me know how to do it. Thanks in advance -- John Pop=c(1406,1862,2333,2664,2685,2521,2726,0,0,0,0,214,289,365,900,952,1002,1205,1041,942,1004,

Hello, I am new to R. I would like others to explain to me how to add absolute values inside the individual stacked bars in a consistent way using the basic R plotting function (R base). I tried to plot a stacked bar graph using R base but the values appear in an inconsistent/illogical way in such a way that its supposed to be 100% for each village but they don't sum up to 100%. Here is the

Hallo to everybody, I would like to perform an analysis but I don't know how to proceed and whether R packages are available for my purpose or not. Therefore I'm here to request your support. *The idea is the following:* I noticed that the names of the towns and villages in northern Italy most of the time sound differently from names of cities based on southern Italy. Just to give you an

Dear All, I am using the function vuong from pscl package to compare 2 non nested models NB1 (negative binomial I ) and Zero-inflated model. NB1 <-  glm(, , family = quasipoisson), it is an object of class: "glm" "lm" zinb <- zeroinfl( dist = "negbin") is an object of class: "zeroinfl"   when applying vuong function I get the following: vuong(NB1,

I would like to test in R what regression fits my data best. My dependent variable is a count, and has a lot of zeros. And I would need some help to determine what model and family to use (poisson or quasipoisson, or zero-inflated poisson regression), and how to test the assumptions. 1) Poisson Regression: as far as I understand, the strong assumption is that dependent variable mean = variance.

I want to compare the poisson and the zero-inflated poisson distribution on describing the data. So, after using the GLM and the ZEROINFL function, I used the Voung function to compare them. Here is my code: library(pscl) glm1 <- glm(nmer9_1[, 1] ~ 1, family = poisson) zip <- zeroinfl(nmer9_1[, 1] ~ 1) vuong(glm1, zip) However, R returns the following error message: "Error: cannot

Hi, im new to RoR so sorry if the questions sounds stupid to you. Im trying to build a tool for a browsergame. The game allows access to its database and the db looks like this: table villages: `id` int(9) unsigned NOT NULL default ''0'', `x` smallint(3) NOT NULL default ''0'', `y` smallint(3) NOT NULL default ''0'', `tid` tinyint(1)

Hi all, I wrote a simple code to set up the "pit_id" from a village table from 1 to 2, from 2 to 3 and so on. It works fine without .delay . But i need it so whats wrong with that code? #village_controller def post_info @village = Village.find(params[:id]) @village.delay.upgrade flash[:notice] = "Upgrading" redirect_to village_url(@village) end

Good Morning. Below I like statement like j<-grep(".r\\b",colnames(mydata),value=TRUE); j with the \\b option which I read long time ago which Ive found useful. Are there more or these options, other than ? grep? Thanks. dstat is just my own descriptive routine. > x ?[1] "age"????????? "sleep"??????? "primary"????? "middle" ?[5]

Hi, I have two questions. First, I'd like to compare a ZINB model to a negativ binomial model with the Vuong test, but I can't find how to performe it from the zicount package. Does a programm exist to do it ? Second, I'd like to know in which cases we have to use a double hurdle model instead of a zero inflated model. Many thanks, St??phanie Payet REES France R??seau

?s 02:10 de 02/08/2024, Steven Yen escreveu: > Good Morning. Below I like statement like > > j<-grep(".r\\b",colnames(mydata),value=TRUE); j > > with the \\b option which I read long time ago which Ive found useful. > > Are there more or these options, other than ? grep? Thanks. > > dstat is just my own descriptive routine. > > > x > ?[1]

..::: Drunk Young :::... http://210.101.95.254/drunk-young ^^^^ ^^^^ Click here to Join Drunk Young will take you to the new level of excitement. It is a guarantee of Your pleasure, because we choose models very seriously & carefully, depending on you offers and suggestions. 'Drunk Young' is regularly updates to bring you more and more sweet nymphets 100%

Dear All, I am using the function vuong of the package pscl to compare 2 non nested glm models with a numeric response. I did the following m1<-glm(y ~x ,data=xxx) m2<-glm(y ~z , data=xxx) When calling the vuong function I get the following message: > vuong(m1,m2) Error in predprob.glm(m1) : your object of class glm is unsupported by predprob.glmyour object of class lm is unsupported

Hello, I've been using a pre-release version of R v 2.8.0 for Windows for the last couple months. I think that there have been consistent problems with subsetting data sets, but I had usually been able to find work-arounds or was unable to confirm this as a bug. I think now I have, and would love advice on what to do if I've made some error. The data set in question ("c") has

Hi, I've been using Debian for a few years, and there was one nifty little app that made installing and updating so much easier: apt-proxy. Most of the time, I'm taking care of small LANs with an average of five client PCs. But this is a very remote place in South France, so most villages only have 512 kbps DSL. One major update for openoffice.org-*, and I have to wait the whole day

Dear All, I am using the function vuong of the package pscl to compare 2 non nested glm models with a numeric response. I did the following m1<-glm(y ~x ,data=xxx) m2<-glm(y ~z , data=xxx) When calling the vuong function I get the following message: > vuong(m1,m2) Error in predprob.glm(m1) : your object of class glm is unsupported by predprob.glmyour object of class lm is unsupported

Dear all, We are doing some research about the time series analysis of NDVI, and we found the NDVITS package which is a very great tool. Unfortunately when we run it, after TimeSeriesAnalysis it asks to enter "Village or Country". library("ndvits", lib.loc="/home/vahe/R/i686-pc-linux-gnu-library/2.15")

Hi R experts: I am new to mixed model commodity. I am tryping to specify a model using lmer in lme4 package. I am not sure if I am doing right, so I need your help....please...... Treatment / factor structure Year: level 1:3, the whole the experiment was repeated in three years, random factor village: level 1:2 # the level is much higher just three are shown as example, random factor Farm