similar to: random effects using lmer

I have a dataframe imported from csv file below: Houseid,Personid,Tripid,taz 1,1,1,4 1,1,2,7 2,1,1,96 2,1,2,4 2,1,3,2 2,2,1,58 There are three groups identified based on the combination of first and second columns. How do I split this data frame? I tried aa <- split(inpfil, inpfil[,1:2]) but it has problems. Output desired is aa[1] Houseid,Personid,Tripid,taz 1,1,1,4 1,1,2,7 aa[2]

I needed to compute a complicated cross tabulation to show weighted means and standard deviations and the only method I could get that worked uses a series of nested for next loops. I know that there must be a better way to do so, but could use some assistance pointing the way. Here is my working, but inefficient script: library(Hmisc) rm(list=ls()) load('NHTS.Rdata') day.wt <-

I have cobbled together the following logic. It works but is very slow. I'm sure that there must be a better r-specific way to implement this kind of thing, but have been unable to find/understand one. Any help would be appreciated. hh.sub <- households[c("HOUSEID","HHFAMINC")] for (indx in 1:length(hh.sub$HOUSEID)) { if ((hh.sub$HHFAMINC[indx] == '01')

Hi all, can the glm package be used to estimate a logit model to panel data? I am asking this because stata has a standard logit model and then an xtlogit for longitudinal data. Is there something similar in R? Thank you jcm [[alternative HTML version deleted]]

Hi all, If you're looking at a Windows server to Linux & Samba migration (and let's face it, who isn't nowadays :-), the EU (god bless their little cotton socks :-) has published a wonderful migration how-to document here : http://europa.eu.int/ISPO/ida/jsps/index.jsp?fuseAction=showDocument&parent=news&documentID=1647 Cheers, Jeremy.

Dear User R It's been a pleasure talking with you. I am newcomer use R. Would you please help me how to translate the script below to "R" script? * Area under receiver operating characteristic (AU-ROC) predict r1m1p, p roctab malaria r1m1p, graph summary * Area under receiver operating characteristic (AU-ROC) curve predict r1m2p, mu roctab malaria r1m2p, graph summary

*Dear Bert, * Thank you very much for your feedback and the useful link https://rseek.org/ and https://www.r-bloggers.com/calculating-auc-the-area-under-a-roc-curve/. Actually, I want to know different performance between Stata and R, in multilevel logistic regression. For this purposes, I replicate ".do" file use Stata in

DEPARTMENT OF ZOOLOGY Web developer Grade 7: Salary range:£26,666 - £32,796 Three year post A web developer is required to join the Malaria Atlas Project (MAP, http://www.map.ox.ac.uk). The job will be to lead in the improvement of the on-line presence and thus public engagement activities of MAP. Specific responsibilities will be to implement an online database of epidemiological information

This is a bit obscure but caused me some grief... R1.5.1/MS Windows There seems to be a peculiarity in legend associated with the use of expression text. It seems as if expression text is used in (eg) an mtext call then if the call to legend also includes expression text and a subscript then the legend box and line spacing depends on the value of cex used in the previous call to mtext,

Hello R users I'm using MuMIn but for some reason I'm not getting the adjusted confidence interval and uncoditional SE whe I use model.avg(). I took into consideration the steps provided by Grueber et al (2011) Multimodel inference in ecology and evolution: challenges and solutions in JEB. I created a global model to see if malaria prevalence (binomial distribution) is related to any

Hi all, Quick question: What function can I use to draw a line in the margin of a plot? segments() and lines() both stop at the margin. In case the answer depends on exactly what I'm trying to do, see below. I'm using R v. 2.8.1 on Windows XP. Cheers, Alan I'm trying to make a horizontal barplot with a column of numbers on the right side. I'd like to put a line between the

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the

Hi R friends, I've been attempting to create plots with multiple alpha values using Cairo to save them on a windows (32b XP) platform as it doesn't support more than 3 alpha values. This worked well until I wanted a postscript file (unsupported) and as a attempted work around I installed RGtk2. So far so good, however now when I try to use a >CairoPDF("alpha.pdf", 6, 6,

A while ago (April 2002) there was a short thread on software for generalized linear mixed models. Since that time, has anyone written or found R code to use a glmm to analyze binomial data? I study CWD in white-tailed deer, and I'd like to do a similar analysis as Kleinschmidt et al. (2001, Am. J. Epidemiology 153: 1213-1221) used to assess control for spatial structure in malaria

I am a graduate student who's just starting to use more advanced statistics and is completely new to R. ?I'm looking for whether parasite prevalence varies with region. ?I'd like to try using just a subset of the data defined by one these covariates, sex, looking at data only from adults, without using it as a covariate in the analyses. ?I tried using subset=sex=='A' in the

Dear all, at the Division of Clinical Pharmacology at the University of Cape Town (South Africa) we are looking for highly motivated science/engineering PhD graduates to join our pharmacometrics group as a post-doctoral fellow. Pharmacometrics is an emerging discipline that applies mathematical modelling and simulation to pharmacology, aiding the interpretation of data from clinical trials,

Hi I hope someone can shed some light on this: For some reason when I read.table("bfx.txt") R decides to only give back the first character from each column in each row as one single column. Like this: V1 1 ÿþr 2 \n 3 r 4 1 5 0 6 A 7 G 8 \n 9 r 10 1 11 0 12 T 13 C 14 \n The data should be:

First, apologies for no example data but I don't think it's needed in this case, Q: Can (and if so how ) the line end style be changed for 'cloud' plots? I've tried par(lend=2), trellis.par.set(add.line = list(lend=2)) and much googling but to no avail Thanks in advance Dan P.S. the reason for this is that the round end looks bad at lwd=3 or more Daniel Alcock Malaria

Hi, It's unfortunate that with recent revisions of R 2.13 (this appeared in revision 54640, March 2), 'R CMD build' now removes empty dirs in the package. People might have good reasons for having empty dirs in their packages. For example, in Bioconductor, we have some tools to automatically generate annotation packages and those tools are implemented in software packages that use

Dear list, I will shortly have some data that contains numeric dates in the Persian / Jalali calendar format, which I would like to convert to gregorian. At the moment there doesn't seem to be a function for this in R, but it would be great if someone could come up with same - I would attempt it but the algorithm is very complex and this is also way beyond my fairly rudimentary knowledge of