search for: day's

Hi, I have a vector in a data frame that looks something like this: day<-c('Day -1','Day 6','Day 10') This vector specifies the order in which several panel will appear in a lattice plots. But the order in which such plots will appear will be the following: Day -1, Day 10, Day 6. Which makes sense, but I cannot name the Days like this: D...

##### BEGIN CODE ###### #!/usr/bin/perl ###### # # --start, -s = The date you would like to start generating regressors #--end, -e = When to stop generating holiday regressros # --scope, -c = D, W for Daily or Weekly respectively (e.g. Does this week have a particular holiday) # --file, -f = Ummm where to write the output silly! # # **NOTE** The EOM holiday is "End of Month" for computer systems this may be important for # extra processing and what n...

Here is the correct version. The old version is the redirect only version of the script. ### BEGIN SCRIPT #### #!/usr/bin/perl ###### # --start, -s = The date you would like to start generating regressors #--end, -e = When to stop generating holiday regressros # --scope, -c = D, W for Daily or Weekly respectively (e.g. Does this week have a particular holiday) # --file, -f = Ummm where to write the output silly! # # **NOTE** The EOM holiday is "End of Month" for computer systems this may be important for # extra processing and what n...

...wall, but?I don't seem to be able to get this right. ?Would somebody?be able to?give me a hint ? This is my code so far, the particle starts at some point and moves towards the wall, but I don't get the "bouncing off" part right?. ? Any help would be much appreciated. ? Juliane ? days=10 circularspace=data.frame(day=c(0:days),xcoord=1,ycoord=1,xvelocity=1,yvelocity=1,xdistwall=0, ydistwall=0, wallxvel=0, wallyvel=0,stochasticxvel=0,stochasticyvel=0) xmax=10 xmin=-10 ymax=10 ymin=-10 mindist=8 plot(xmin:xmax, ymin:ymax, type = "n") circularspace radius=10 timesteplengt...

...the wall, but I don't seem to be able to get this right. Would somebody kindly give me a hint ? This is my code so far, the particle starts at some point and moves towards the wall, but I don't get the "bouncing off" part right . Any help would be much appreciated. Juliane days=10 circularspace=data.frame(day=c(0:days),xcoord=1,ycoord=1,xvelocity=1,yvelocity=1,xdistwall=0, ydistwall=0, wallxvel=0, wallyvel=0,stochasticxvel=0,stochasticyvel=0) xmax=10 xmin=-10 ymax=10 ymin=-10 mindist=8 plot(xmin:xmax, ymin:ymax, type = "n") circularspace radius=10 timesteplengt...

