similar to: Odds Ratio and Logistic Regression

I'm resending this message because I did not include a subject line in my first posting. Apologies for the inconvenience! Tanja > Hello, > > I'm trying to fit a generalized linear mixed model to estimate diabetes prevalence at US county level. To do this I'm using the glmer() function in package lme4. I can fit relatively simple models (i.e. few covariates) but when

Hello, I'm trying to fit a generalized linear mixed model to estimate diabetes prevalence at US county level. To do this I'm using the glmer() function in package lme4. I can fit relatively simple models (i.e. few covariates) but when expanding the number of covariates I usually encounter the following error message. gm8 <-

Hi all, I've searched high and low on this and found nothing of help. I'm using v2.6.2 and trying to write a function that will count how many people from a dataset fall under a poverty line of 50% of the mean income. a9 is my 100-element vector of incomes. I want pa9 to be my vector that counts how many of these are classed as in poverty. My problem is that my command if(a9[i] <

Dear R users, Consider the first two columns of a data frame like this: > z[,1:2] x y 1 1 1 2 2 2 3 3 3 4 1 4 Imagine that y represents the times that the value x happens in a population. But z is not exactly a frequency table, because in z we have x=1 twice. So, the x=1 in the first line and the x=1 in the fourth are not the same, differing according to a third variable in

I'm on too many lists already, so haven't joined this one. If you have things that I really need to know, please reply or CC me off-list. I often work through masqueraded (NAT) links, and find that the following patch is able to keep an idle ssh link up where nothing else will. Please consider making it a part of the main ssh package.

Dear All, I am economist and working on poverty / income inequality. I need descriptive statitics like the ratio of education expentitures between different income quintiles where each household has a different weight. After a bit of google search I found 'Hmisc' and 'quantreg' libraries for weighted quantiles. The problem is that these packages give me only weighted quintiles;

Hello R-Helpers, I was looking into how to test whether the beta coefficient from a regression would be the same for two different groups contained in the dataset for the regression. When I put that question into google, AI returned a very nice looking answer (and a couple of variations on it). library(car) data <- data.frame(income = c(30, 45, 50, 25, 60, 55), education =

FAO Statistical Yearbook consultant The Food and Agriculture Organization (www.fao.org) is a United Nations agency that leads international efforts to reduce hunger in the world. This is achieved through three main goals: the eradication of hunger, food insecurity and malnutrition; the elimination of poverty and the driving forward of economic and social progress for all; and, the sustainable

My sincerest apologies, as this is a very elementary problem, but I have searched through archives and FAQs and on the web, and I am at the end of my own resources. I need to do analysis of a spatial point process occuring with four counties, using spatstat. I've had no trouble importing the shapefiles, but I can't seem to figure out what I need to do in order to use the boundary as an

Dear useRs, my last announcement is an update of the ineq package for measuring inequality, concentration and poverty. The current version is now 0.2-4. Thanks to suggestions from Rein Halbersma the Pen() function for plotting Pen's parade was improved and now allows for much more flexibility. See the help page for examples. Best wishes, Z ------------- Package: ineq Version: 0.2-4 Date:

Hello, I've started using R few months ago and I really like it. I need to estimate standard deviation of certain statistics (some measures of poverty). I found a really simple program, and I just need to check whether it's OK and really calculates what it's supposed to. Let's suppose e. g. head-count index (as one of the simplest measures of poverty) which is calculated as a

Dear useRs, my last announcement is an update of the ineq package for measuring inequality, concentration and poverty. The current version is now 0.2-4. Thanks to suggestions from Rein Halbersma the Pen() function for plotting Pen's parade was improved and now allows for much more flexibility. See the help page for examples. Best wishes, Z ------------- Package: ineq Version: 0.2-4 Date:

Dear Amy, I have the same question as Mark. Is it possible to work from abroad - remotely? Best regards, Martins On 6 January 2015 at 19:32, Mark LoPresti <marklopresti at yahoo.com> wrote: > Hi Amy, > > one key question: can this position work remotely? > Thank you, > Mark LoPresti > From: "Heyman, Amy (ESS)" <Amy.Heyman at fao.org> > To:

Hi, I'm sure this is simple, but I haven't been able to find this in TFM, say I have some data in R like this (pasted here: http://pastebin.com/raw.php?i=sjS9Zkup): > head(df) gender age smokes disease Y 1 female 65 ever control 0.18 2 female 77 never control 0.12 3 male 40 state1 0.11 4 female 67 ever control 0.20 5 male 63 ever state1 0.16

Hi, I'm moving (slowly) to R from STATA. I often have need to move through a set of tasks across a series of years. In this case, you can see that I'm mimicking -reshape- in STATA, but I'm less interested in the task than in programming R. library(foreign) mydata<-read.dta("z:\example.dta") for (y in 2000:2002) {

Hello list, I'm having some difficulty getting rJava to load in littler. Even after a R CMD javareconf and a reinstall of littler, I get this: jlaing at xenon:~$ r -e "require(rJava)" Loading required package: rJava Loading required package: methods Error : .onLoad failed in loadNamespace() for 'rJava', details: call: dyn.load(file, DLLpath = DLLpath, ...) error: unable

Hi all, Could somebody be so kind to explain to me what is the saturated model on which deviance and degrees of freedom are calculated when fitting a binomial glm? Everything makes sense if I fit the model using as response a vector of proportions or a two-column matrix. But when the response is a factor and counts are specified via the "weights" argument, I am kind of lost as far as

I don?t understand why you don?t include the full text of the error. ? David Sent from my iPhone > On Jan 19, 2025, at 10:00?AM, Sparks, John via R-help <r-help at r-project.org> wrote: > > ?Hello R-Helpers, > > I was looking into how to test whether the beta coefficient from a regression would be the same for two different groups contained in the dataset for the

I've been translating an economic model from Python into R, and I thought members of the list would like to see a presentation I've written about it. I've blogged this at http://www.j-paine.org/blog/2018/03/r-taxben-a-microsimulation-economic-model-in-r.html , and the presentation itself is a slideshow at http://www.j-paine.org/rtaxben/R/reveal/rtaxben.html . The slideshow is written

Dear All, I'am looking for examples showing that correlation does not imply causality, the targeted audience consists of undergraduate students (their first year at the university but in the BioMathStat track). All practicals are under R. I was able to extract this from R datasets: ### begin data(sunspots) data(lynx) spots <- window(sunspots, freq = 1, start = 1880, end = 1900) lnx <-