Given my reproducible example:
new.ex<-structure(list(TEC = c(0.21, 0.077, 0.06, 0.033, 0.014, 0.007,
0.21, 0.077, 0.01, 0.033, 0.05, 0.014), LR = c(FALSE, FALSE,
TRUE, FALSE, TRUE, FALSE, FALSE, FALSE, TRUE, FALSE, TRUE, TRUE
), group = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L), .Label = c("1", "2"), class = "factor")), .Names = c("TEC",
"LR", "group"), row.names = c(NA, -12L), class = "data.frame")
And this few lines of code:
library(NADA)
out<-with(new.ex,
cenfit(TEC, LR, group)
)
out
giving the following output:
n n.cen median mean sd
group=1 6 2 0.033 0.05827778 0.08357853
group=2 6 3 0.033 0.06983333 0.07925407
I would like to add one more result for each group to the above output,
namely “sum”, computed as the product of “n” times “mean”
This is pretty much similar (a slight variation) to a question I posted
earlier in:
http://r.789695.n4.nabble.com/How-to-modify-a-S4-function-in-the-package-NADA-td4649586.html
But in this case I have some problems in modifying the cenfit() function
dealing with group as a factor
My objective is to modify the original function cenfit() so that to get also
the computation of “sum” as product of “n” times “mean”
For some reasons I can not properly understand, I’m not able to successfully
modify my earlier attempt (which was not accounting for groups)
mycenfit <- function(x) {
s = summary(x)
c(n = nrow(s), n.cen = nrow(s) - sum(s$n.event), median = median(x),
mean = mean(x)[["mean"]], sd = sd(x), sum=mean(x)[["mean"]]*length(x))
}
how to change it in order to properly deal with groups?
Thank you for any help
max
--
View this message in context:
http://r.789695.n4.nabble.com/function-customization-tp4649711.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
