similar to: survey package question

svymean w/ svyby is working for me... > svyby(~visitcnt, ~agegrp3.f, svymean, design=svydes) agegrp3.f visitcnt se.visitcnt 18-44 18-44 8.755552 0.4953235 45-64 45-64 10.131555 0.5347806 65+ 65+ 9.588802 0.4323629 svyquantile is working for me... > svyquantile(~visitcnt, quantiles=c(.25, .5, .75), ties="rounded", design=svydes) 0.25

I am having trouble understanding (i.e. getting) confidence intervals from the survey package. I am using R version 2.8.1 (2008-12-22) and survey package (3.11-2) on FC7 linux. To simplify my question I use an example from that package: R> data(api) R> dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc) R> (tst <- svyby(~api99, ~stype,

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

Good morning! I have 4 questions which trouble me: 1. I want to test the hypothesis that the 2 proportions (the mean of a binomial) which come from 2 different samples are equal. I want to use the following function z= (p1-p2)/ sqrt((p1(1-p1)/n1)+(p2(1-p2)/n2)) which is one of the standard formulas for this case. Is there such a function in R? p1=the proportion from the first sample n1=the

Hello, Is there a way to get complex sample variances in the survey package on summary statistics other than means? If not, can they be added to a future version? It would be be great to have them on totals, quantiles, ratios, and tables (eg row percent, columns percent, etc). Thanks. Fred --------------------------------- [[alternative HTML version deleted]]

Hi, I want to view the code for the function confint in the stats package. Typing just confint doesn't work. How can I view the code? Thanks, Walt ________________________ Walter R. Paczkowski, Ph.D. Data Analytics Corp. 44 Hamilton Lane Plainsboro, NJ 08536 ________________________ (V) 609-936-8999 (F) 609-936-3733 walt at dataanalyticscorp.com www.dataanalyticscorp.com

Dear List: I have a data frame prepared in the couting process style for including a binary time-dependent covariate. The first few rows look like this. PtNo Start End Status Imp 1 1 0 608.0 0 0 2 2 0 513.0 0 0 3 2 513 887.0 0 1 4 3 0 57.0 0 0 5 3 57 604.0 0 1 6 4 0 150.0 1 0 The outcome

Hello everybody I'm having some difficulties to design the decision tree algorithm J48. I am using the following code and when I run it gives me the following message plot(m1) Error in plot.Weka_tree(m1) : Plotting of trees with multi-way splits is currently not implemented. #The code library(RWeka) library(randomForest) library(party) if(require(mlbench, quietly = TRUE) &&

Hi all, Appologies for the rather basic IO question but I am rather new to R... Migrated from IDL/Matlab recently. I have a rather simple Fortran control file (sigh...) that I am trying to parse and read using R. My problem is that the file's format is somewhat flexible. Imagine: --- 1> 39 1901 2> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 3> 22 24 26 28 30 32 34 36

Hello all, here's a real-world example: I'm measuring a quantity (d) at five sites (site1 thru site5) on a silicon wafer. There is a clear site-dependence of the measured value. To find out if this is a measurement artifact I measured the wafer four times: twice in the normal position (posN), and twice rotated by 180 degrees (posR). My data looks like this (full, self-contained code at

Hi, I have a question about estimating IRT models by using nlme, not just rasch model, but also other models. Behavior Research Methods <http://www.springerlink.com/content/1554-351x/> Volume 37, Number 2 <http://www.springerlink.com/content/1554-351x/37/2/>, 202-218, DOI: 10.3758/BF03192688 Using SAS PROC NLMIXED to fit item response theory models (2005). Ching-Fan