library(utils)
# many series at once
data(gp.rwl)
gp.pt <- powt(gp.rwl)
hist(summary(gp.rwl)$skew)
hist(summary(gp.pt)$skew)
# single series
gp01A <- gp.rwl[, "01A"]
names(gp01A) <- rownames(gp.rwl)
gp01A.pt <- powt.series(gp01A,rescale=TRUE)
plot(gp01A.pt,gp01A)
abline(c(0,1))
Run the code above in your browser using DataLab