Or just approximately? –Ben Bolker Nov 29 '11 at 2:39 A bit of a stab in the dark, but what does dividing by sqrt(data2$REPLICATE_VAR) (rather than sqrt(sum(data2$REPLICATE_VAR)) do ... Sander Greenland , Model-based Estimation of Relative Risks and Other Epidemiologic Measures in Studies of Common Outcomes and in Case-Control Studies, American Journal of Epidemiology 2004;160:301-305 4. The model runs when I say type=ind which, if I understand correctly, means that the repeated measure are not correlated. This resolved the issue.Hopefully this will be helpful to others who might encounter similar issues as this type of error seems to have relatively little coverage.

Similarly, the use of the identity link does not ensure valid mean values for response distributions having bounded ranges of their means such the binomial, Poisson, gamma, and others. Can I release a pattern without releasing the whole held expression? Estimating the Relative Risk in Cohort Studies and Clinical Trials of Common Outcomes. The request for odds ratios can be dealt with once we figure out how to get the model to stop throwing errors.

Not the answer you're looking for? How exactly does mining software work? The system returned: (22) Invalid argument The remote host or network may be down. I waited as long as 20 minutes but nothing.

Here's my syntax for calculating the rate of hospitalisation. more stack exchange communities company blog Stack Exchange Inbox Reputation and Badges sign up log in tour help Tour Start here for a quick overview of the site Help Center Detailed Are the cells with sparse counts also typified by extreme values of the continuous covariates?2. I would think that either different response variable should be used or different DIST and LINK...

Both methods use proc genmod. Using a Poisson model without robust error variances will result in a confidence interval that is too wide. Parameter Information Parameter Effect carrot Prm1 Intercept Prm2 carrot 0 Prm3 carrot 1 Criteria For Assessing Goodness Of Fit Criterion DF Value Value/DF Deviance 98 132.3665 1.3507 Scaled Deviance 98 132.3665 Mitt kontoSÃ¶kMapsYouTubePlayNyheterGmailDriveKalenderGoogle+Ã–versÃ¤ttFotonMerWalletDokumentBloggerKontakterHangoutsÃ„nnu mer frÃ¥n GoogleLogga inDolda fÃ¤ltSÃ¶k efter grupper eller meddelanden ERROR The requested URL could not be retrieved The following error was encountered while trying to retrieve the URL: http://0.0.0.9/

PROC GENMOD DATA = temp1; ODS OUTPUT ParameterEstimates = results; CLASS bmi_cat id; MODEL hosp_flag = bmi_cat age_year / DIST = poisson LINK = log OFFSET = logpyr TYPE3 SCALE = Analysis Of Initial Parameter Estimates Standard Wald 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept 1 -0.8873 0.2182 -1.3150 -0.4596 16.53 <.0001 carrot 0 1 0.4612 There is no single cause for this error message that can be readily seen by examining the data or model specification. Here it is specified as log instead of logit: proc genmod data = eyestudy descending; class carrot; model lenses = carrot/ dist = binomial link = log; estimate 'Beta' carrot 1

Got a question you need answered quickly? Any chance you can post your data somewhere, or we can work with a subset? Add your answer Question followers (9) Robert A Yaffee New York University Stefan K Lhachimi Leibniz Institute for Prevention Research and Epidemiology â€“ BIPS Joachim Rosenbauer German Diabetes All that needs to be changed is the link function between the covariate(s) and outcome.

The number of rows would be 177,050,435, taken from the sum of REPLICATE_VAR. Apparently this would be the number of males aged 9 years and up in the US population. –Michelle Nov 29 '11 at 1:26 add a comment| Your Answer draft saved If they are approximately equal, change to a Poisson distribution. Check to see if the sample mean is roughly equal to the sample variance for the cells in your cross-tabulation.

Here gender and latitude will be added to the model: proc genmod data = eyestudy; class carrot gender id; model lenses = carrot gender latitude/ dist = poisson link = log; Analysis Of GEE Parameter Estimates Empirical Standard Error Estimates Standard 95% Confidence Parameter Estimate Error Limits Z Pr > |Z| Intercept -0.6521 0.4904 -1.6134 0.3091 -1.33 0.1836 carrot 0 0.4832 0.1954 Analysis Of Parameter Estimates Standard Wald 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept 1 -0.3567 0.2845 -0.9143 0.2010 1.57 0.2100 carrot 0 1 0.9892 0.4136 Scale 0 1.0000 0.0000 1.0000 1.0000 NOTE: The scale parameter was held fixed.

