On Jan 17, 2010, at 8:17 PM, maiya wrote:
There must be a very basic thing I am not getting...
I'm working with some entropy functions and the convention is to use
log(0)=0.
I suppose the outcome of that effort may depend on whether you have
assumed the needed godlike capacities to change the laws of
mathematics. But I suppose that as the Earth mother that might occur
to you. Go ahead, define a new mathematics.
--
David
So I wrote a function:
llog<-function(x){
if (x ==0) 0 else log(x)
}
which seems to work fine for individual numbers e.g.
llog(0/2)
[1] 0
but if I try whole vectors or tables:
p<-c(4,3,1,0)
q<-c(2,2,2,2)
llog(p/q)
I get this:
[1] 0.6931472 0.4054651 -0.6931472 -Inf
What am I missing?
Thanks!
Maja
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
______________________________________________
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.