Thank you!
Suhaila
> From: ggrothendi...@gmail.com
> Date: Sun, 13 May 2012 07:58:19 -0400
> Subject: Re: [R] plot.ts error: cannot plot more than 10 series
> To: bell_beaut...@hotmail.com
> CC: r-help@r-project.org
>
> On Sun, May 13, 2012 at 5:02 AM, Suhaila Haji Mohd Hussin
> <bell_beaut...@hotmail.com> wrote:
> >
> > Hello.
> > I'm doing Principal Component Analysis. So I'm trying to plot the new data
> > as time series like this:
> > plot.ts(pr2$x)
> > But it gives me error: Error in plotts(x = x, y = y, plot.type = plot.type,
> > xy.labels = xy.labels, : cannot plot more than 10 series as "multiple"
> > Here's my data dput it's really long so I just put in some essential
> > information only for you to see:
> > structure(list(sdev = c(163.444257295666, 104.206376841326,
> > 97.9593394161277, 85.13008012602, 78.672940030061, 73.7044913582009,
> > 68.9919940515061, 66.4914391167912, 63.5364862493486, 53.0067010985736,
> > 50.3850160673695, 41.0505058921616, 36.6847894205371, 34.4922997539369,
> > 32.8160496561452), rotation = structure(c(-0.341624025780871,
> > -0.465361625382171, -0.254733743609541, -0.501601736275946,
> > -0.32909803586058, 0.10346905261603, 0.0826920791661065, 0.04180674828032,
> > 0.0317514777236264, -0.041673145049323, 0.00552096153263742,
> > 0.302880866160923, 0.0611515661970536, -0.19666860832608,
> > -0.0512002313448844, 0.0532192307778435, -0.279691445123637,
> > -0.0330375822617488, 0.0296862808186525, 0.00754185914289494,
> > -0.0483037204914515, -0.145633004945105, 0.0663131380872606,
> > ...................................0.018728860610997, -0.867607192081349,
> > -0.00746836023138116, -0.197845760644588, 0.354006015245109,
> > -0.0234568336704251, 0.17996661095123, 0.0910352921083386,
> > 0.118413710329641)!
, .D!
> > im = c(20L, 15L), .Dimnames = list( c("er", "pr", "ar", "l1", "l2",
> > "b1", "b2", "h1", "h2", "h3", "p1", "p2", "o1", "o2", "o3", "o4", "o5",
> > "o6", "o7", "o8" ), c("PC1", "PC2", "PC3", "PC4", "PC5", "PC6", "PC7",
> > "PC8", "PC9", "PC10", "PC11", "PC12", "PC13", "PC14", "PC15"))),
> > center = structure(c(94.4986072423398, 94.9637883008357,
> > 71.1838440111421, 179.793871866295, 180.57938718663, 12.6490250696379,
> > 10.5292479108635, 10.8690807799443, 38.2479108635098, 137.543175487465,
> > 0.612813370473538, 45.150417827298, 10.5459610027855, 91.3259052924791,
> > 116.885793871866, 59.0055710306407, 172.598885793872, 14.3676880222841,
> > 5.0974930362117, 3.8857938718663), .Names = c("er", "pr", "ar", "l1",
> > "l2", "b1", "b2", "h1", "h2", "h3", "p1", "p2", "o1", "o2", "o3", "o4",
> > "o5", "o6", "o7", "o8")), scale = FALSE, x =
> > structure(c(-38.7640736739586, -85.3987965131485, 279.953119866405,
> > 121.49850702536, 348.387859356463, -149.756720052989, 1!
14.1!
> > 99587790773, -258.129158201887, 245.434330475019, -74.518625315929
> > , 114.858457403311, 246.989390591038, 264.221327740388,
> > -59.7706080070487, -66.9519268984565, -202.267758377027,
> > -61.2743296889704, -63.3534745591336,
> > .................................................
> > 403375504553, 41.2280423622485, 78.1688899057297, -17.9229807680497,
> > 1.892930692659, -0.121464277679548, -10.9161618391865,
> > 9.77739319210259, -30.4117535843182), .Dim = c(359L, 15L), .Dimnames =
> > list(c("9", "12", "13", "14", "18", "22", "24", "29", "30", "31", "32",
> > "33", "34", "44", "49", "53", "54", "58", "59", "60", "62", "66", "69",
> > "70", "75", "76", "80", "87", "88", "93", "94", "96", "100", "101",
> > "102", "106", "114", "115", "116", "126", "128", "129", "131", "133",
> > "135", "137", "139", "141", "142", "144", "147", "151", "154", "157",
> > "158", "160", "161", "164", "169", "170", "171", "175", "179", "180",
> > "181", "186", "193", "195", "200", "205", "206", "209", "218", "220",
> > "221", "222", "224", "225", "226", "228", "231", "232", "233", "235",
> > "239", "243", "244", "253", "256", "260", "264", "265", "266", "268",
> > "269", "272", "275", "280",
> > ......................................................... "1025", "1032",
> > "1!
033"!
> > , "1035", "1036", "1038", "1039", "1041", "1042", "1046", "1047",
> > "1048", "1051", "1055", "1059", "1060", "1064", "1072", "1076"),
> > c("PC1", "PC2", "PC3", "PC4", "PC5", "PC6", "PC7", "PC8", "PC9",
> > "PC10", "PC11", "PC12", "PC13", "PC14", "PC15")))), .Names = c("sdev",
> > "rotation", "center", "scale", "x"), class = "prcomp")
> > Please help,Suhaila
>
>
> The dput output seems to have gotten messed up in transmission so lets
> use an artificial example:
>
> set.seed(123)
> DF <- as.data.frame(matrix(rnorm(600), nc = 12)) # 12 columns
> p <- prcomp(DF)
>
> # error - cannot plot more than 10
> plot.ts(p$x)
>
> # try it several ways with plot.zoo
>
> library(zoo)
> plot(as.zoo(p$x)) # ok
>
> # or - all on one panel
>
> k <- ncol(p$x)
> plot(as.zoo(p$x), screen = 1, col = 1:k)
>
> # or - 1st two on 1st panel, 2nd two on 2nd panel, etc.
>
> plot(as.zoo(p$x), screen = seq(0, len = k) %/% 2, col = 1:2)
>
> --
> Statistics & Software Consulting
> GKX Group, GKX Associates Inc.
> tel: 1-877-GKX-GROUP
> email: ggrothendieck at gmail.com
[[alternative HTML version deleted]]
______________________________________________
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.
