[18734] in s-news-athena
[S] Error in function...
daemon@ATHENA.MIT.EDU (Luca Boni)
Sat Jul 24 11:34:44 1999
Message-Id: <3.0.6.32.19990724172907.009db100@ermes.cba.unige.it>
Date: Sat, 24 Jul 1999 17:29:07 +0200
To: s-news@wubios.wustl.edu
From: Luca Boni <boni@ermes.cba.unige.it>
Mime-Version: 1.0
Content-Type: text/plain; charset="us-ascii"
I wrote this simple function:
"ct.planning"<-
function( tev, tref, mc, hratio, accrual, tmin, nc.c = 0, nc.i = 0 )
{
totp<-uniroot( function(x) e.ev( tref, mc, hratio, x, accrual, tmin, nc.c,
nc.i ) - tev, c(1,50000))$root
totp
}
(Note that e.ev is another function that I wrote and it works well).
When I make a call to ct.planning, an error occurs:
Error in f(1): Object "mc" not found
I include the function e.ev also:
"e.ev"<-
function(tref, mc, hratio, yr.acc, accrual, tmin, nc.c = 0, nc.i = 0)
{
label <- c("Control", "Intervention")
#Find exponential hazard for control group
lamc <- - log(1 - mc)/tref
lamc.c <- lamc + nc.c/100
#Find exponential hazard and mortality for intervention group
lami <- lamc * hratio
lami.c <- lami + nc.i/100
mi <- 1 - exp( - lami * tref)
#Find expected deaths in control group
edc <- ((yr.acc * accrual * lamc)/(2 * lamc.c)) * (1 - exp( - lamc.c *
tmin)/(lamc.c *
accrual) * (1 - exp( - lamc.c * accrual)))
edi <- ((yr.acc * accrual * lami)/(2 * lami.c)) * (1 - exp( - lami.c *
tmin)/(lami.c *
accrual) * (1 - exp( - lami.c * accrual)))
ed <- c(edc, edi)
names(ed) <- label
cat("\nNumber of expected deaths\n")
print(ed)
cat("\nNumber of total expected deaths\n")
c(Total = edc + edi)
}
Any and all suggestions are welcome.
Thanks in advance to all.
Luca Boni
\\\|///
\\ - - //
( @ @ )
----o00o-(_)-o00o----------------------------------------------------
Luca Boni, MD
Unit of Clinical Epidemiology and Trials
Advanced Biotechnology Center-National Institute for Cancer Research
GENOVA-ITALY
Tel.: +39 10 5737477
Portable Tel.: +39 335 6777997
E-mail: boni@ermes.cba.unige.it
Ooooo
----ooooO----( )---------------------------------------------------
( ) ) /
\ ( (_/
\_)
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news
