[R] Instantaneous hazard of an event i.e. death

From: Andrew McFadden <Andrew.McFadden_at_maf.govt.nz>
Date: Thu, 17 Apr 2008 11:09:11 +1200


Hi all

I am trying to calculate the instantaneous hazard (rather than cumulative hazard) for an event "bird death".

Some sample data is as follows:

year<-c(rep(2003,10),rep(2004,10))
time<-c(5,10,5,15,20,10,15,15,25,25,5,15,15,5,25,10,25,15,35,25)
event<-c(1,0,1,1,0,1,1,0,1,0,1,0,0,1,0,0,0,0,1,0)

Ideally I would like to know if the hazard is greater at specific time periods for birds dying within each year ie plotting confidence limits around the hazard for each year.

I would also like to present this data graphically.

I think in theory I should be able to use: the "hazard.ratio.plot" but lack the experience to know how to code this

Anyone be able to help?

Regards

Andrew

Andrew McFadden MVS BVSc
Incursion Investigator
Investigation & Diagnostic Centres - Wallaceville Biosecurity New Zealand Ministry of Agriculture and Forestry

Phone 04 894 5600 Fax 04 894 4973 Mobile 029 894 5611 Postal address: Investigation and Diagnostic Centre- Wallaceville Box 40742 Ward St Upper Hutt

########################################################################
This email message and any attachment(s) is intended solely for the addressee(s) named above. The information it contains is confidential and may be legally privileged. Unauthorised use of the message, or the information it contains, may be unlawful. If you have received this message by mistake please call the sender immediately on 64 4 8940100 or notify us by return email and erase the original message and attachments. Thank you.

The Ministry of Agriculture and Forestry accepts no responsibility for changes made to this email or to any attachments after transmission from the office.
########################################################################

        [[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. Received on Thu 17 Apr 2008 - 00:45:34 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 17 Apr 2008 - 01:30:29 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive