similar to: new behavior in model.response

On Thu, 28 Jun 2018, Therneau, Terry M., Ph.D. via R-devel wrote: > I now understand the issue, which leads to a different and deeper issue > which is "how to assign a proper length to Surv objects". > > > Surv(c(1,2,3), c(1,0,1)) > [1] 1? 2+ 3 > > The above prints as 3 elements and is conceptually 3 elements. But if I > give it length method to return a 3

To the R development team; I found an unusual behavior in zoo when I upgraded to R 2.11.0 - it abruptly terminated when I performed certain operations on large zoo objects. I sent an e-mail to Achim Zeileis and he said this was a potential bug that I should report to the R development team. The details are given below in the thread below. Basically, I can crash R with this code: library(zoo)

I want to estimate the following fixed effect model: y_i,t = alpha_i + beta_1 x1_t + beta_2 x2_i,tx2_i,t = gamma_i + gamma_1 x1_t + gamma_2 Z1_i + gamma_3 Z2_i I can use ivreg from AER to do the iv regression. fm <- ivreg(y_i,t ~ x1_t + x2_i,t | x1_t + Z1_i + Z2_i, data = DataSet) But, I'm not sure how can I add the fixed effects. Thanks! [[alternative HTML

Hi everyone, I wanted to draw your attention to a new post on the developer.R-project.org blog: https://developer.R-project.org/Blog/public/2019/04/01/hcl-based-color-palettes-in-grdevices/ A new function grDevices::hcl.colors() greatly extends the color palette functionality available in base R. Also, the defaults in the heatmap functions image() and filled.contour() have been adapted

On Fri, 7 Apr 2006 16:09:53 +0200 Brandt, T. (Tobias) wrote: > Hi > > Is there a reason why as.matrix.zoo doesn't set the rownames to the > index like as.data.frame.zoo does? Yes, oversight, I think ;-) Added to the devel-version of zoo. I'll try to get it out to CRAN in the next days. Z > Thanks > > Tobias > > > > library(zoo) > > x.Date

It's two days I try to answer to a mail of Achim Zeileis, the maintainer of strucchange package, using the mail address: zeileis at ci.tuwien.ac.at. I got the following error: ----- The following addresses had permanent fatal errors ----- <zeileis at ci.tuwien.ac.at> (reason: 550 5.7.1 <zeileis at ci.tuwien.ac.at>... Access denied) Can anyone tell me how to reach him?

At the risk of ridicule for my deficient linear algebra skills, I ask for help using the solve.QP function to do portfolio optimization. I am trying to following a textbook example and need help converting the problem into the format required by solve.QP. Below is my sample code if anyone is willing to go through it. This problem will not solve because it is not set up properly. I hope I

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

On Thu, Jan 22, 2015 at 12:45 PM, Achim Zeileis <Achim.Zeileis at uibk.ac.at> wrote: > On Thu, 22 Jan 2015, Max Kuhn wrote: > >> I've had a lot of requests for additions to the reproducible research >> task view that fall into a grey area (to me at least). >> >> For example, roxygen2 is a tool that broadly enable reproducibility >> but I see it more as

Dear all, I was intrigued by the ctree command and wanted to check it out. I first ran the demo with example(ctree) and did get the survival graphs in the end. Upon doing this with my own data and yielding a "Invalid operation on a survival time" I tried to rerun example(ctree) and now I also get "Invalid operation on a survival time" after the example runs plot(GBSG2ct)...

I have NEWS, WISHLIST and THANKS files in the 'dyn' package in the same directory as the DESCRIPTION file but I noticed that they did not survive the move to CRAN, at least on Windows. How do I incorporate them so that they are not omitted?

Hi, we came across the following unexpected (for us) behavior in terms.formula: When determining whether a term is duplicated, only the order of the arguments in function calls seems to be checked but not their names. Thus the terms f(x, a = z) and f(x, b = z) are deemed to be duplicated and one of the terms is thus dropped. R> attr(terms(y ~ f(x, a = z) + f(x, b = z)),

Hello, I need to classify (i.e., export a vector with terminal node id's) new cases using a ctree (party package) model based on different cases (learning sample). I tried the where command with the following syntax: > where(tree, newdata=data2) expecting to get terminal nodes of data2 cases based on rules of tree model (data1 as learning sample). However it returned the following error

On Thu, Jan 22, 2015 at 1:05 PM, Achim Zeileis <Achim.Zeileis at uibk.ac.at> wrote: > On Thu, 22 Jan 2015, Max Kuhn wrote: > >> On Thu, Jan 22, 2015 at 12:45 PM, Achim Zeileis >> <Achim.Zeileis at uibk.ac.at> wrote: >>> >>> On Thu, 22 Jan 2015, Max Kuhn wrote: >>> >>>> I've had a lot of requests for additions to the

Hi everyone and Achim, Achim, I appreciate your help about the function "NCOL". When I use "NCOL" instead of "ncol", I can find out the number of columns (number of time series) in the presence of only one time series (one variable, one column). Now I want to know how I can find out the column names of the zoo objects? In case of more than one time series, the

Dear Zeileis: Thank you very much abou On Tue, Sep 19, 2017 at 4:13 AM, Achim Zeileis <Achim.Zeileis at uibk.ac.at> wrote: > > > On Tue, 19 Sep 2017, AbouEl-Makarim Aboueissa wrote: > > Dear All: good morning >> >> I am trying to graph the function y=f(x)=1/x over the interval (-5,5). But >> I am getting an error message. Please see below. >>

On Thu, 23 Feb 2023, Duncan Murdoch wrote: > On 23/02/2023 6:09 a.m., Achim Zeileis wrote: >> Duncan, >> >> thanks for your feedback. I just received your response after sending out >> mine. You came to the same conclusion. Should I prepare a patch and send >> it to you so that you can also have a look? Or view Bugzilla? > > Copying Sigbert's message to

Dear Achim, Thanks for the comment. Here is the publisher's style guideline (AGU) "Because the Internet is dynamic environment and sites may change or move, treat World Wide Web, ftp files, and electronically archived data stored at data centers other than World or National Data Centers as unpublished, i.e., in text only." http://www.agu.org/pubs/AuthorRefSheet.pdf So by this

Hi,   A question on something which normally should be easy !   I perform a linear regression using lm function:   > reg1 <- lm (a b+c+d, data = database1)   Then I try to perform the Chow (1960) test (structural change test) on my regression. I know the breakpoint date. I try the following code like it is described in the “Examples” section of the “strucchange” package :   > sctest(reg1,

Dear useRs, the program for the first R user conference useR! 2004 (taking place in Vienna, May 20-22) is now available online at the conference web page: http://www.ci.tuwien.ac.at/Conferences/useR-2004/program.html We received many exciting abstracts about a wide spectrum of applications in which R is used. This is reflected in the program and I think we can look forward to a very