Mark Donoghoe
2008-May-12 07:11 UTC
[R] Standard error of combination of parameter estimates
Hi, Is there a simple command for computing the standard error of a combination of parameter estimates in a GLM? For example: riskdata$age1 <- riskdata$age riskdata$age2 <- ifelse(riskdata$age<67,0,riskdata$age-67) model <- glm(death~age1+age2+ldl, data=riskdata,family=binomial(link=logit)) And we want to find the standard error of the partial linear predictor for the combination of the two age parameters, age1+age2, at some age X, without needing to separately use the covariance matrix? Thanks, Mark Donoghoe ##################################################################################### This e-mail message has been scanned for Viruses and Content and cleared by MailMarshal ##################################################################################### #################################################################################################################### IMPORTANT NOTICE: This e-mail and any attachment to it are intended only to be read or used by the named addressee. It is confidential and may contain legally privileged information. No confidentiality or privilege is waived or lost by any mistaken transmission to you. The CTC is not responsible for any unauthorised alterations to this e-mail or attachment to it. Views expressed in this message are those of the individual sender, and are not necessarily the views of the CTC. If you receive this e-mail in error, please immediately delete it and notify the sender. You must not disclose, copy or use any part of this e-mail if you are not the intended recipient. ##################################################################################################################### [[alternative HTML version deleted]]
Søren Højsgaard
2008-May-12 09:28 UTC
[R] Standard error of combination of parameter estimates
You can use the esticon function in the doBy package. Regards S?ren ________________________________ Fra: r-help-bounces at r-project.org p? vegne af Mark Donoghoe Sendt: ma 12-05-2008 09:11 Til: r-help at r-project.org Emne: [R] Standard error of combination of parameter estimates Hi, Is there a simple command for computing the standard error of a combination of parameter estimates in a GLM? For example: riskdata$age1 <- riskdata$age riskdata$age2 <- ifelse(riskdata$age<67,0,riskdata$age-67) model <- glm(death~age1+age2+ldl, data=riskdata,family=binomial(link=logit)) And we want to find the standard error of the partial linear predictor for the combination of the two age parameters, age1+age2, at some age X, without needing to separately use the covariance matrix? Thanks, Mark Donoghoe ##################################################################################### This e-mail message has been scanned for Viruses and Content and cleared by MailMarshal ##################################################################################### #################################################################################################################### IMPORTANT NOTICE: This e-mail and any attachment to it are intended only to be read or used by the named addressee. It is confidential and may contain legally privileged information. No confidentiality or privilege is waived or lost by any mistaken transmission to you. The CTC is not responsible for any unauthorised alterations to this e-mail or attachment to it. Views expressed in this message are those of the individual sender, and are not necessarily the views of the CTC. If you receive this e-mail in error, please immediately delete it and notify the sender. You must not disclose, copy or use any part of this e-mail if you are not the intended recipient. ##################################################################################################################### [[alternative HTML version deleted]] ______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.