similar to: named arguments in formula and terms

Dear Achim, >>>>> Achim Zeileis <Achim.Zeileis at r-project.org> >>>>> on Fri, 10 Mar 2017 15:02:38 +0100 writes: > 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

On Mon, 4 Apr 2016, varin sacha via R-help wrote: > Hi Deepak, > > In econometrics there is another test very often used : the white test. > The white test is based on the comparison of the estimated variances of > residuals when the model is estimated by OLS under the assumption of > homoscedasticity and when the model is estimated by OLS under the > assumption of

JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 1968 NA 10.7 10.0 9.3 7.4 8.1 9.3 9.5 8.5 10.0 10.0 13.0 1969 13.0 9.9 7.0 5.9 NA 6.5 7.3 6.6 NA NA NA NA 1970 10.7 8.9 8.1 NA NA 7.1 7.7 NA 6.5 NA 8.3 NA 1971 NA 11.6 NA NA 7.8 NA NA 6.2 NA NA 8.4 NA x<-read.zoo("textconnection", sep=",", format= "%Y", header=TRUE) plot(x) #this plots all of the years by each

A pragmatic solution could be to create a simple linear regression example with variables in the global environment and then another example with a data.frame. The latter might be somewhat more complex, e.g., with several regressors and/or mixed categorical and numeric covariates to illustrate how regression and analysis of (co-)variance can be combined. I like to use MASS's whiteside

Dear R users, I want to estimate a Cox PH model with time-dependent covariates so I am using a counting process format with the following formula: Surv(data$start, data$stop, data$event.time) ~ cluster(data$id) + G1 + G2 + G3 + G4 + G5 +G6 Gs represent a B-spline basis functions so they sum to 1 and I can't estimate the model as is without getting the last coefficient to be NA, which

Hi, everybody. I am working processing EEG data from 1000 pacients. I have a specific syntax to perform the Spectral Analysis and a loop to analyse all subjects. each subject data are in separate folders (P1, P2 P3...) My question is: in some cases, some errors can appear in one subject. I want to know if is possible to jump to the next subject and perform the same syntax , exibiting an error

On Sat, 15 Dec 2018, frederik at ofb.net wrote: > I agree with Steve and Achim that we should keep some examples with no > data frame. That's Objectively Simpler, whether or not it leads to > clutter in the wrong hands. As Steve points out, we have attach() > which is an excellent language feature - not to mention with(). Just for the record: Personally, I wouldn't recommend

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! I noticed that there is a (minor) bug either the command all.equal() or in the "plm" package. I demonstrate this using an example taken from the documentation of plm(): ====================================== R> data("Produc", package="plm") R> zz <- plm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp, + data=Produc,

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

Hi! I have three models. In the first model, everything is fine. However, in the second and third models, I have NA's for standard errors: The hessians also have NaN's (same for m2 and m3). What should I do about it? It there a way to obtain the hessian without transforming my variables? I will greatly appreciate your help! -- View this message in context:

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?

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