Please try the request again. Iteration will be terminated.ERROR: Error in parameter estimate covariance computation.ERROR: Error in estimation routin I f I run the analysis without the modification of: repeated subject=id/type=ind; or repeated subject=id/type=unstr; it works The exact message I am getting is NOTE: Class levels for some variables were not printed due to excessive size.NOTE: PROC GENMOD is modeling the probability that c_othvst='1'.NOTE: Algorithm converged.ERROR: Error Generated Mon, 17 Oct 2016 03:30:37 GMT by s_wx1094 (squid/3.5.20)

What kind of distribution is this? Adjusting the relative risk for continuous or categorical covariates Adjusting the RR for other predictors or potential confounders is simply done by adding them to the model statement as you would The control patients are ones initiated on extended release methylphenidate and the cases are those on immediate release methylphenidate. The outcome generated is called lenses, to indicate if the hypothetical study participants require corrective lenses by the time they are 30 years old.

latitude 1 -0.0100 0.0190 -0.0472 0.0272 0.28 0.5980 Scale 0 1.0000 0.0000 1.0000 1.0000 NOTE: The scale parameter was held fixed. References 1. Type:Usage NotePriority:lowTopic:Analytics ==> Analysis of VarianceAnalytics ==> Categorical Data AnalysisSAS Reference ==> Procedures ==> GENMODAnalytics ==> RegressionDate Modified:2007-10-11 16:46:29Date Created:2002-12-16 10:56:41 This content is presented in an iframe, which your browser If you had been doing this in GLIMMIX, the error would have been "Infinite likelihood in iteration 1"--so it is good to know the equivalent wording in GENMOD.

I bolded a couple of changes that might help, but I would not be surprised if this took several hours, and you got the message that it had not converged. If you're off by a factor of exactly 100 then I don't think that dividing by individual REPLICATE_VAR factor is going to be the way to go ... –Ben Bolker Nov Unfortunately, I think the errors that are currently occurring will still occur under these options, so perhaps some others can help out on this.Steve Denham Message 4 of 18 (1,171 Views) Were you possibly considering the number of visits as an ordinal response?

I am new to proc genmod and am having an issue with this. It started out asking about calculation of sample size, but has morphed into a discussion of analysis methods, and I think two of the recent posts (by Steve Simon and Mark Contrast Estimate Results Standard Chi- Label Estimate Error Alpha Confidence Limits Square Pr > ChiSq Beta 0.9892 0.4136 0.05 0.1786 1.7997 5.72 0.0168 Exp(Beta) 2.6891 1.1121 0.05 1.1956 6.0481 The estimate Also notice that the distribution has been changed to Poisson, but the link function remains log.

If excluding the propensity variables does not work, then we are dealing with a whole other set of problems.Steve DenhamSteve Denham Message 6 of 18 (1,171 Views) Reply 0 Likes Pooja Message 1 of 2 (433 Views) Reply 0 Likes statsplank Frequent Contributor Posts: 77 Re: Proc genmod - how to resolve error messages Options Mark as New Bookmark Subscribe Subscribe to For more enlightenment I'd simplify the problem to one variable and see what happens. –Spacedman Nov 28 '11 at 0:01 add a comment| 2 Answers 2 active oldest votes up vote McNutt LA, Wu C, Xue X, Hafner JP.

Filter columns basing on the value of other column Why does this execution plan have Compute Scalars? Spiegelman, D. share|improve this answer answered Nov 28 '11 at 0:14 Spacedman 53.7k561129 I have tried to use expand.dft() from the vcdExtra package, having removed all objects barring the one being Now this can be contrasted with the two methods of calculating the RR described below.

If this is incorrect, what does appear in the output?5. Does chilli get milder with cooking? How can I make LaTeX break the word at the end of line more beautiful? One estimates the RR with a log-binomial regression model, and the other uses a Poisson regression model with a robust error variance.

Is there a way to overcome this issue?Thanks,Pooja Message 9 of 18 (1,171 Views) Reply 0 Likes SteveDenham Super User Posts: 2,546 Re: Proc genmod error Options Mark as New Bookmark For example, in a binomial response (logistic) model, this error appears if the estimated binomial mean, p, is less than zero or greater than one.