similar to: How to script the script sample into script "OR", please advice

*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

If I wanted to fit a logit model and account for clustering of observations, I would do something like: library(Design) f <- lrm(Y1 ~ X1 + X2, x=TRUE, y=TRUE, data=d) g <- robcov(f, d$st.year) What would I do if I wanted to do the same thing with a probit model? ?robcov says the input model must come from the Design package, but the Design package appears not to do probit? Thanks very

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

Hi, I am trying to run a logistic regression to look at the risk of malaria infection in individuals. I want to account for intra household correlation and so want to include a household level random effect. I have been using the lmer command in lme4 package but am getting some strange results that are completely different to those I get using STATA. Can I just check that this is the correct

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

Dear all, Our pharmacometrics group is looking for bright and motivated candidates to join our modelling team within the Division of Clinical Pharmacology at the University of Cape Town. Our group has established PK/PD modelling expertise in the disease areas of tuberculosis, HIV, and malaria. It is active internationally and enjoys on-going collaborations and exchange programs with some of the

Dear all, We are looking for bright and motivated candidates to join our pharmacometric modelling team at the University of Cape Town. Pharmacometrics is an emerging discipline that uses mathematical and statistical tools to quantify drug, disease, and trial information to aid efficient drug development and/or regulatory decisions. It designs and applies mathematical models to describe the

es simple lo unico que tienes que hacer es instalar la libreria(xlsReadWrite) y direccionar y seguidamente especificar el numero de hoja de excel library(xlsReadWrite) dat.excel <- read.xls("d:/ruta del archivo/doc.xls", sheet=1) espero que te sirva ------------------------------------------ El día 18 de agosto de 2010 06:00, <r-help-es-request en r-project.org>