[R] Problem with augPred?

2010-09-27 Thread Meissner, Tony (DFW)
I have a dataset, example of the data is shown below: Grouped Data: drain_irr ~ irr | crop Year decades crop irrisystem drain_irrirr 1310 1995-96 1990s Citrusvarious 0.400 9.021 1311 1995-96 1990s Citrus drip 0.541 6.468 1312 1995-96 1990s Citrus overhead

[R] Problem with augPred?

2010-09-27 Thread Meissner, Tony (DFW)
I have a dataset, example of the data is shown below: Grouped Data: drain_irr ~ irr | crop Year decades crop irrisystem drain_irrirr 1310 1995-96 1990s Citrusvarious 0.400 9.021 1311 1995-96 1990s Citrus drip 0.541 6.468 1312 1995-96 1990s Citrus

[R] standard error of difference for mixed effects

2010-09-15 Thread Meissner, Tony (DFW)
I have a dataset relating the effects of engineering works on the level of salinity in a river before and after the works. I have modelled this using linear mixed effects models to determine if the significance and level of the response to the works. I am wanting to calculate the se of

[R] termplot for mixed linear effects models

2010-08-12 Thread Meissner, Tony (DFW)
Is there an equivalent package for mixed linear effects models developed using the package nlme as there is for linear models? Tschüß Tony Meissner Principal Scientist (Monitoring/Statistics) Resource Monitoring Science, Monitoring and Information Division Department for Water Imagine © *(ph)

Re: [R] POSIXlt error

2010-07-14 Thread Meissner, Tony (DFW)
OK, I thought I was being explicit enough, but here goes. The out of the first five rows of mu dataframe is: Time hour lev.morgan lev.lock2 lev.lock1 flow direction velocity bearingvelcat 1 2009-07-06 15:00:00 15 3.266 3.274 3.240 1710.6180.28 4.352

Re: [R] POSIXlt error -Oops

2010-07-14 Thread Meissner, Tony (DFW)
OK, I thought I was being explicit enough, but here goes. The out of the first five rows of mu dataframe is: Time hour lev.morgan lev.lock2 lev.lock1 flow direction velocity bearingvelcat 1 2009-07-06 15:00:00 15 3.266 3.274 3.240 1710.6180.28 4.352

[R] taking daily means from hourly data

2010-07-14 Thread Meissner, Tony (DFW)
I have a data frame (morgan) of hourly river flow, river levels and wind direction and speed thus: Time hour lev.morgan lev.lock2 lev.lock1 flow direction velocity 1 2009-07-06 15:00:00 15 3.266 3.274 3.240 1710.6 180.282 4.352 2 2009-07-06 16:00:00

Re: [R] taking daily means from hourly data

2010-07-14 Thread Meissner, Tony (DFW)
is prohibited and may be unlawful. If you have received this e-mail in error, please advise by return e-mail or by telephoning +61 8 8595 2209 From: Allan Engelhardt [mailto:all...@cybaea.com] Sent: Thursday, 15 July 2010 3:09 PM To: Meissner, Tony (DFW) Cc: r-help@r-project.org Subject: Re

[R] Converting POSIXct vales to real values

2010-07-13 Thread Meissner, Tony (DFW)
I have a dataframe that contains time values in the form of -mm-dd hh:mm i.e. 2010-07-14 13:00. When I convert this to numeric via tvec - as.numeric(Time) I get a number that is in seconds. So far so good. When I then divide the numeric value by the number of seconds in a day eg