search for: ahrq

I need some guidance from someone who is familiar/has some experience with the survey package. The data that I am using is from the Medical Expenditure Panel Survey (www.meps.ahrq.gov <http://www.meps.ahrq.gov/> ). The STRATA and PSU variables are varstr01 and varpsu01 respectively. When I try to specify them with the svydesign function I get an error message. An excerpt of my session is as follows... library(foreign) > h60 <- read.xport("h:\\meps\\t...

...her programs. Typically, I provide relatively thorough project documentation in my program files so that if someone else inherits them they can pick-up on them quickly. Thanks. Marc Marc W. Zodet, MS Health Statistician Agency for Healthcare Research and Quality Rockville, Maryland 20841 mzodet@ahrq.gov [[alternate HTML version deleted]]

...in the same tables....anything similar to crosstabs in S? Marc W. Zodet, MS Health Statistician Center for Cost and Financing Studies Division of Statistical Research and Methods 2101 East Jefferson Street, Suite 500 Rockville, Maryland 20852 Phone: 301-594-7072 Fax: 301-594-2166 E-mail: mzodet@ahrq.gov <mailto:mzodet@ahrq.gov> [[alternate HTML version deleted]]

...det, MS Health Statistician Agency for Healthcare Research and Quality Center for Financing, Access, and Cost Trends Division of Statistical Research and Methods 540 Gaither Road, Room 5058 Rockville, Maryland 20850 Phone: 301-427-1563 FAX: 301-427-1276 E-mail: <mailto:mzodet@ahrq.gov> mzodet@ahrq.gov [[alternative HTML version deleted]]

If there are any other users who use AHRQ's SAS code comoanaly2010 and comformat2010 to create comorbidity variables, I thought you might be interested in the following PRELIM code we wrote to mimic its functionality in R. It seems to yield similar results, but may contain errors. Please feel free to comment (kindly) or enhance. I'...

If anyone wants to see a prime example of how inefficient it is to program in SAS, take a look at the SAS programs provided by the US Agency for Healthcare Research and Quality for risk adjusting and reporting for hospital outcomes at http://www.qualityindicators.ahrq.gov/software.htm . The PSSASP3.SAS program is a prime example. Look at how you do a vector product in the SAS macro language to evaluate predictions from a logistic regression model. I estimate that using R would easily cut the programming time of this set of programs by a factor of 4. Fran...

...erent and I'm trying to understand the differences between the two packages (approaches) so I can explain why the final terms might be different. Thanks for your help. Marc Marc W. Zodet, MS Health Statistician Agency for Healthcare Research and Quality E-mail: <mailto:mzodet@ahrq.gov> mzodet@ahrq.gov [[alternative HTML version deleted]]

How do I go about generating a WEIGHTED mean (and standard error) of a variable (e.g., expenditures) for each level of a categorical variable (e.g., geographic region)? I'm looking for something comparable to PROC MEANS in SAS with both a class and weight statement. Thanks. Marc [[alternative HTML version deleted]]

...better to attach newdta? Thanks in advance for any guidance. Marc Marc W. Zodet, MS Senior Health Statistician Agency for Healthcare Research and Quality Center for Financing, Access, and Cost Trends 301-427-1563 (Telephone) 301-427-1276 (Fax) marc.zodet@ahrq.hhs.gov [[alternative HTML version deleted]]

What is the best way to sort a dataframe? For example, how would I go about sorting a dataframe (with variables V1-V5) by ascending V1, V2 and descending V3 while retaining V4 and V5. Also, is there a relatively easy way by which to re-order my columns? Thank you for your time. Marc Marc Zodet [[alternative HTML version deleted]]

I'm attempting to select a model using stepAIC. I want to use a forward selection procedure. I have specified a "scope" option, but must not be understanding how this works. My results indicate that the procedure begins and ends with the "full" model (i.e., all 17 independent variables)...not what I expected. Could someone please point out what I'm not

>I am interested in R as an alternative for a statistical tool >at our firm. Ditto... I have recently moved to this agency from a company where I had access to Splus. There is also a coworker here who had used Splus at a previous employer. We both would like some access to the S language. We are considering either begging loud and long to try to get the agency to purchase two copies of

I'm attempting to use stepAIC to select a model through a forward procedure. I want to consider up to all 3-way interactions. I've attempted to use the following code: m2.Fwd3way <- stepAIC(m1.Ionly, direction="forward", scope=list(upper=~(var1 + var2 + var3 + var4)^3, lower=~1)) When I submit this the trace indicates that only 2-way

How do I rotate 3D plots/surfaces generated by either cloud or wireframe? Thanks. Marc Marc W. Zodet, MS Health Statistician Agency for Healthcare Research and Quality [[alternative HTML version deleted]]

All: Is there any way to run an R script file (i.e., *.R) from the command prompt in the console window. Ultimately, I'm looking to put such code in a script file so that it can set off other R scripts/programs as needed. Thanks. Marc [[alternative HTML version deleted]]

This is a relatively recent article that is somewhat accessible. Jensen, D. R., and Solomon, Herbert (1994), "Approximations to joint distributions of definite quadratic forms", Journal of the American Statistical Association, 89 , 480-486 It has references to previous work. I also have an old paper that is so old I can't tell what journal it came out of:( Grad, Arthur and Solomon,

I'm trying to simulate trend data over a five year period. I want different trend profiles...the simplest being a linear trend. I've been using the following code: patBdta1 <- NULL for(i in 1:100) patBdta1 <- rbind(patBdta1,c(yr1= mean(rbinom(50,1,.50)), yr2 =mean(rbinom(50,1,.51)), yr3 =mean(rbinom(50,1,.52)),

...gt;>>> program in SAS, take a look at the SAS programs provided by the US >>>>>>> Agency for Healthcare Research and Quality for risk adjusting and >>>>>>> reporting for hospital outcomes at >>>>>>> http://www.qualityindicators.ahrq.gov/software.htm . The >>>>>>> >>>>>> PSSASP3.SAS >>>>> >>>>>> program is a prime example. Look at how you do a vector product in >>>>>>> >>>>>> the >>>>>> >>>...

Dear R users, I?m a beginner user of R and I?ve a problem with permutations that I don?t know how to solve. I?ve 12 elements in blocks of 3 elements and I want only to make permutations inter-blocks (no intra-blocks) (sorry if the terminology is not accurate), something similar to: 1 2 3 | 4 5 6 | 7 8 9 | 10 11 12 ----------1st permutation 1 3 2 | 4 5 6 | 7 8 9 | 10 11 12 NO - - 3 2 1

< Does anyone know why they're called random deviates, as opposed to random numbers?> Others will probably give you some technical reason about random numbers can be considered as random deviates from a mean (I think at least the 1875 Galton paper at http://www.mugu.com/galton/ uses similar terminology (I'm not claiming this is the earliest use - just the easiest to access at the