On Jun 25, 2010, at 9:31 PM, Julia Choi wrote:
> Dear All:
>
> I have been trying to remove "Inf" values from a series of data
> (multiple ratio operations). I have tried the na.rm() function (and
> tried this inside mean(), median(), and var()) but this does not
> always work. Though the !is.finite() or is.infinite() functions
> appear to be very promising, this would require the creation of an
> intermediary variable and would make the direct creation of a
> data.frame() from a series of ratios much less convenient.
>
> Would anyone have a hint on a function similar to na.rm(), except it
> removes Inf values? Would anyone in the R project be possibly
> interested on adding this functionality?
>
> Minimal example:
>
> TissueRatios = rbind(
> rbind(c(median(XmodA$IDPG[XmodA$Organ=="Tumor" &
XmodA$Time=="48h"]/
> XmodA$IDPG[XmodA$Organ=="Blood" &
XmodA$Time=="48h"], na.rm=TRUE),
> var(XmodA$IDPG[XmodA$Organ=="Tumor" &
XmodA$Time=="48h"]/XmodA
> $IDPG[XmodA$Organ=="Blood" & XmodA$Time=="48h"],
na.rm=TRUE), 1, 1,
> 1),
> c(median(XmodA$IDPG[XmodA$Organ=="Tumor" &
XmodA$Time=="72h"]/XmodA
> $IDPG[XmodA$Organ=="Blood" & XmodA$Time=="72h"],
na.rm=TRUE),
> var(XmodA$IDPG[XmodA$Organ=="Tumor" &
XmodA$Time=="72h"]/XmodA
> $IDPG[XmodA$Organ=="Blood" & XmodA$Time=="72h"],
na.rm=TRUE), 2, 1,
> 1),
> c(median(XmodA$IDPG[XmodA$Organ=="Tumor" &
XmodA$Time=="96h"]/XmodA
> $IDPG[XmodA$Organ=="Blood" & XmodA$Time=="96h"],
na.rm=TRUE),
> var(XmodA$IDPG[XmodA$Organ=="Tumor" &
XmodA$Time=="96h"]/XmodA
> $IDPG[XmodA$Organ=="Blood" & XmodA$Time=="96h"],
na.rm=TRUE), 3, 1,
> 1)), . . .)
>
That is not really a minimal example. For one thing it will not parse
correctly with the three dots at the end. I think it is unreasonable
to expect us to parse out some pattern from that curious mixture of
logical and numeric operators. It is a fairly simple operation to use
functions like is.infinite in subset(), but I cannot tell whether that
is what would suffice since I cannot see the pattern in the reference
to dataset that we cannot see.
>
> Thank you in advance for any help!
> Julia
>
> ______________________________________________
> R-help at 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.