search for: death

I came across this one: http://www.nysun.com/article/32787 which says that the violent death rate in Iraq (which presumably includes violent deaths from the war) is lower than the violent death rate in major American cities. Does anyone have any insights from statistics on how to interpret this?

...the calc is to be made but I can't figure out how to tell R to take the current value and divide it by the next value. The formula is P0=L1/LO The following is an example of the lifetable I'm constructing: example=as.vector(c(8,2,1,5,6,7,7,0,8,10,13,8,11,11,11,2,7,1,5,6,8,6)) #prime k=1 #Deaths# deaths=numeric(k) for(k in 0:max(example)) { deaths[k]=sum(example==k)} #adjust for no zero!!# deaths=c(0,deaths) #Alive, Kx# alive=sum(deaths)-cumsum(deaths) #Day, or age class,x# day=seq(from=0,to=length(deaths)-1) #mortality, lx# lx=alive/sum(deaths) #proportion surviving to next interval, p...

...urceCancel [==========] Running 1 test from 1 test case. [----------] Global test environment set-up. [----------] 1 test from AddressSanitizerMac [ RUN ] AddressSanitizerMac.GCDSourceCancel /sw/src/fink.build/llvm32-3.2-0/llvm-3.2/projects/compiler-rt/lib/asan/tests/asan_test.cc:2052: Failure Death test: TestGCDSourceCancel() Result: died but not with expected error. Expected: Shadow byte and word Actual msg: [ DEATH ] ================================================================= [ DEATH ] ==88865== ERROR: AddressSanitizer: heap-buffer-overflow on address 0x00010cf70e4a at pc...

Dear all, I am having a (really) hard time getting pyears to work together with a ratetable to give me the number of expected events (deaths). I have the following data: dos, date of surgery, as.Date dof, date of last follow-up, as.Date dos, date of surgery, as.Date sex, gender, as.factor (female,male) ev, event(death), 0= censored at time point dof, 1=death at time point dof Could someone please help. I am searching the web 5 days...

I am interning in a computer science lab and I'm very new to R. The language basics are clear, but this particular problem isn't: I have a very large dataframe called "data" which holds data from Halo matches. I'm trying to analyze a certain window such that data$deaths>20 and data$deaths<=27. When I enter the line kills = data$kills[data$deaths>20] or any single condition, it selects the region I want. However, when I input kills = data$kills[data$deaths>20 && data$deaths<=27] "kills" becomes the entire vector data$kills, wi...

Hi everyone, I am not new to R, but new to running survival models in R. I am trying to create some basic KM curves, using the following code: library(survival) library(KMsurv) (import data etc - basic right censored, with continuously observed time of death) sleepfit <- survfit(Surv(timeb, death), data = sleep) Here timeb is measured is survival in years, death is a 1/0 indicator (1 = event = death), and the dataset is called sleep. But I get the following error Error in UseMethod("survfit", formula) : no applicable method for '...

...one of my current issue. Although my working setting implies a different parametric model (which cannot be framed in the glm), I guess that what I'll get from the following example it would help for the next steps. Anyway here it is. Firstly I simulated from a series of exposures, a series of deaths (given a model, Gompertz, and a probability distribution, Poisson). Then a multiply both deaths and exposures by a constant. Finally I fitted with glm the simulated data sets and I compared the deviance residuals. They're different by a certain constant, but I could not get a meaningful relati...

...147,0):dlm_send_remote_unlock_request:356 ERROR: status = -107 Mar 8 07:20:50 groupwise-1-mht last message repeated 255 times Mar 8 07:20:53 groupwise-1-mht kernel: (4377,0):dlm_send_remote_convert_request:398 ERROR: status = -107 Mar 8 07:20:53 groupwise-1-mht kernel: (4377,0):dlm_wait_for_node_death:371 2062CE05ABA246988E9CCCDAE253F458: waiting 5000ms for notification of death of node 2 Mar 8 07:20:58 groupwise-1-mht kernel: (4377,0):dlm_send_remote_convert_request:398 ERROR: status = -107 Mar 8 07:20:58 groupwise-1-mht kernel: (4377,0):dlm_wait_for_node_death:371 2062CE05ABA246988E9CCCDAE25...

Dear friends. I have these well known data on horsekicks: years <- c(109, 65, 22, 3 , 1, 0) deaths <- 0:5 and get a nice but wrong fit from summary(z1 <- glm(years~deaths,family=poisson)) Can I take away the intercept ? Best wishes Troels Ring, Aalborg -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~...

