similar to: Inverse Gaussian Distribution

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

I have a set of data that includes various data columns. One if the survival time and another if a continuous variable of ages. I want to put the ages into intervals so that I can then perform the Kalpan Meier test. I am trying to use the following code to build a column with the age group numbers in agecatagory<-c( ) for (i in 1:137) { { if(age[i]<=46) {agecat[i]<-1} if(age[i]>46

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)+

Dear Sir/Madam, I am using R version 1.8.1. I am doing following tast: First generate 100 Gaussion(3,1) numbers, then write the likelihood function to estimate the parameters of Gaussian distribution by direct maximizing the likelihood function. My likelihood function is: >fn<-function(x) >(-50*log((sd(x))^2))-50*log(sqrt(2*pi))-(1/2*((mean(x))^2))*(sum((x-(mean(x))^2)) After I

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 trying to create a power curve to show how the power of a t-test varies depending on the mean. Any ideas how I should go about this? Louisa _________________________________________________________________ [[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

Hi, I am trying to write R code to produce a power curve to show how the power of a Wilcoxon-test varies depending on the mean, with data generated from a uniform distribution. Any ideas how I should go about this?Louisa _________________________________________________________________ Amazing prizes every hour with Live Search Big Snap [[alternative HTML version deleted]]

Dear, I'm very new to R Gui and I have to make an assignment on Gamma Regressions. Surfing on the web doesn't help me very much so i hope this forum may be a step forward. The question sounds as follows: The data set is in the library MASS first install library(MASS) then type data(mammals) attach(mammals) Assignment: Fit the gamma model and lognormal model for the mammals data.

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

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 am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponentiate the result." This is because R is unable to calculate high powers such as 545.501. As

I have read a frequency table in to R called "temp." I now want to create a histogram table from it, but I obviously first have to expand the data - to the sample size of 100. I want to use the command rep(), but I'm not sure how to go about it..I tried using the code: temp1<-rep(temp$Chest,100) hist(temp1) But this creates a v. odd histogram so I know it must be wrong!

Hi, I?m having much trouble on configuring Samba to work on an Active Directory environment. Using getent password I?m able to see AD?s users. wbinfo -u and wbinfo -g also work fine. When someone from a Windows try to access my Samba server, the smd password window is shown (I think that the autehntication would be transparent, wouldn't it ?), any password I provide is rejected: I tried AD

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