matrix(age1, ncol=length(age1), nrow=length(age2), byrow=TRUE) + age2 can be replaced by
outer(age2, age1, "+") and the same for 'p'. Best, Dimitris ---- Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://med.kuleuven.be/biostat/ http://www.student.kuleuven.be/~m0390867/dimitris.htm ----- Original Message ----- From: "Rainer M Krug" <[EMAIL PROTECTED]> To: "r-help" <[EMAIL PROTECTED]> Sent: Tuesday, January 15, 2008 1:24 PM Subject: [R] Looking for simpler solution to probabilistic question > Hi > > I have two processes which take with a certain probability (p1 and > p2) x > number of years to complete (age1 and age2). As soon as thge first > process is completed, the second one begins. I want to calculate the > time it takes for the both processes to be completed. > > I have the following script which gives me the answer, butI think > there > must be a more elegant way of doing the calculations between the > ##### > > Any ideas welcom, > > Rainer > > > age1 <- 1:20 > p1 <- dnorm(age1,10,2) > > age2 <- 10:30 > p2 <- dnorm(age2,15,1) > > ###### > age <- matrix(age1, ncol=length(age1), nrow=length(age2), > byrow=TRUE) + age2 > > p <- matrix(p1, ncol=length(p1), nrow=length(p2), byrow=TRUE) + p2 > p <- as.numeric(tapply(p, age, sum)) > p <- p / sum(p) > ##### > > age <- sort(unique(as.numeric(age))) > > plot(age, p) > > lines(age1, p1) > lines(age2, p2) > > ______________________________________________ > R-help@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. > Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm ______________________________________________ R-help@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.