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]]

