Displaying 9 results from an estimated 9 matches for "ddat".
Did you mean:
dat
2010 Jun 14
1
recursively Merging a list a zoo objects
The zoo package as a merge function which merges a set of zoo objects
result<-merge(zoo1,zoo2,...)
Assume your zoo objects are already collected in a list
# make a phony list to illustrate the situation. ( hat tip to david W for
constructing a list in a loop)
ddat <- as.list(rep("", 20))
ytd<-seq(3,14)
for(i in 1:20) {
+ ddat[[i]] <- zoo(data,ytd )
+ }
ddat
[[1]]
1 2 3 4 5 6 7 8 9 10 11 12
3 4 5 6 7 8 9 10 11 12 13 14
[[2]]
1 2 3 4 5 6 7 8 9 10 11 12
3 4 5 6 7 8 9 10 11 12 13 14
[[3]]
1 2 3 4 5 6...
2011 Jun 13
1
maintaining row connections during aggregate
...is a value for each day of multiple years. In this particular data set it goes up to 1974. I am am looking to obtain the minimum and maximum values for each year, but also know on which julian day ("jday") they occurred.
I can get the maximum value for each year with:
> mx = aggregate(ddat$avg_m3s, list(Year=ddat$year), max, na.rm=T)
> colnames(mx) <- c("year","max_daily")
year max_daily
1 1960 60.24615
2 1961 73.90000
3 1962 56.40000
...
But I want to output the max with the corresponding day on which it occurred, such as:
year jday avg_m3s
1...
2017 Aug 19
4
My very first loop!! I failed. May I have some start-up aid?
...%H:%M"), tz="GMT")
mysamplexy <- project (cbind (mysample$Longitude, mysample$Latitude),
"+proj=utm +zone=32 +ellps=WGS84")
colnames(mysamplexy) <- c ("xCord", "yCord")
ID <- mysample$ID
datltraj <- as.ltraj(mysamplexy, mysampletime, id=ID)
Ddat <- BRB.D(datltraj, Tmax=21600, Lmin=36)
BRBdat <- BRB(datltraj, D= Ddat,type=c("UD"),Tmax=21600,Lmin=36, hmin=100)
kernel.area(BRBdat, unout=c("km2"))
# unfortunately my data are not a very good example. Sorry about that.
But I guess you know what I mean.
# Because I wi...
2017 Aug 19
0
My very first loop!! I failed. May I have some start-up aid?
...]
)
, "+proj=utm +zone=32 +ellps=WGS84"
)
colnames( mysamplexy ) <- c( "xCord", "yCord" )
datltraj <- as.ltraj( mysamplexy, mysample$mysampletime, id=mysample$ID )
Ddat <- BRB.D( datltraj, Tmax=21600, Lmin=36 )
BRBdat <- BRB( datltraj, D = Ddat, type = "UD", Tmax = 21600, Lmin = 36, hmin = 100 )
result[ , i ] <- kernel.area( BRBdat, unout = "km2" )
}
#> Warning in kernel.area(BRBdat, unout = "km2"): The grid is too s...
2009 Jan 22
1
Problem with cex=0.1 when making jpegs
I am using the following script to make .jpg files.
jpeg('plotx.jpg')
ddat <-read.table("file",header=T)
attach(ddat)
tdat<-read.table("file1")
plot(xx1,yy1,type='p',pch=1,col="blue",cex=0.2,xlim=c(0,3.5),ylim=c(-75,75))
points(tdat,col="green",pch=1,cex=0.2)
dev.off()
The problem is that I want the points to be very...
2017 Aug 19
0
My very first loop!! I failed. May I have some start-up aid?
...uot;)
> mysamplexy <- project (cbind (mysample$Longitude, mysample$Latitude),
> "+proj=utm +zone=32 +ellps=WGS84")
> colnames(mysamplexy) <- c ("xCord", "yCord")
> ID <- mysample$ID
> datltraj <- as.ltraj(mysamplexy, mysampletime, id=ID)
> Ddat <- BRB.D(datltraj, Tmax=21600, Lmin=36)
> BRBdat <- BRB(datltraj, D= Ddat,type=c("UD"),Tmax=21600,Lmin=36, hmin=100)
> kernel.area(BRBdat, unout=c("km2"))
>
> # unfortunately my data are not a very good example. Sorry about that. But
> I guess you know what...
2008 Oct 21
5
how to plot the histogram and the curve in the same graph
i want to plot the histogram and the curve in the same graph.if i have a set
of data ,i plot the histogram and also want to see what distribution it
was.So i want to plot the curve to know what distribution it like.
--
View this message in context: http://www.nabble.com/how-to-plot-the-histogram-and-the-curve-in--the-same-graph-tp20082506p20082506.html
Sent from the R help mailing list archive at
2012 Aug 05
4
find date between two other dates
Hi,
I am trying to assign "Season" values to dates depending on when they occur.
For example, the following dates would be assigned the following "Season"
numbers based on the "season" intervals detailed below in the code:
ddate Season
29/12/1998 20:00:33 1
02/01/1999 05:20:44 2
02/01/1999 06:18:36 2
02/02/1999 07:06:59 3
02/03/1999 07:10:56 4
02/03/1999 07:57:18 4
My approach so far doesnt work because of the time stamps and is probably
very long winded....
2013 Feb 20
0
Bayesian mixing model
...Na2O Al2O3
Topsoils 1.8552515 0.1135672 0.06212094 1.491125
ChannelBanks 9.8400162 0.1401057 0.08599080 2.708710
FieldDrains 2.3896499 0.1961217 0.02545431 4.300644
RoadRunoff 0.7780579 0.1749869 0.02848264 1.116747", header =TRUE)
...and absolute stdev on the target called adssdDat:
abssdDat<-read.table(text=" CaO MgO Na2O Al2O3
1 0.877 0.531 0.264 2.439", header=TRUE)
Whilst the model results are OK, there is one problem. The model fails to take into consideration covariance between the geochemistry in each source area (Rat) and therefore the uncertaint...