R notacji naukowej działkami

Question

mam prosty wykres:

#!/usr/bin/Rscript                      

png('plot.png') 

y <- c(102, 258, 2314)                   
x <- c(482563, 922167, 4462665) 

plot(x,y) 
dev.off() 

R wykorzystuje 500, 1000, 1500, itp na osi y. Czy istnieje sposób, w jaki mogę użyć notacji naukowej dla osi y i umieścić * 10^3 na górze osi, jak na poniższym rysunku?

Answer 1

Jest to swego rodzaju hacky sposób, ale nie ma nic złego.

plot(x,y/1e3, ylab="y /10^3") 

Answer 2

Jak masz etykiety na swoją oś zależy od używanego systemu wykreślania (podstawa, ggplot2 lub kratownica) można użyć funkcji z pakietu scales formatowanie liczb osi:

library(scales) 
x <- 10^(1:10) 
scientific_format(1)(x) 
[1] "1e+01" "1e+02" "1e+03" "1e+04" "1e+05" "1e+06" "1e+07" "1e+08" "1e+09" "1e+10" 

tu przykład stosując ggplot2:

library(ggplot2) 
dat <- data.frame(x = c(102, 258, 2314),                  
        y = c(482563, 922167, 4462665)) 

qplot(data=dat,x=x,y=y) + 
    scale_y_continuous(label=scientific_format(digits=1))+ 
    theme(axis.text.y =element_text(size=50)) 

EDIT The OP ma specyficzne potrzeby. Oto kilka pomysłów, których użyłem tutaj, aby to osiągnąć:

1. Można dostosować etykiety wydruku za pomocą funkcji axis.
2. Zastosowanie mtext umieścić tekst w zewnętrznym obszarze działki
3. używając wyrażenia zysków dzięki funkcji plotmath ...

y <- c(102, 258, 2314)                   
x <- c(482563, 922167, 4462665) 
plot(x,y,ylab='',yaxt='n') 
mtext(expression(10^3),adj=0,padj=-1,outer=FALSE) 
axis(side=2,at=y,labels=round(y/1000,2)) 

Answer 3

Podobna technika jest użycie eaxis (rozszerzony/oś inżynieryjna) z pakietu sfsmisc.

To działa tak:

library(sfsmisc) 

x <- c(482563, 922167, 4462665) 
y <- c(102, 258, 2314) 

plot(x, y, xaxt="n", yaxt="n") 

eaxis(1) # x-axis 
eaxis(2) # y-axis