Hello; I am having a problems with the interpretation of models using ordered or unordered predictors. I am running models in lmer but I will try to give a simplified example data set using lm. Both in the example and in my real data set I use a predictor variable referring to 3 consecutive days of an experiment. It is a factor, and I thought it would be more correct to consider it ordered. Below is my example code with my comments/ideas along it. Can someone help me to understand what is happening? Thanks a lot in advance; Catarina Miranda y<-c(72,25,24,2,18,38,62,30,78,34,67,21,9...

CRAN reached 12,000 packages [1] today (December 15, 2017). A huge thank you to the CRAN team! Milestones: 2017-12-15 12000 pkgs (+6.1/day over 165 days) 6910 mnts (+3.2/day) 2017-07-04 11000 pkgs (+6.3/day over 159 days) 6377 mnts (+3.3/day) 2017-01-27 10000 pkgs (+6.3/day over 158 days) 5845 mnts (+3.5/day) 2016-08-22 9000 pkgs (+5...

...n lme4 syntax so I can follow what is happening in the text I'd be grateful. Cheers Paul Cossens #Pixel # P&B(2000) p40-45 Pixel<-read.csv("Pixel.csv",header=TRUE); Pixel$Side<-as.factor(Pixel$Side) Pixel$Dog<-as.factor(Pixel$Dog) (fm1Pixel <- lmer(pixel ~ day + I(day^2) +(day|Dog)+(1|Side), data = Pixel)) (fm2Pixel <- lmer(pixel ~ day + I(day^2) +(day|Dog), data = Pixel)) (fm3Pixel <- lmer(pixel ~ day + I(day^2) +(1|Dog:Side), data = Pixel)) or should I do it this way? Pixel$DS<-with(Pixel,Dog:Side)[drop=TRUE] (fm3Pixel <- lmer(pixel ~ day +...

How many of you are using Rails: 1. As the primary development tool at your Day Job 2. As a small part of your Day Job 3. Not using Rails in your Day Job, but are using Ruby 4. For OS/outside work/hobby -- sam http://www.magpiebrain.com/

I know that this is correct: library(chron) x = dates("01-03-04", format="d-m-y", out.format="day mon year") print(x) It gives me the string "01 Mar 2004" which is correct. I also know that I can say: print(day.of.week(3,1,2004)) in which case he says 1, for today is monday. My question is: How do I combine these two!? :-) I have a data file which is being parsed n...

...------------------------- When data imported from Excel 2007 into R (2.13) all time variables, depending on date (summer or winter) get (un-asked for it!) a time zone addition CEST (for summer dates) or CET (for winter dates). > Dataset Start End1 End2 days2End1.from.Excel days2End2.from.Excel days2End1.in.R days2End2.in.R 1 2010-01-01 2011-01-01 2012-01-01 365 730 365 days 730.0000 days 2 2010-02-01 2011-02-01 2012-01-01 365...

Hello! I have a csv file of intra-day financial data (5-min closing prices) that looks like this: (obs - the dates are formated as day/month/year, as is usual here in Brazil) Date;Time;Close 01/09/2009;10:00;56567 01/09/2009;10:05;56463 01/09/2009;10:10;56370 ##(goes on all day) 01/09/2009;16:45;55771 01/09/2009;16:50;55823 01/09/2009...

Continuing the tradition to post millennia milestones on CRAN: So, it happened. Today (January 27, 2017 PCT) CRAN reached 10,000 packages [1]. Needless to say, the rate with which new packages are added to CRAN keeps increasing and so does the number of contributors (maintainers). Somewhere out there, there are ~3 persons who are about to submit their first packages to CRAN today a...

All, I have: sales <- c(2300,900,155,102,42,10) names(sales) <- c("Christmas","Valentine's Day", "Mother's Day","Father's Day", "Thanksgiving","New Year's Day") barplot(sales,ylim=c(0,2500)) But it doesn't place all of the name labels on the plot. So I tried: sales <- c(2300,900,155,102,42,10) names(sales...

...000(PARTIAL) 33.8 at interval beginning 19720517 144507 19730629 080000(PARTIAL) 32.2 at interval beginning 19720910 135747   and so on...   The original data set are given below:       Year Ending           Maximum Value (Rain mm)  Gap from 19710509 090000 to 19710607 080000 of   28.96 Days 19710629 080000(PARTIAL)       39.3 at interval beginning 19701126 010326  Gap from 19710809 110221 to 19710816 080000 of    6.87 Days  Gap from 19720207 070000 to 19720303 080000 of   25.04 Days  Gap from 19720424 070000 to 19720504 080000 of   10.04 Days 19720629 080000(PARTIAL)       33.8 at in...

...ping question however, given my limited exposure to R, I have not found a solution yet despite my research efforts and wild attempts at what I thought "might" produce some sort of result. I have a very simple list of integers that range between 1 and 24. These correspond to hours of the day. I am trying to create a grouping of Day and Night with Day = 6 to 17.99 Night = 1 to 5.59 and 18 to 24 Using the Cut command I can create the segments but I have not found a "combine" type of command to merger the two "night" segments. No luck with if/else either. Any he...

...----------------------------- When data imported from Excel 2007 into R (2.13) all time variables, depending on date (summer or winter) get (un-asked for it!) a time zone addition CEST (for summer dates) or CET (for winter dates). > Dataset Start End1 End2 days2End1.from.Excel days2End2.from.Excel days2End1.in.R days2End2.in.R 1 2010-01-01 2011-01-01 2012-01-01 365 730 365 days 730.0000 days 2 2010-02-01 2011-02-01 2012-01-01 365...

Dear list member, My question is related to input file format to an Anova from car package. Here is an example of what I did: My file format is like this (and I dislike the idea that I will need to recode it): Hormone day Block Treatment Plant Diameter High N.Leaves SH 23 1 1 1 3.19 25.3 2 SH 23 1 1 2 3.42 5.5 1 SH 23 1 2 1 2.19 5.2 2 SH 23 1 2 2 2.17 7.6 2 CH 23 1 1 1 3.64 6.5 2 CH 23 1 1 2 2.8 3.7 2 CH 23 1 2 1 3.28 4 2 CH 23 1 2 2 2.82 5.2 2 SH 23 2 1 1 2.87 6.4 2 SH 23 2 1 2 2.8 6 2 SH 23 2 2 1 2.02 4.5 2 SH 23...

I have a set of data with 12 readings for temperature per day, with 180 days. I want to find the average temperature of each day. I am able to do this one by one, but with that many days to calculate the average for, it will get very long. I'm sure there is a faster way to do this, I just don't know how. What i have so far is: av1 <- sub...

...randomly divided into 3 groups ("c", "t1", "t2"). There is no need for inclusion of the random effect of the intercept, because the nine sample units are homogeneous among each other (samples taken from the same culture). The algal concentration was measured every two days for 10 days. The goal is to test differences between treatments. I estimated a model with only the intercept and the interaction Group * Day to test which is the best: library(nlme) library(lme4) Day = rep(c(0,2,4,6,8,10),each=9) Group = rep(c("c","c","c","t1&q...