similar to: Help with for/if loop

Hi everyone, my apologies in advance if I'm overlooking something simple in this question. I am trying to use R's survey package to make a direct method age-adjustment to some complex survey data. I have played with postStratify, calibrate, rake, and simply multiplying the base weights by the correct proportions - nothing seems to hit the published numbers on the nose. I am trying to

I am new to R. I have a data set like this (given below is a fictional dataset): AgeCat FINWT 1 98 2 62 1 75 3 39 4 28 2 47 2 66 4 83 1 19 3 50 I need to calculate the weighted distribution of the variable AgeCat. In SAS i can do: proc freq data=ageval; tables agecat; weight finwt; run; What or is there an equivalent in R? TIA, Krishnan -- Krishnan Viswanathan 1101 High Meadow Dr

Dear, I want to fit an inverse gaussion distribution to a data set. The predictor variables are gender, area and agecategory. For each of these variables I've defined a baseline e.g. #agecat: baseline is 3 data<-transform(data, agecat=C(factor(agecat,ordered=TRUE), contr.treatment(n=6,base=3))) The variable 'area' goes from A to F (6 areas: A,B,C,D,E,F) How can i

R 2.8.1 windows XP I am getting an error message that I don't understand when I try to run GLM. The error only occurs when I have all independent variables in the model. When I drop one independent variable, the model runs fine. Can anyone help me understand what the error means and how I can correct it? Thank you, John > fit11<-glm(AAMTCARE~BMI+BMIsq+SEX+jPHI+jMEDICAID+factor(AgeCat)+

I am doing group wise plots by using the following commands; it shows errors that I do not know how to fix it. Please help. xyplot(within.2.special.care ~ agecat| mco.cms.ind, neuro, panel = function(x,y){ + panel.grid() + panel.xyplot(x,y) + panel.loess(x,y, span =1)}) Error: couldn't find function "xyplot" [[alternative HTML version deleted]]

Hi! I am quiet new with R and I have some problems to perform a posthoc test with an lme model. My model is the following: >lme1<-lme(eexp~meal+time, random=~1|id,na.action=na.omit) and then i try to get a post hoc test: >summary(glht(lme1,linfct=mcp(meal="Tukey))) but I get a warning message: Erreur dans as.vector(x, mode) : argument 'mode' incorrect Thank you for your

I have a variable which is roughly age categories in decades. In the original data, it came in coded: > str(xxx) 'data.frame': 58271 obs. of 29 variables: $ issuecat : Factor w/ 5 levels "0 - 39","40 - 49",..: 1 1 1 1... snip I then defined issuecat as ordered: > xxx$issuecat<-as.ordered(xxx$issuecat) When I include issuecat in a glm model, the result

R Fans-- I am having problems with the following code. It worked under R 2.6.0 but not in 2.7.0. > library(Epi) > df <- read.table( "c:/Documents and Settings/Troy S/My Documents/debug_chisq_080613b.txt") > summary(df) cvd agecat Min. :0.0000 (0,40] :1 1st Qu.:0.0000 (40,60]:2 Median :0.0000 Mean :0.3333 3rd

Good point. Now this returns 0.04062184. Hmmm..... On Mon, Oct 2, 2017 at 6:30 PM, Hollie Johnson (PGR) < h.a.johnson at newcastle.ac.uk> wrote: > Hi Eric, > > > Thanks for having a look into this. I think you have a small typo... > > B <- matrix(x, nrow=3, byrow = TRUE) should read B <- matrix(y, nrow=3, > byrow = TRUE) > > > Regards, Hollie >

Hi, I'm working on some social networks and I managed to create the graphs with labels and edges weight, but I would also like to change the size of the vertices according to the age of the persons in the network and the shape according to the gender. Now for the age, I have people with ages between 20 and 78, and I would like to have 4 categories (sizes): 20-35, 36-50, 50-65, >65. I have

Due to missing data with some of my predictor variables I first do multiple imputation as follows: library(foreign) library(Amelia) library(norm) set.seed(666) M=10 impdat <-

Dear Colleagues, We are attempting to create trees using R with our Ruby on Rails application. However, we are running into a problem involving the creation of the graphic. We would like them to be in either jpg or png format so that users can save, but due to a lack of control over our sever we are unable to start X11 server. Is there a way to create these images without using X11 server?

Hi, I'm quite new at R and I haven't found the answer to my question anywhere on the net, so either it is trivial or not documented. So, bare with be. I am using the quadprog package and its solve.QP routine to solve and quadratic programming problem with inconsistent constraints, which obviously doesn't work since the constraint matrix doesn't have full rank. A way to solve this

Currently doing some work on local maxima on a random field and have encountered an issue with the Miwa algorithm used with the pmvnorm function in the mvtnorm R package. Based on recommendations by Mi et al., we ran the mvtnorm package using the Miwa algorithm, since we have a maximum of 4 dimensions with non-singular matrices. However, running the estimation procedure in this way, we obtained

Rather specialized. As this appears to be primarily a statistical, not an R programming question, you may do better posting on a statistical site like http://stats.stackexchange.com/ if you don't get a satisfactory reply here . Alternativey, if you think this is a package bug, perhaps contact the package maintainer directly, as (s)he may not monitor this list. -- Bert Bert Gunter

I am new in using R. I want to know all about building decision tree model in R. Few options which I searched are rpart and rattle to build a decision tree.Both the functions are giving me splits which are statistically appropriate. But I am not able to figure out how to change those splits as per my business requirement. for example : the automatic split of Age by using rattle is > 30 and

Hello R-experts, I have a set of data as follows: age time 1 28 1 2 53 2 3 53 3 4 36 4 5 54 4 6 46 4 7 45 5 8 31 6 9 53 7 10 35 7 11 62 8 12 19 8 13 43 2 14 51 3 15 45 0 16 48 2 17 49 3 18 57 2 19 45 3 20 27 10 21 33 12 22 29 14 23 46 16 24 45 19 25 43 2 26 28 6 27 50 7 28 28 4

On Sep 4, 2010, at 11:21 AM, Borja Ferrer wrote: > I've noticed this pattern happening with other operators aswell, but used xor in this example. As i said before, i tried with different register allocation orders, but it will produce always the same result. GCC is emitting longer code, but since LLVM is so nearer to the optimal code sequence i wanted to reach it. In LLVM, copies are

Hi Eric, Thanks for having a look into this. I think you have a small typo... B <- matrix(x, nrow=3, byrow = TRUE) should read B <- matrix(y, nrow=3, byrow = TRUE) Regards, Hollie ________________________________ From: R-help <r-help-bounces at r-project.org> on behalf of Eric Berger <ericjberger at gmail.com> Sent: 02 October 2017 16:16:13 To: Bert Gunter Cc: r-help at

hi, i am using par(mrow=c(6,6)) function to get 6x6 plots on one screen. the problem that i am having is that the axis tick labels are far away from the ticks and going into previous plots (see attached figure). i need to know how can i reduce the distance between the ticks and their values (y axis values).the part of the code that i am using (after reading in the data) to create the top row is: