search for: 105l

...c(50, 0, 55, 12, 2, 0, 37, NA, 17, 3, 31 ), min = c(1L, 10L, 11L, 15L, 16L, 18L, 18L, NA, 20L, 22L, 22L ), hour = c(12L, 12L, 12L, 12L, 12L, 12L, 12L, NA, 12L, 12L, 12L), mday = c(13L, 13L, 13L, 13L, 13L, 13L, 13L, NA, 13L, 13L, 13L), mon = c(5L, 5L, 5L, 5L, 5L, 5L, 5L, NA, 5L, 5L, 5L), year = c(105L, 105L, 105L, 105L, 105L, 105L, 105L, NA, 105L, 105L, 105L), wday = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, NA, 1L, 1L, 1L), yday = c(163L, 163L, 163L, 163L, 163L, 163L, 163L, NA, 163L, 163L, 163L), isdst = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, -1L, 1L, 1L, 1L)), .Names = c("sec", "min", &quo...

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

...380), class = c("POSIXct", "POSIXt"), tzone = ""), Score = c(80L, 11L, 81L, 152L, 130L, 122L, 142L, 20L, 1L, 31L, 93L, 136L, 128L, 112L, 48L, 57L, 92L, 108L, 100L, 107L, 81L, 37L, 47L, 70L, 114L, 125L, 99L, 46L, 108L, 106L, 111L, 75L, 75L, 136L, 36L, 13L, 35L, 71L, 105L, 113L, 116L, 116L, 94L, 130L, 102L, 19L, 1L, 33L, 78L, 89L, 115L, 130L, 117L, 105L, 0L, 36L, 80L, 130L, 124L, 124L, 66L, 75L, 57L, 94L, 124L, 165L, 170L, 158L, 127L, 116L, 107L, 118L, 150L, 96L, 101L, 59L, 38L, 46L, 77L, 83L, 54L, 61L, 60L, 55L, 110L, 127L, 129L, 87L, 61L, 70L, 85L, 110L, 131L,...

...(c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L), .Label = c("Well", "Mild", "Moderate", "Impaired"), class = c("ordered", "factor")), Freq = c(64L, 94L, 58L, 46L, 57L, 94L, 54L, 40L, 57L, 105L, 65L, 60L, 72L, 141L, 77L, 94L, 36L, 97L, 54L, 78L, 21L, 71L, 54L, 71L)), .Names = c("ses", "mental", "Freq"), row.names = c(NA, -24L), class = "data.frame") -- Michael Friendly Email: friendly AT yorku DOT ca Professor, Psychology Dept. York Universi...

...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, 68L, 62L, 11...

...740), class = c("POSIXct", "POSIXt"), tzone = ""), Score = c(80L, 11L, 81L, 152L, 130L, 122L, 142L, 20L, 1L, 31L, 93L, 136L, 128L, 112L, 48L, 57L, 92L, 108L, 100L, 107L, 81L, 37L, 47L, 70L, 114L, 125L, 99L, 46L, 108L, 106L, 111L, 75L, 75L, 136L, 36L, 13L, 35L, 71L, 105L, 113L, 116L, 116L, 94L, 130L, 102L, 19L, 1L, 33L, 78L, 89L)), .Names = c("Date", "Score"), row.names = c(NA, 50L), class = "data.frame") I'm trying to round the time down to the nearest half hour, without losing the date. For example, 01/01/2009 1:29 would roun...

...2L, 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/21/08 11:21", "1/21/08 14:18", "1/22/07 13:59", &quot...

..."), 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) plot(augPred(...

...botto at hotmail.com> To: "smartpink111 at yahoo.com" <smartpink111 at yahoo.com> Sent: Wednesday, February 27, 2013 8:47 AM Subject: RE: matrix multiplication there 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...

Hello, I have an input file: http://r.789695.n4.nabble.com/file/n3700031/testOut.txt testOut.txt where col 1 is chromosome, column2 is start of region, column 3 is end of region, column 4 and 5 is base position, column 6 is total reads, column 7 is methylation data, and column 8 is the strand. I would like a summary output file such as:

...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 labels, whi...

..., 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("", "10/1/2008...

...154L, 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...

..., 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, 48L, 52L, 84L, 21L, 60L, 72L, 56L, 124L, 57L, 12L, 112L, 13L, 4L, 26L, 89L, 116L, 27L, 128L, 20L, 117L, 40L, 33L, 28L, 37L, 25L), class = "data.frame") I use this codes (I made them more simple for this post) # per grou...

...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", "L", &q...

I have data in a file named hands.dat, which is given at the end of this question. (It's from a stats textbook example on anova). I'd like to do an aov on this, which I guess would be d <- read.table("~/hands.dat", header=TRUE) aov(Bacterial.Counts ~ Water + Soap + Antibacterial.Soap + Alcohol.Spray, data=d) but this fails. Do I need to break d$Method up into columns for

..., 111L, 142L, 187L, 43L, 52L, 75L, 64L, 91L, 78L, 196L, 314L, 265L, 44L, 39L, 98L, 197L, 273L, 274L, 89L, 91L, 74L, 91L, 112L, 98L, 140L, 90L, 121L, 120L, 161L, 83L, 230L, 266L, 282L, 35L, 53L, 57L, 315L, 332L, 202L, 90L, 79L, 89L, 67L, 116L, 109L, 44L, 68L, 75L, 29L, 52L, 52L, 253L, 203L, 87L, 105L, 234L, 152L, 247L, 243L, 144L, 167L, 165L, 95L, 300L, 128L, 125L, 84L, 183L, 88L, 153L, 185L, 175L, 226L, 216L, 118L, 118L, 94L, 224L, 259L, 176L, 175L, 147L, 197L, 141L, 176L, 187L, 87L, 92L, 148L, 86L, 139L, 122L), country = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,...

Hello, I'm a new R user. I have two different dummy tables with the variable name tb1 and tb2. tb1< v1 v2 v3 v4 "chr1" 22 23 3 "chr1" 36 37 1 "chr1" 54 55 0 "chr1" 77 78 1 "chr2" 80 81 4 "chr2" 85 86 0 "chr2" 99 100 1

Hi there, I have a huge data set with multiple firms years and other firm characteristics. I want to run a regression on the dependent variable and other explanatory variables and calculate the residual terms by grouping the firms in same year and same industry. What I want to do is to divide my obseravtion into sub sample that contains the observation with same fiscal year(FYEAR=1990) and same

...e", "pf_mcl", "year"), row.names = c(NA, 6L), class = "data.frame") myDF2 <- structure(list(studyno = c("J740/4", "J1000/9", "J895/7", "J931/6", "J609/1", "J941/3"), pf_mcl = c(101L, 102L, 103L, 104L, 105L, 106L)), .Names = c("studyno", "pf_mcl"), row.names = c(NA, 6L), class = "data.frame") m <- merge(myDF1, myDF2, by="studyno", all.x=TRUE, all.y=FALSE, suffixes=c(".raw", ".curated")) The results are: > myDF1 studyno date...