R在法线轴上绘制logit转换数据（逻辑转换）

Question

我想以易于理解的方式绘制 logit 转换数据（最初是百分比），即在带有逻辑转换 y 轴的图中。

library(tidyverse)

set.seed(123)

# Logit transformation function
logit <- function(x){
  x <- x/100
  log(x/(1-x))
}

# Create sample data
df <- data.frame(x = seq(1,10,1), 
           y = runif(10, 0, 10)) %>% 
  mutate(y = logit(y))

# Plot
ggplot(df, aes(x,y)) +
  geom_point()

我注意到scale_y_continuous 函数有一个trans 参数，它应该能够完成这项工作。但是，logit 函数是内置的，但与 logistic 相对的函数不是。

如何正常绘制 logit 转换数据（即，使用逻辑转换 y 轴）？

Answer 1

如果数据是百分比，并且您希望 y 轴进行 logit 转换，并且 y 轴标签保持为百分比，您可以从 scales 创建一个新的 trans 类，在下面的示例中，我使用与logit 函数等效的 qlogis：

library(scales)
logit_perc <- trans_new("logit perc",
transform = function(x)qlogis(x/100),
inverse = function(x)100*plogis(x)
)

set.seed(1001)
df <- data.frame(x = seq(1,10,1),y = runif(10, 0, 10)) 

ggplot(df,aes(x,y)) + geom_point() + coord_trans(y = logit_perc)

如果您已经有 logit 中的值，并且您想使用 logit 比例绘制百分比标签：

df = df %>% mutate(y=logit(y))

# Plot
ggplot(df, aes(x,y)) + geom_point() + 
scale_y_continuous(labels=function(x)round(100*plogis(x),digits=1))

而原来的情节是这样的：

ggplot(df, aes(x,y)) + geom_point()