[R] adding results to plot
Rui Barradas
ru|pb@rr@d@@ @end|ng |rom @@po@pt
Fri Oct 8 09:15:47 CEST 2021
Hello,
Thanks for the compliment.
The R Core team, to which we must be very grateful for their great work
along so many years, is, for good and obvious reasons, known to be
change resistant and a new method would overload them even more with
maintenance worries so I guess text.htest won't make it to core R.
(And that would be to open a precedent. Is, for instance, text.lm next?)
If it does make it to base R, which I doubt, then the base R package
should be package graphics, right?
In the mean time I have found a small bug, near the end of print.htest
there's
print(x$estimate, digits = digits, ...)
corresponding to my
paste("sample estimates:", round(ht$estimate, digits = digits), sep = "\n")
The bug is that x$estimate/ht$estimate is a named vector and with paste
the names attribute is lost. It doesn't plot "mean of x". I will
probably try to sort this out but make no promises.
Thanks once again,
Rui Barradas
Às 14:47 de 07/10/21, PIKAL Petr escreveu:
> Hallo Rui.
>
> I finally tested your function and it seems to me that it should propagate
> to the core R or at least to the stats package.
>
> Although it is a bit overkill for my purpose, its use is straightforward and
> simple. I checked it for several *test functions and did not find any
> problem.
>
> Thanks and best regards.
>
> Petr
>
>> -----Original Message-----
>> From: Rui Barradas <ruipbarradas using sapo.pt>
>> Sent: Friday, September 17, 2021 9:56 PM
>> To: PIKAL Petr <petr.pikal using precheza.cz>; r-help <r-help using r-project.org>
>> Subject: Re: [R] adding results to plot
>>
>> Hello,
>>
>> *.test functions in base R return a list of class "htest", with its own
>> print method.
>> The method text.htest for objects of class "htest" below is a hack. I
>> adapted the formating part of the code of print.htest to plot text().
>> I find it maybe too complicated but it seems to work.
>>
>> Warning: Not debugged at all.
>>
>>
>>
>> text.htest <- function (ht, x, y = NULL, digits = getOption("digits"),
>> prefix = "", adj = NULL, ...) {
>> out <- list()
>> i_out <- 1L
>> out[[i_out]] <- paste(strwrap(ht$method, prefix = prefix), sep = "\n")
>> i_out <- i_out + 1L
>> out[[i_out]] <- paste0("data: ", ht$data.name)
>>
>> stat_line <- NULL
>> i_stat_line <- 0L
>> if (!is.null(ht$statistic)) {
>> i_stat_line <- i_stat_line + 1L
>> stat_line[[i_stat_line]] <- paste(names(ht$statistic), "=",
>> format(ht$statistic, digits =
>> max(1L, digits - 2L)))
>> }
>> if (!is.null(ht$parameter)) {
>> i_stat_line <- i_stat_line + 1L
>> stat_line[[i_stat_line]] <- paste(names(ht$parameter), "=",
>> format(ht$parameter, digits =
>> max(1L, digits - 2L)))
>> }
>> if (!is.null(ht$p.value)) {
>> fp <- format.pval(ht$p.value, digits = max(1L, digits - 3L))
>> i_stat_line <- i_stat_line + 1L
>> stat_line[[i_stat_line]] <- paste("p-value",
>> if (startsWith(fp, "<")) fp else
>> paste("=", fp))
>> }
>> if(!is.null(stat_line)){
>> i_out <- i_out + 1L
>> #out[[i_out]] <- strwrap(paste(stat_line, collapse = ", "))
>> out[[i_out]] <- paste(stat_line, collapse = ", ")
>> }
>> if (!is.null(ht$alternative)) {
>> alt <- NULL
>> i_alt <- 1L
>> alt[[i_alt]] <- "alternative hypothesis: "
>> if (!is.null(ht$null.value)) {
>> if (length(ht$null.value) == 1L) {
>> alt.char <- switch(ht$alternative, two.sided = "not equal to",
>> less = "less than", greater = "greater than")
>> i_alt <- i_alt + 1L
>> alt[[i_alt]] <- paste0("true ", names(ht$null.value), " is ",
>> alt.char,
>> " ", ht$null.value)
>> }
>> else {
>> i_alt <- i_alt + 1L
>> alt[[i_alt]] <- paste0(ht$alternative, "\nnull values:\n")
>> }
>> }
>> else {
>> i_alt <- i_alt + 1L
>> alt[[i_alt]] <- ht$alternative
>> }
>> i_out <- i_out + 1L
>> out[[i_out]] <- paste(alt, collapse = " ")
>> }
>> if (!is.null(ht$conf.int)) {
>> i_out <- i_out + 1L
>> out[[i_out]] <- paste0(format(100 * attr(ht$conf.int, "conf.level")),
>> " percent confidence interval:\n", " ",
>> paste(format(ht$conf.int[1:2], digits =
>> digits), collapse = " "))
>> }
>> if (!is.null(ht$estimate)) {
>> i_out <- i_out + 1L
>> out[[i_out]] <- paste("sample estimates:", round(ht$estimate,
>> digits = digits), sep = "\n")
>> }
>> i_out <- i_out + 1L
>> out[[i_out]] <- "\n"
>> names(out)[i_out] <- "sep"
>> out <- do.call(paste, out)
>> if(is.null(adj)) adj <- 0L
>> text(x, y, labels = out, adj = adj, ...)
>> invisible(out)
>> }
>>
>>
>> res <- shapiro.test(rnorm(100))
>> plot(1,1, ylim = c(0, length(res) + 1L))
>> text(res, 0.6, length(res) - 1)
>> res
>>
>> res2 <- t.test(rnorm(100))
>> plot(1,1, ylim = c(0, length(res2) + 1L))
>> text(res2, 0.6, length(res2) - 1L)
>> res2
>>
>>
>> Hope this helps,
>>
>> Rui Barradas
>>
>>
>>
>> Às 15:12 de 16/09/21, PIKAL Petr escreveu:
>>> Dear all
>>>
>>> I know I have seen the answer somewhere but I am not able to find it.
>> Please
>>> help
>>>
>>>> plot(1,1)
>>>> res <- shapiro.test(rnorm(100))
>>>> res
>>>
>>> Shapiro-Wilk normality test
>>>
>>> data: rnorm(100)
>>> W = 0.98861, p-value = 0.5544
>>>
>>> I would like to add whole res object to the plot.
>>>
>>> I can do it one by one
>>>> text(locator(1), res$method)
>>>> text(locator(1), as.character(res$p.value))
>>> ...
>>> But it is quite inconvenient
>>>
>>> I could find some way in ggplot world but not in plain plot world.
>>>
>>> Best regards
>>> Petr
>>>
>>>
>>> ______________________________________________
>>> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> 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.
>>>
More information about the R-help
mailing list