similar to: tail.matrix returns matrix, while tail.mts return vector

Hi My goal is to do a (multiple) regression, just knowing that my Y variables will be the say k first variables of a matrix/data frame. I thought I should do it with eval(parse)) but encounter a strange problem. See: lm(y~.-y, data=freeny) #that's what I want to do in the one equation case #Problem is I don't know name of the variable... only that it is the first one... #so idea is to

Hi Say I want to add manually an intercept in the function lm. Even if almost all results will be identical, few stats are different as DF counting will be different as intercept will not be included in "automatic" case, while it will be in "manual" case. See: ###usual lm on freeny fr<-lm(freeny.y~freeny.x) ###manual lm on freeny man<-cbind(1,freeny.x)

can you recommend a good manual for R that starts with a data set and gives demonstrations on what can be done using R? I downloadedR Langauage definition and An introduction to R but haven't found them overly useful. I'd really like to be able to follow some tutorials using a dataset or many datasets. The datasets I have available on R are Data sets in package 'datasets':

Hello! I would love to be able to include an external variable to a lm call, I mean something: if(TRUE) a<-freeny.x[,4] else a<-NULL lm(freeny.y~freeny.x[,-4] +a) but it does not work with a<-NULL, whereas lm(freeny.y~freeny.x[,-4] +NULL) I don't understand why and did not find an answer in the manuals... do you see it? Any idea? Thanks!!

Hi, A recent thread provided a (working) construct for lm: lm(as.matrix(freeny[ix]) ~., freeny[-ix]) Can someone explain what is meant by the formula in that expression, that is, what does "mymatrix~." do? I couldn't find any such example in the lm() or formula() help pages. thanks Carl

I am getting the following: *** Weave Errors *** Error in driver$runcode(drobj, chunk, chunkopts) : Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" *** Source Errors *** Error in eval(expr, envir, enclos) : couldn't find function "pacf.mts" make[1]: *** [checkVignettes] Error 1 I don't really understand the new namespace mechanism,

There appears to be a bug in plot.mts() in R-1.3.0 when attempting to do a multiple point plot (this worked in 1.2.2). Reproduce by: > version _ platform sparc-sun-solaris2.8 arch sparc os solaris2.8 system sparc, solaris2.8 status major 1 minor 3.0 year 2001

Dear all, it looks like fix() changes the class of mts objects, here is a reproducible example (tested both on WinXP and Linux): > x <- ts(cbind(1:5,1:5)) > x Time Series: Start = 1 End = 5 Frequency = 1 Series 1 Series 2 1 1 1 2 2 2 3 3 3 4 4 4 5 5 5 > class(x) [1] "mts" "ts" > edit(x)

Full thread here: https://github.com/tidyverse/broom/issues/287 Reproducible example: is.object(freeny$y) # [1] TRUE attr(freeny$y, 'class') # [1] "ts" class(freeny$y) # [1] "ts" # ts attribute wiped by model.frame class(model.frame(y ~ ., data = freeny)$y) # [1] "numeric" attr(model.frame(y ~ ., data = freeny)$y, 'class') # NULL # but still:

This one showed up while looking at one of Ripley's other reports: > data(freeny) > model.frame(y~1,data=freeny,subset=1:10) y 1962.25 8.79236 1962.5 8.79137 1962.75 8.81486 1963 8.81301 1963.25 8.90751 1963.5 8.93673 1963.75 8.96161 1964 8.96044 1964.25 9.00868 1964.5 9.03049 > model.frame(y~1,data=freeny,subset=1:10)$y Warning: Replacement length not a

Hi! I am new to R. I need your help. I have got the time series of fifteen variables in data file. I would like to plot it in R in separate ps pages, not in same ps page. I was reading about mts, but I could not figure out how to do it. Can anyone help me out? with regards; Santosh -------------------------------------------------------------------------- Santosh Kumar

MTS is a file extension for an AVCHD video clip format for high-definition video, s supporting 1080i (a high definition video format with 1080 horizontal scan lines, interlaced) and 720p (720 horizontal scan lines, progressive scan, rather than interlaced) in a relatively small file size. Recently more and more people have updated their Mac computer to Mac OS X 10.9 Mavericks, but video

Panasonic HC-X800 is expert HD camcorder with 3MOS system Pro, shooting 1080 60p/50p AVCHD footage at bitrates of up to 28Mbit/sec and exhibiting outstanding HD images with minimum noise. The HC-X800 records AVCHD MTS videos on P2 card. You may want to edit AVCHD MTS videos you shoot with iMovie. The problem is that iMovie often cannot import and recognize the 60p/50p footages, nor can it accept

Dear all, I have discovered an object named 'freeny.x' which seems in-built. Can somebody tell me why it is there and what purpose it is serving? 'ls()' syntax does not seem to catch this object in the environment. Thanks,

I'm trying to use garchFit from fSeries, with Student or Skewed Student conditionnal distribution. Let's say that eps (vector) is my series of daily log-returns: data(EuStockMarkets) eps = diff(log(EuStockMarkets[,"CAC"])) library(fSeries) g = garchFit(series = eps, formula.var = ~garch(2,2), cond.dist = "dstd") s = g at fit$series All the coefficients are ok

Dear R-help, Has anyone implemented sparklines in the strips of a lattice plot? What I have in mind is, say, highlighting that part of a time series that one is examining in more detail in a set of lattice plots. Regads,. Mark Difford. PS: (Andreas Loffler has implemented a simple but functional version for TeX/LaTeX: http://www.tug.org/tex-archive/help/Catalogue/entries/sparklines.html)

Hi, I’d like to fit one GLM model for each possible combination of inputs (i.e. exhaustive search). The package leaps can help me to generate all possible variable subsets, but I’ll appreciate your guidance as of how to generate one model for each of those possible subsets. I’m new in R! Thanks in advance for your help. library(leaps) subsets <- regsubsets( y ~ . , data=freeny, nbest=100,

Dear Sir, ???????????????? I am unable to work MTS MBLAZE ULTRA in CentOS 6.4. Can any one suggest a way to active MTS MBLAZE ULTRA in CentOS 6.4 properly? ????????????????????????????????????????????????????????????????????????????????????????????????? Yours faithfully ????????????????????????????????????????????????????????????????????????????????????? ?? ASOK UMAR MALLICK

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

Hi useRs, I am trying to calculate the Sharpe ratio with "sharpe" of the library "tseries". The documentation requires the univariate time series to be a portfolio's cumulated returns. In this case, the example given data(EuStockMarkets) dax <- log(EuStockMarkets[,"FTSE"]) is however not the cumulated returns but rather the daily returns of the FTSE