search for: 103l

...2L, 2L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("Fri", "Mon", "Thu", "Tue", "Wed"), class = "factor"), SITE = c(101L, 101L, 102L, 102L, 103L, 103L, 104L, 104L, 101L, 101L, 102L, 102L, 103L, 103L, 104L, 104L, 101L, 101L, 102L, 102L, 103L, 103L, 104L, 104L, 101L, 101L, 102L, 102L, 103L, 103L, 104L, 104L, 101L, 101L, 102L, 102L, 103L, 103L, 104L, 104L), TOD = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1...

...4 [3,] 4 7 2 [4,] 0 3 8 a1 <- structure(c(7, 4, 4, 0, 6, 2, 7, 3, 8, 4, 2, 8), .Dim = c(4L, 3L)) > a2 [,1] [,2] [,3] [1,] 101 102 103 [2,] 101 102 103 [3,] 101 102 103 [4,] 101 102 103 a2 <- structure(c(101L, 101L, 101L, 101L, 102L, 102L, 102L, 102L, 103L, 103L, 103L, 103L), .Dim = c(4L, 3L)) ------ I want to get a3: > a3 [,1] [,2] [,3] [1,] 102 101 103 [2,] 102 101 103 [3,] 103 101 102 [4,] 101 102 103 where the rows of a3 are the rows of a2 sorted according to the rows of a1. I can get the necessary sorting index: > app...

...and seq is the sequence variable that I've created. Thanks very much for the help. david freedman, atlanta ds=data.frame(list(id = c(1L, 1L, 1L, 1L, 8L, 8L, 16L, 16L, 16L, 16L, 16L, 19L, 32L, 64L, 64L, 64L, 64L, 64L, 64L, 64L, 64L, 64L, 79L, 79L, 80L, 80L, 80L, 80L, 85L, 86L, 96L, 96L, 96L, 103L, 103L, 106L, 106L, 106L, 106L, 106L, 106L, 106L, 140L, 140L, 144L, 144L, 144L, 144L, 144L, 144L, 144L, 146L, 146L, 146L, 146L, 160L, 160L, 160L, 160L, 160L, 160L, 164L, 164L, 176L, 176L, 176L, 176L, 176L, 176L, 176L, 176L, 181L, 190L, 192L, 192L, 192L, 192L, 192L, 192L, 197L, 197L, 197L, 224L, 224L...

