similar to: calibration question

Good afternoon! As mentioned in the subject, my question regards more the methodological part that accompanies survey design and the statistical part that is involved. So, I have the following data: a<-data.frame (id_hh=c(1:5), strata=c(1,1,2,2,1), Nhstrata=c(100,100,200,200,100), Nrmemb=c(2,4,2,5,4)) a$ocmemb1<-c("wk","jl","st","jl","st")

Hi, I would like to do coloring of map regions based on the region values "weight". The approach I am taking is first to break regions into equal intervals, classIntervals(spdf$weight,4)$brks #4 intervals in this case and coloring all regions within the interval with the same color col = brewer.pal(4,"RdYlGn")) The max "weight" is as well the boundary of the

VR, James James T. Durant, MSPH, CIH Statistics and Data Analysis Team (proposed) Science Support Branch (proposed) Division of Community Health Investigations (proposed) US Agency for Toxic Substances and Disease Registry Mailing Address: 4770 Buford Highway, MS F-62 Atlanta, Georgia 30341 Telephone: 770-488-0668 (office) _____________________________________________ From: Lewis,

Hello List, Could someone tell why I can't install the Himsc and survey packages for R version 2.13.0 (2011-04-13)? What am I doing wrong here? Thanks, Pradip > install.packages ("Hmisc", dependencies=TRUE) --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository http://watson.nci.nih.gov/cran_mirror/bin/windows/contrib/2.13

Hi Everyone, I am a new R user and trying to run R jobs in batch mode. Robert Muenchen (2009), in his book "R for SAS and SPSS Users", has suggested writing a small batch file like mR.bat as shown below: "C:\Program File\R\R-2.10.0\bin\Rterm.exe" --no-restore --no-save < %1 > %1.Rout 2>& 1 Could anyone tell me in which directory or subdirectory I should

Hello, Using the svyplot () function, I have plotted four graphs that are saved in four different .png files. I am looking for examples how to redraw the same four graphs within grid viewports so that they stay together on a page. The goal is to create one .png file that will include all four graphs (2 rows, 2 columns). Any help would be appreciated. Thanks, Pradip Pradip K. Muhuri,

I am trying to import Indian National Sample Survey Data. It is ASCII flat and an example is below 001460306025001117101011117111201*01*000000211 270104070204093 002460306025001117101011117111201*02*0000008011116104910519572 022 2600 1150 30 55 170 4005 003460306025001117101011117111201*03* 0000111039204112 222 Previous post don't seem to provide much help. I have an document

Hi List, My goal is to force R not to print in scientific notation in the sixth column (rel_diff - for the p-value) of my data frame (not a matrix). I have used the format.pval () and printCoefmat () functions on the data frame. The R script is appended below. This issue is that use of the format.pval () and printCoefmat () functions on the data frame gives me the desired results, but coerces

Hello, I was trying to draw histograms of age at death and got the following 2 error messages: 1) Error in tapply(1:NROW(x), list(factor(strata)), function(index) { : arguments must have same length 2) Error in findInterval(mm[, i], gx) : 'vec' contains NAs In addition: Warning messages: 1: In min(x) : no non-missing arguments to min; returning Inf 2: In max(x) :

Hello List, My goal is to create a 5 category variable (p1_st_data$ob_mrj_cat), based on the p1_st_data$obt_mrj_p variable, using the following code for 50 States and District of Columbia (N=51). p1_st_data$ob_mrj_cat <- cut (p1_st_data$obt_mrj_p, quantile (p1_st_data$obt_mrj_p, (0:5/5), include.lowest=TRUE)) The issue is that, for Utah, I am getting an <NA> instead of (42,48.7] in

Hello List, My goal is to apply a user-defined function on several columns of a data frame. When testing the code on a reproducible example below, I get the following error message. > #now Write a new function using the above cut ()/quantile function to apply on different columns of the data frame > > CutQuintiles <- function(x) { + cut (test1$x,quantile (test1$x,

Dear R community, ? I am happy to announce the publication on CRAN of the revengc package: https://cran.r-project.org/web/packages/revengc/index.html The statistical package revengc was designed to reverse engineer censored, decoupledcensus data into a likely hhs x area uncensored contingency table for estimating interior residentialoccupancy.? ? For basic examples:

Hi there, I'm pleased to announce the availability (CVS only for now) of pymp3cut. pymp3cut is a GPLed Python command line tool designed to automatically cut very big mp3 files according to a timeline, and postprocess each subfile with any tool of your choice in a single run. It's useful if you stream very long events (10 hours or so) like congresses using IceCast and want to cut the

Dear R community, ? I am happy to announce the publication on CRAN of the revengc package: https://cran.r-project.org/web/packages/revengc/index.html The statistical package revengc was designed to reverse engineer censored, decoupledcensus data into a likely hhs x area uncensored contingency table for estimating interior residentialoccupancy.? ? For basic examples:

Hello, Although my R code for the svymean () and svyquantile () functions works fine, I am stuck with the svyby () and make.formula () functions. I got the following error messages. - Error: object of type 'closure' is not subsettable # svyby () - Error in xx[[1]] : subscript out of bounds # make.formula () A reproducible example is appended below. I would appreciate if

I have 2 server including the PDC. This is going to grow to 4 servers total. I have one PDC and i have 1000+ user accounts setup. Logins work fine from all computers. Currently I have 2 user groups one on each server in order to reduce the load because I am using old 300mhz computers. The problem that is see is that when i login and decide to set permissions on the samba "server" not the

I'm trying to create a function which will allow me to subset a data set based on values of various specified variables. I also want to then apply some other function(s) (e.g., summary). This is what I've tried so far.... > test.fx <- function(dta, expvar, expval) { + newdta <- subset(dta, eval(expvar)>expval) + summary(newdta$eval(expvar)) + } > >

Hello, I'd need to compute the calibration error of posterior class probabilities p(y|x) estimated by using rpart as classification tree. Namely, I train rpart on a dataset D and then use predict(... type="prob") to estimate p(y|x). I've found the possibility to do that in the ROCR package, but I cannot find a link to a paper/book which explains the details of the

Dear all, I used fit.mult.impute in Dr. Harrell's Design package to fit a cox ph regression model on five imputed datasets, where all missing predictors were filled by multiple imputation using R package Mice. Are there any functions able to do bootstrapping or cross-validation for the aggregated model? I tried function 'validate' and 'calibrate' in Design package, but

Is there a way to calibrate gamepads for wine? I have managed to get my PS3 sixaxis working, but when I play a game, everything works properly but the analogue stick. If I move it all the way to the right or down, it goes in the opposite direction, respectively. Anyone know how I might fix this?