Displaying 20 results from an estimated 2000 matches similar to: "Pointwise division of two zoo objects?"
2008 Jun 18
1
Pointwise Confidence Bounds on Logistic Regression
Hi all. I hope I have my terminology right here...
For a simple lm, one can add ?pointwise confidence bounds? to a fitted line
using something like
>predict(results.lm, newdata = something, interval = "confidence")
(I'm following DAAG page 154-155 for this)
I would like to do the same thing for a glm of the logistic regression type,
for instance, the example in MASS pg
2007 Jun 08
1
pointwise confidence bands or interval values for a non parametric sm.regression
Dear all,
Is there a way to plot / calculate pointwise confidence bands or
interval values for a non parametric regression like sm.regression?
Thank you in advance.
Regards,
Martin
2005 Dec 29
1
use of predict() with confidence/prediction bands
To my understanding, a confidence interval typically covers a single
valued parameter. In contrast, a confidence band covers an entire line
with a band. In regression, it is quite common to construct confidence
and prediction bands. I have found that many people are connecting
individual confidence/prediction interval values produced with
predict(object,sd.fit=T,type="conf/pred") and
2012 Dec 06
1
bootstrap based confidence band
I'm trying to find a bootstrap based confidence band for a linear model.
I have created a data set with X and Y
X=runif(n,-1.25,1.25)
e=rnorm(n,0,1)
Y=exp(3*X)+5*sin((30*X)/(2*pi))+2*e
fit=lm(Y~X)
summary(fit)
I define a bootstrap function named PairedBootstrap which is not listed here. Than I try many ways to find the confidence band. One way is to predict Y using the model I get above for
2013 Jun 05
0
[R-pkgs] bpcp package for pointwise confidence intervals for a survival distribution
Hi all,
I just uploaded a new version of the bpcp package. It calculates confidence intervals for a survival distribution for right-censored data using the newly developed beta product confidence procedure. Previously developed methods can have substantial error rate inflation for the lower limit, especially at the right end of the curves when there are small numbers of events. The bpcp method
2011 Sep 27
2
Coercing a character zoo to a numeric
Dear R-helpers,
It seems to me that a character zoo cannot be coerced to a numeric zoo.
Below is a minimal example. Can someone tell me what I have done wrong?
> z<-zoo(1:4,order.by=1:4)
> coredata(z)<-as.character(coredata(z))
> str(z)
‘zoo’ series from 1 to 4
Data: chr [1:4] "1" "2" "3" "4"
Index: int [1:4] 1 2 3 4
>
2009 Jun 19
1
(FULL) Need help to optimize a piece of code involving zoo objects
(Sorry, sent the message before I finished it)
Hello, everyone
I have a long script that uses zoo objects. In this script I used
simple moving averages and these I can very efficiently calculate with
filter() functions.
Now, I have to use special "exponential" moving averages, and the only
way I could write the code was with a for-loop, which makes everything
extremely slow.
I don't
2009 Jun 19
1
Need help to optimize a piece of code involving zoo objects
Hello, everyone
I have a long script that uses zoo objects. In this script I used
simple moving averages and these I can very efficiently calculate with
filter() functions.
Now, I have to use special "exponential" moving averages, and the only
way I could write the code was with a for-loop, which makes everything
extremely slow.
I don't know how to optimize the code, but I need to
2006 Nov 23
1
Problem with as.ts(zoo-object)
Dear all,
I have an error message, when I try to convert a zoo object (called
test) to ts (on R 2.4.0, Package zoo version 1.2-1, Windows XP)
> test
1994-05-10 1994-06-09 1994-07-09
0.0024943889 0.0024881824 0.0006955831
> str(test)
atomic [1:3] 0.002494 0.002488 0.000696
- attr(*, "index")=Class 'Date' num [1:3] 8895 8925 8955
> is.regular(test)
[1] TRUE
2006 Nov 29
2
problem with indexing a zoo object
My problem is the following : I create 2 zoo objects and then I try to
subset one of them using logic. indicesthatpass is a vector of trues and
falses but when I send it into bckret, it returns an empty bckret.
Obviously it has something to do with bckret being a zoo object and if I
do the same subsctripting off of coredata(bckret), I'm confident it will
work. But, I need to keep the minute
2008 Oct 22
1
R 2.8.0 qqnorm produces error with object of class zoo?
Dear list-reader,
by running the following script:
library(zoo)
sessionInfo()
search()
packageDescription("zoo")
data(EuStockMarkets)
dax <- as.zoo(EuStockMarkets[1:10, "DAX"])
daxr <- diff(log(dax))
identical(as.vector(qnorm(daxr)), qnorm(coredata(daxr)))
qqnorm(coredata(daxr))
qqnorm(daxr)
qqnorm() produces an error:
> qqnorm(daxr)
Fehler in if (xi == xj) 0L
2006 Nov 03
1
as.zoo behavior (
hi all : the code pasted below runs but then, a dput on
rollmeandifflogbidask gives me what is below the code. the structure of
rollmeandifflogbidask is a zoo object but with a "frequency"
so it's not the same structure as the original actual diff and this
really causes things to blow up in later code. i'm sure gabor and achim
know what to do but in the case that they
are not
2012 Jun 06
2
package zoo, function na.spline with option maxgap -> Error: attempt to apply non-function?
Hello,
I'm trying to use na.spline (package zoo) to fill some missing data in a time series.
this works fine, however, if I apply the 'maxgap' argument, I always get the error:
<------
Error in na.spline.vec(x., coredata(object.), xout = xout., ...) : attempt to apply non-function
------>
I couldn't find a similar error for this case in the mailing lists and zoo vignette,
2011 Jan 24
2
how to slice a zoo object
Hi
Would anyone have any pointers on how to slice up a large zoo table. I
have the following structure: -
> str(ZOO_OBJ)
?zoo [1:632, 1:83] 30.4 30.4 30.4 30.4 30.3 ...
?- attr(*, "dimnames")=List of 2
??..$ : NULL
??..$ : chr [1:83] "COL1" "COL2" "COL3" "COL4" ...
?- attr(*, "index")= POSIXct[1:632], format: "2009-05-01
2010 Dec 22
3
How to integrate a function with additional argument being a vector or matrix?
Dear expeRts,
I somehow don't see why the following does not work:
integrand <- function(x, vec, mat, val) 1 # dummy return value
A <- matrix(runif(16), ncol = 4)
u <- c(0.4, 0.1, 0.2, 0.3)
integrand(0.3, u, A, 4)
integrate(integrand, lower = 0, upper = 1, vec = u, mat = A, val = 4)
I would like to integrate a function ("integrand") which gets an "x" value (the
2003 Jun 04
2
gam()
Dear all,
I've now spent a couple of days trying to learn R and, in particular, the
gam() function, and I now have a few questions and reflections regarding
the latter. Maybe these things are implemented in some way that I'm not yet
aware of or have perhaps been decided by the R community to not be what's
wanted. Of course, my lack of complete theoretical understanding of what
2010 May 16
3
Vector recycling and zoo
I am a bit confused about the different approaches taken to recycling in
plain data frames and zoo objects. When carrying out simple arithmetic,
dataframe seem to recycle single arguments, zoo objects do not. Here is an
example
> x <- data.frame(a=1:5*2, b=1:5*3)
> x
a b
1 2 3
2 4 6
3 6 9
4 8 12
5 10 15
> x$a/x$a[1]
[1] 1 2 3 4 5
> x <- zoo(x)
> x$a/x$a[1]
1
1
>
2010 Aug 25
2
Repeat the first day data through all the day. Zoo
down vote favorite
Hello
I have a zoo series. It lasts 10 years and its frequency is 15min.
I'd like to get a new zoo series (or vector) with the same number of
elements, whith each element equal to the first element of the day. That's,
The first element everyday is repeated throughout the wole day.
This is not same as aggregate(originalseries,as.Date,head,1) because this
gives a
2004 Mar 12
1
confidence interval in local polynomial regression
Dear all,
Is there any package or function can do the pointwise confidence interval
and confidence band for the local polynomial regression? Maybe the local
linear regression is enough. Thanks.
Regards,
Zhen
2005 Aug 18
1
display of a loess fitted surface
Good morning,
I am Marta Colombo,student at Politecnico,Milan. I am studying local regression models and I am using loess function. My problem is that when I have a loess object I don't know how to display the fitted surface; in fact, while in S when you have a loess object you can see it writing plot(object), in R this dosen't work. Also I'd like to know if there is something like the