..."B", "C", "D", "E", "F", "G"), class = "factor"), total = c(135L, 118L, 121L, 64L, 53L, 49L, 178L, 123L, 128L, 127L, 62L, 129L, 126L, 99L, 183L, 45L, 57L, 45L, 72L, 30L, 71L, 123L, 89L, 102L, 60L, 44L, 59L, 124L, 145L, 126L, 103L, 67L, 97L, 66L, 76L, 108L, 36L, 48L, 41L, 69L, 47L, 57L, 167L, 136L, 176L, 85L, 36L, 82L, 222L, 149L, 171L, 145L, 122L, 192L, 136L, 164L, 154L, 46L, 57L, 57L, 70L, 55L, 102L, 111L, 152L, 204L, 41L, 46L, 103L, 156L, 148L, 155L, 103L, 124L, 176L, 111L, 142L, 187L, 43L, 52L, 75L, 64L, 91L, 78L, 19...

Hi All, I have a scatter plot produced using ggplot2 and I want to add the regression line to this scatter plot. I suppose I can use geom_smooth() to do this, but for the sake of learning ( I am new both to R and ggplot2), I want to try and add it as a function (something that curve() does in the standard R plotting). I did some search and found that stat_function() can be used for this. But

..., 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 111L, 112L, 113L, 114L, 115L, 116L, ...... ...... ..... V2 = structure(c(278L, 23L, 11L, 169L, 14L, 86L, 94L, 51L, 37L, 43L, 22L, 169L, 49L, 120L, 18L, 60L, 42L, 41L, 38L, 64L, 38L, 32L, 140L, 146L, 106L, 26L, 46L, 65L, 17L, 106L, 20L, 33L,...

..., 15L, 25L, 39L, 1L, 71L, 95L, 19L, 56L, 62L, 76L, 85L, 99L, 111L, 42L, 4L, 52L, 61L, 67L, 91L, 13L, 24L, 36L, 38L, 7L, 81L, 82L, 57L, 64L, 100L, 113L, 21L, 33L, 44L, 78L, 54L, 73L, 83L, 108L, 16L, 28L, 2L, 58L, 65L, 86L, 101L, 32L, 43L, 77L, 93L, 104L, 10L, 48L, 53L, 97L, 110L, 18L, 27L, 75L, 92L, 103L, 31L, 9L, 47L, 70L, 72L, 106L, 14L, 26L, 40L, 79L, 98L, 112L, 68L, 88L, 105L, 20L, 11L, 49L, 63L, 5L, 80L, 50L, 59L, 89L, 114L, 22L, 34L, 45L), .Label = c("1/14/08 12:26", "1/15/08 15:25", "1/15/08 16:32", "1/16/08 11:36", "1/16/08 15:29", "1/2...

...le below, the behavior is the opposite of the help page, if I am understanding it correctly, and a warning is thrown when it does seem to work as desired. c1 <- structure(list(Contr = c(0.028, 0.043, 0.064, 0.097, 0.146, 0.219 ), Correct = c(34L, 57L, 94L, 152L, 160L, 160L), Incorrect = c(126L, 103L, 66L, 8L, 0L, 0L)), .Names = c("Contr", "Correct", "Incorrect" ), row.names = c("13", "15", "17", "19", "21", "23"), class = "data.frame") q25 <- rep( qlogis( 0.25 ), nrow(c1) ) # offset defined...

...ot;, "y"), row.names = c(1L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 53L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L), class = "data.frame") library(nlme) grpdat = groupedData(y ~ age | ptno, data = tmpdat, FUN = mean, labels = list(x="Age", y="Y")) fit.lme = lme(grpdat) p...

...you go >dput(eliza) structure(c(428L, 359L, 353L, 498L, 169L, 46L, 29L, 444L, 194L,? 439L, 284L, 79L, 226L, 105L, 106L, 326L, 234L, 196L, 417L, 190L,? 480L, 275L, 478L, 341L, 197L, 395L, 85L, 295L, 70L, 18L, 471L,? 218L, 19L, 43L, 280L, 427L, 229L, 5L, 205L, 494L, 393L, 307L,? 465L, 383L, 265L, 103L, 348L, 166L, 58L, 320L, 14L, 107L, 440L,? 24L, 347L, 419L, 135L, 202L, 294L, 173L, 121L, 358L, 51L, 231L,? 118L, 109L, 182L, 219L, 42L, 73L, 435L, 389L, 431L, 81L, 149L,? 388L, 415L, 256L, 312L, 446L, 45L, 162L, 485L, 154L, 55L, 153L,? 137L, 290L, 93L, 40L, 110L, 167L, 315L, 266L, 317L, 314L, 264L,...

..., -69L, -57L, -59L, -67L, -73L, -81L, -90L, -87L, -82L, -80L, -80L, -81L, -83L, -90L, -89L, -107L, -122L, -128L, -113L, -124L, -93L, -89L, -92L, -85L, -79L, -84L, -83L, -148L, -142L, -127L, -130L, -123L, -131L, -85L, -87L, -90L, -81L, -65L, -81L, -119L, -121L, -94L, -90L, -117L, -95L, -69L, -122L, -103L, -85L, -106L)), .Names = c("mine", "dCu65", "dD"), class = "data.frame", row.names = c(NA, -130L))

I have been working with the pls procedure and have problems getting the procedure to work with matrix or frame data. I suspect the problem lies in my understanding of frames, but can't find anything in the documentation that will help. Here is what I have done: I read in an 10000 x 8 table of data, and assign the first four columns to matrix A and the second four to matrix B pls <-

