Displaying 1 result from an estimated 1 matches for "repfinalvalue".
2009 May 29
1
final value of nnet with censored=TRUE for survival analysis
...[2,] 1 0 0 0 0 0
[3,] 0 0 1 1 1 1
[4,] 0 1 0 0 0 0
...
where 0,6,17,... are the time values and record 1 has an event in time period 37 and record 3 is censored in time period 20.
I tried to reproduce the final value with this term:
tmp <- rowSums(cat*summary(net)$fitted.values)
repFinalValue <- sum(-log(tmp))
but got the value 21.65438 instead of 21.71108. Is this only a rounding error or what´s wrong with my formula?
My whole test example is printed below.
Thanks,
Andrea.
# Some tests ...
#
###############################################################################
req...