Displaying 20 results from an estimated 500 matches similar to: "substitution in a function"
2010 Jul 14
1
Converting POSIXct vales to real values
I have a dataframe that contains time values in the form of yyyy-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 tvec/(60*60*24) I only get integer values and not fraction of a day which is what I want.
What do I need to
2001 Oct 26
3
question about anova() output
Hello,
I am getting output from anova() and summary(aov()) that depends on the
order of the factors in the fitted model object, and this has me baffled. I
see this dependency with the data.frame below but not with an example (table
6.4) from Montgomery's DOE book. This is with R 1.3.0 on Debian GNU-Linux.
Where have I gone wrong?
> centerpts
run sample CH50mg
1 day1 dev126 0.56
2
2010 Aug 12
0
termplot for mixed linear effects models
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) (08) 8595 2209
*(mob) 0401 124 971
*(fax) (08) 8595 2232
* 28 Vaughan Terrace,
2006 May 07
1
Anyone care for a braindump?
I have this problem with records that have to be alligned on end- and
startdate. I came up with the following. but am not convinced that this
is the best way to tackle this problem. Anyone care for a braindump?
def head_to_tail
# Remove days that have startdate >= self.startdate AND enddate <=
self.enddate
#
# before after
#
# |=====| |=====|
2010 Jul 15
2
taking daily means from hourly data
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 16 3.268 3.272 3.240 1441.8 192.338 5.496
3 2009-07-06 17:00:00 17 3.268
2008 Jul 31
2
S 3 generic method consistency warning please help
I would like to include this in a package. The S3 methods on R CMD check
says
* checking S3 generic/method consistency ... WARNING
window:
function(x, ...)
window.chron:
function(data, day1, hour1, day2, hour2, ...)
See section 'Generic functions and methods' of the 'Writing R Extensions'
manual.
I have looked and can not figure it out. This function is for convience.
What
2009 Jan 21
1
finding row and column indices of date in multiple columns of a data frame
Hi,
I have a data.frame SAMPLES with columns:
Site Site# Season Day1 Day2 Day3
Day1, Day2, Day3 are class "Date", the other columns are numeric or
factor.
I have a date "mydate" that may or may not be listed in my data.frame
and I need to find that out.
If "mydate" is there, I want to get the number of the data.frame row
where it occurs.
2009 Sep 11
1
help with plotting
HI all,
raw_urine =
read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\mz.spot.sam.dat.new",
header = TRUE )
pvalue =
read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\all.urine.features.t.test.result",
header = TRUE )
library(compositions)
p = function(a,b){
y = pvalue[,a]
if(y<0.01){
index = which(y, arr.ind=TRUE)
day1 = raw_urine[index,3:7]
day2 =
2008 Jul 22
0
loop for multiple regressions
Dear all,
I have the following data in excel:
day1 y 1 2 3 2 3 x1 0.2 0.3 0.4 0.3 0.2 x2 7 3.4 2 8 6 day2 y 2 4 3
2 2 x1 0.4 0.5 0.3 0.3 0.2 x2 7 8 9.1 6 5
I have the following problems:
first of all, when I ask R to read the file (with the package xlsReadWrite
and the command read.xls) it has a problem with the fact that the left most
corner is labelled the same way, so in order for it
2010 May 21
4
indexing problem
Dear group,
Here is my environment :
> ls()
[1] "l" "PLglobal" "Pos100415" "Pos100416" "Pos100419"
"Pos100420" "position" "select" "Trad100415" "Trad100416"
"Trad100419" "Trad100420" "trade" "y"
With objects :
> l
[1]
2007 Aug 17
1
finding the row(s) for a date in a data frame
Hi,
If I have a data frame A with the following format:
Day1 Day2 Day3 Day4
1 1979-11-02 1979-11-03 1979-11-04 <NA>
2 1979-12-06 <NA> <NA> <NA>
3 1979-12-13 1979-12-14 1979-12-15 1979-12-16
4 1979-12-20 <NA> <NA> <NA>
And a date "1979-12-14", for
2010 Sep 16
2
standard error of difference for mixed effects
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 difference at several points. I am not sure of how to calculate the difference and which estimate of se I
2012 Nov 01
2
Name assignment in for loop
Dear helpeRs-
I'm using a for loop to create a series of models.
I'm trying to assign a name to each model created,
using the loop index. The loop gets stuck at the name
of the model, giving the error "target of assignment
expands to non-language object". The linear model runs
without error; only the name is problematic.
Here is the current loop syntax. The use of dat
2009 Jun 01
1
Reshaping Data
Hi,
i did a mistake with my first post.
I have to reshape data from this matrix:
id x1 x2 x3 x4 day1 day2 day3 day4 day5 day6 day7 day8 day9
1 0.129 0.797 0.231 0.615 4 4 1 1 1 1 3 3 3
2 0.420 0.376 0.501 0.282 4 4 4 4 5 4 2 5 5
3 0.377 0.486
2006 May 31
0
Ruby on Rails Workshop in Kansas City !
Ruby on Rails Workshop
------------------------------------
Kansas City Ruby User''s Group (http://kcrug.org/) proudly presents the first "Ruby on Rails Workshop" for those eager to get started into Ruby on Rails. This one-and-half-day workshop will offer a quick introduction to Ruby and the Ruby on Rails web application framework.
Sponsored by Reevik, Inc., it will
2012 Mar 15
1
Bar graph with 2 Y axis
Dear R users,
I need to draw a barplot with 2 Y axis. I have 3 days each of wich having 2 groups (and error bar for each of them). The height of the 3rd day is too tall compared to others. That's why I have to use a second Y axis for that. I am using "barplot2" function of "gplots" library (to be able to add error bars as well). Data and codes currently I am using is
2002 Nov 01
1
Reshape function
Can someone help me with the proper usage of the
reshape function?
Let's say I have a dataset with columns like this
(wide format):
Id Sex Group Test Day1 Day 2 Day 3...
And I want to transpose this into something like this
(thin):
Id Sex Group Test Time
Where the new column labeled time contains all the
time variables (Day 1, Day 2, Day 3...) that were in
the wide
2011 Nov 01
1
Counting entries to create a new table
Hi,
I am an R novice and I am trying to do something that it seems should be fairly simple, but I can't quite figure it out and I must not be using the right words when I search for answers.
I have a dataset with a number of individuals and observations for each day (7 possible codes plus missing data)
So it looks something like this
Individual A, B, C, D
Day1 1,1,1,1
Day 2 1,3,4,2
Day3
2006 Jan 25
0
Interpolating spline problems and akima
Hi everyone
I was using spline to interpolate single or two consecutive missing data points
in time series. However, when it comes to longer gaps in the data the spline
function generate new data for both my known and unknown data (see below).
Aside from not understanding why this happens, I thought thought I might try
function "aspline" in library (akima). However, I cannot install or
2013 Jan 17
1
plotting from dataframes
thanks to your guys help I am closer to solving my problem but I have some
small problem. So let's say I start with
>data
number day hour
1 17 10
2 17 11
3 17 6
4 18 4
5 18 10
6 19 8
7 19 8
I want to split to odd days, which I am able to do, I call this object
frames, which looks like:
> frames
$`1`
c1 day1 hour1
1 1 17 10
2 2 17 11
3 3 17 6
$`2`
c1 day1