...umeric vector of length either one or equal to the number of cases. ... but with the following example, I get an error if the offset is of length 1 c1 <- structure(list(Contr = c(0.028, 0.043, 0.064, 0.097, 0.146, 0.219 + ), Correct = c(34L, 57L, 94L, 152L, 160L, 160L), Incorrect = c(126L, + 103L, 66L, 8L, 0L, 0L)), .Names = c("Contr", "Correct", "Incorrect" + ), row.names = c(NA, 6L), class = "data.frame") glm(cbind(Correct, Incorrect) ~ Contr - 1, binomial, data = c1, offset = qlogis(0.25)) Error in model.frame.default(formula = cbind(Correct, In...

...= c(1L, 2L, 4L, 5L, 6L, 7L, 8L, 12L, 14L, 16L, 17L, 20L, 21L, 23L, 24L, 26L, 27L, 28L, 30L, 31L, 35L, 36L, 37L, 38L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 55L, 56L, 58L, 59L, 61L, 62L, 65L, 66L, 67L, 69L, 71L, 72L, 76L, 78L, 79L, 80L, 88L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 112L), class = "data.frame") The data (in the "structure") is removed, but the labels for all LithClass remain. If I then produce a boxplot or some other graphic where LithClass is involved, blank plots are produced for the deleted LithClass...

Folks, I have a data frame with columns 200401, 200402, ..., 201207, 201208. These represent years/months. What would be the best way to sum these columns by year? What about by quarter? Thanks for your time, KW -- [[alternative HTML version deleted]]

...5L, 46L, 47L, 49L, 50L, 51L, 52L, 53L, 54L, 55L, 56L, 57L, 58L, 60L, 61L, 68L, 79L, 90L, 93L, 94L, 95L, 96L, 97L, 98L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 88L, 89L, 91L, 92L, 99L, 110L, 121L, 123L, 124L, 125L, 126L, 127L, 128L, 100L, 101L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 111L, 112L, 113L, 114L, 115L, 116L, 117L, 118L, 119L, 120L, 122L, 2L, 13L, 24L, 27L, 28L, 29L, 30L, 31L, 32L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 25L, 26L, 33L, 34L, 35L, 36L), .Label = c("", &qu...

...155L, 153L, 151L, 152L, 148L, 150L, 149L, 147L, 145L, 146L, 142L, 143L, 144L, 139L, 140L, 141L, 138L, 137L, 136L, 135L, 133L, 134L, 132L, 130L, 131L, 127L, 129L, 128L, 126L, 124L, 125L, 121L, 122L, 123L, 120L, 118L, 119L, 117L, 115L, 116L, 112L, 114L, 113L, 109L, 110L, 111L, 106L, 108L, 107L, 105L, 103L, 104L, 102L, 100L, 101L, 99L, 97L, 98L, 94L, 96L, 95L, 93L, 91L, 92L, 90L, 88L, 89L, 85L, 86L, 87L, 84L, 82L, 83L, 81L, 79L, 80L, 76L, 77L, 78L, 75L, 73L, 74L, 72L, 70L, 71L, 69L, 67L, 68L, 66L, 64L, 65L, 61L, 62L, 63L, 60L, 58L, 59L, 55L, 56L, 57L, 54L, 52L, 53L, 51L, 49L, 50L, 46L, 48L, 47L, 43L,...

...ot;color", "ordered", "pchloc"), row.names = c(98L, 99L, 71L, 100L, 70L, 51L, 68L, 83L, 34L, 35L, 82L, 50L, 42L, 46L, 63L, 62L, 43L, 6L, 44L, 87L, 5L, 31L, 61L, 79L, 55L, 36L, 67L, 2L, 32L, 108L, 47L, 66L, 86L, 88L, 122L, 11L, 16L, 94L, 95L, 115L, 123L, 24L, 80L, 78L, 111L, 103L, 3L, 22L, 23L, 10L, 107L, 97L, 75L, 96L, 106L, 102L, 118L, 7L, 110L, 30L, 59L, 76L, 114L, 74L, 104L, 1L, 64L, 73L, 54L, 58L, 19L, 69L, 92L, 14L, 119L, 18L, 127L, 41L, 15L, 8L, 49L, 91L, 65L, 120L, 113L, 93L, 17L, 126L, 29L, 105L, 90L, 77L, 81L, 45L, 38L, 85L, 39L, 53L, 9L, 121L, 125L, 109L, 101L, 4...

...SILIKATG = c(0.0152, 0.0161, 0.0141, 0.0131, 0.0108, 0, 0, 0, 0, 1e-04, 0.0056, 0.0078, 0.0072, 0.0047, 0.0091, 3e-04, 5e-04, 0, 0.0013, 0.0012)), .Names = c("DIET", "SILIKATG"), row.names = c(1L, 2L, 3L, 4L, 5L, 31L, 32L, 33L, 34L, 35L, 61L, 62L, 63L, 64L, 65L, 101L, 102L, 103L, 104L, 105L), class = "data.frame") 1. The first plot I want to do is: plot(test$SILIKATG, ylim=c(0,0.02), type="l") It looks the way I want to, but my question is:* How to change the labels on the x-axis? I would like to have the groups ("G", "GG", &quot...

...r"), V4 = c(3L, 1L, 3L, 2L, 0L, 2L, 1L, 3L, 2L, 1L, 4L, 0L, 1L, 2L, 4L, 2L, 3L, 4L, 1L, 2L ), V5 = c(0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 1L), V6 = c(1831L, 536L, 760L, 472L, 752L, 503L, 576L, 839L, 767L, 257L, 80L, 223L, 120L, 175L, 96L, 200L, 192L, 103L, 119L, 191L)), .Names = c("V1", "V2", "V3", "V4", "V5", "V6"), row.names = c(NA, 20L), class = "data.frame"), ? ? structure(list(V1 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, ? ? 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), V2 =...