...e raw form and convert it to Lx (% survival) for a number of treatments. The raw data has the number of days lived for each individual for the respective treatment. Thus, for example, when R selects the data for a single treatment I end up with the following vectors: >day=seq(from=0,to=6) >deaths=c(0,0,2,0,0,1,6) where deaths is the number of deaths on a given day. Now I need to create a new vector with the number alive for each day and this is where I'm stuck... I've tried to work various for and while loops but haven't had success. The vector should be: >Alive=c(9,9,7,7...

...) Error writing 4 bytes to client. -1. (Connection reset by peer) Anyone know why? My configuration has been in http://www.nerdgroup.org/doc/samba+ldap+qmail.txt Thanks -- Fernando Ribeiro - GPG-KEY: 0x8D7255F4 Linux Counter: #273768 - ICQ: 175630330 Linux Professional Institute - LPIC-1 Death the graph! Death the mouse! Death patents! Death closed standards! http://www.nerdgroup.org http://musb.nerdgroup.org -------------------------------------- "Grandes mentes discutem id?ias; Mentes medianas discutem eventos; Mentes pequenas discutem pessoas."

...pdf My data (used in my study) is 'right censored' i.e. my variable corresponding to 'event' indicates whether a person is alive (0) or dead (1) at date last seen and my 'time' indicates time from transplant to date of last contact (where this is time from transplant to death if person has died or time from transplant to date last seen if person is still alive). Now I am using function, rcorr.cens http://lib.stat.cmu.edu/S/Harrell/help/Hmisc/html/rcorr.cens.html This function involves use of Surv. Now here is a section of my syntax: > time<-data$ovsrecod &gt...

On a fresh Ubuntu, sudo apt-get install wine wine iexplore yields a white screen of death with a title bar for Wine. I tried winecfg and tweaking settings, and I can get a desktop with IE WSOD again, so no real progress. The shell outputs a ton of errors/debug info, pasted below... What have I done wrong? wine iexplore fixme:ole:CoResumeClassObjects stub fixme:shdocvw:go_home stub f...

Dear all, I'm struggling with predicting "expected time until death" for a coxph and survreg model. I have two datasets. Dataset 1 includes a certain number of people for which I know a vector of covariates (age, gender, etc.) and their event times (i.e., I know whether they have died and when if death occurred prior to the end of the observation period). Dat...

Hello, I would like to use the spatstat package to generate realizations from a spatial birth death process. The death rate is dependent on the pairwise interaction function. 1 - I'd like to know how to define my own interaction potential function 2 - I'm not sure which model to use for a birth death process (my parameters are the birth and death rates) Any help would be greatly appreciat...

...alysis in industry) uses parametric almost exclusively, since the tail of the distribution is of greatest interest. Your use of predict(,type='response') is almost correct, there is just the math detail that the Weibull fits on a log scale, so the returned value is a geometric mean time to death rather than an arithmetic mean. The suggestion to use ordinary regression on the observed times is wrong. Censored data is more complex than that. Terry Therneau

Dear R users, using dotplot (R1.2.0, WinNT4.0), I am trying to change the character size of the labels of the points: > # example > data(VADeaths) > dotplot(VADeaths, main = "Death Rates in Virginia - 1940") > # I'd like to have smaller character size of the labels (for age and population groups) > ?dotplot > # for argument "cex", this says: "Setting cex to a value smaller than one can be a useful wa...

...per limits of the age-range categories (typically used by paleodemographers), and the third column is the number of skeletons assigned to each group. I've co-opted some R script for fitting a Gompertz-Makeham hazard model to this data... ############################## GM.naga <- function(x,deaths=naga) { a2=x[1] a3=x[2] b3=x[3] shift<-15 nrow<-NROW(deaths) S.t<-function(t) { return(exp(-a2*(t-shift)+a3/b3*(1-exp(b3*(t-shift))))) } d<-S.t(deaths[1:nrow,1])-S.t(deaths[1:nrow,2]) obs<-deaths[,3] lnlk<-as.numeric(crossprod(obs,log(d)))...

Hello all, I am hoping to use survival analysis to examine whether parasite attack increases nest death in a species of social wasp. I therefore have data for 1. Whether the nest "died" in the 6 week census period ("Status", where 1=died, 0=survived) 2. The day number of death/last recorded day it was observed alive. 3. Whether the nest was attacked by the parasite (0/1 as with 1...

...on-negative at all times. What is the right way to do this in R? Thanks, Jeremy P.S., Below is a simplified version of the code I use to try to do this, but I am not sure that it is theoretically right dynmodel <- function(t,y,p) { ## Initialize parameter values birth <- p$mybirth(t) death <- p$mydeath(t) recover <- p$myrecover beta <- p$mybeta vaxeff <- p$myvaxeff vaccinated <- p$myvax(t) vax <- vaxeff*vaccinated/100 ## If the state currently has negative quantities (shouldn't have), then reset to reasonable values for computing meaningful derivatives...