[R] How does R determine best value of alpha in forecasting

2014-01-30 Thread neel shah
Hi I need to implement Holt-Winter's Exponential Smoothing in my project. I was going through implementation of Holt-Winter's Exponential Smoothing in forecast package. I would appreciate if someone can explain how R finds best value of alpha and beta given the dataset. It ll be great if someone c

[R] Change points in R

2014-01-30 Thread Benjamin Ward (ENV)
Hi R helpers, I have a set of data best shown in this below graph. Each coloured line represents a statistic calculated across pairs of DNA sequences. And for each coloured line, I would like to identify breakpoints - so identify the chunks where the values are high, for example, in the light

[R] error in assetsFit function of fPortfolio package

2014-01-30 Thread Mohammad Nikzad
Dear Sir/Madam I get an error when I run a assetsFit function of fPortfolio package. You just need to load the fPortfolio package and run this command to get the same error that I get. The only fitting method that works is "norm" function.the "snorm" and "st" methods give error. library(fPortfol

Re: [R] reversed variables in stats::reshape()

2014-01-30 Thread Doug Morrison
Thanks Arun! I see now that the description of the "varying" argument for reshape includes the following: "This is canonically a list of vectors of variable names" Originally, I saw in the Details section: "Notice that the order of variables in varying is like x.1,y.1,x.2,y.2." Best, Doug On

[R] Generalized Ordered Logit in R

2014-01-30 Thread B. Longhurst
Dear community, I am replicating a paper (The Effect of Structures and Power on State Bargaining Strategies, http://tinyurl.com/oyk289c) for a class project. The author used “gologit” in STATA, and I need the equivalent function in R. She used a “Generalized Ordered Logit” (see below model de

Re: [R] print.matrix

2014-01-30 Thread David Winsemius
On Jan 30, 2014, at 6:46 AM, Göran Broström wrote: > In the documentation of 'prmatrix' (base) I read, under Details: > > ‘prmatrix’ is an earlier form of ‘print.matrix’ > > but 'print' doesn't seem to have a 'matrix' method. And in the 'Examples' > section: > > chm <- matrix(... > chm # uses

Re: [R] Controlling font size on code chunk outputs using Knitr

2014-01-30 Thread Jeff Johnson
Thanks Yihui and Jeff. I've retrieved the default CSS file and made a tweak to it (changing a header 1 size just to test it) and saved it to the same local directory as my .Rmd file using the name 'mymarkdown.css' for testing. I've added: options(rstudio.markdownToHTML = function(inputFile, out

Re: [R] Controlling font size on code chunk outputs using Knitr

2014-01-30 Thread Yihui Xie
Exactly. Please see RStudio documentation: https://support.rstudio.com/hc/en-us/articles/200552186-Customizing-Markdown-Rendering Regards, Yihui -- Yihui Xie Web: http://yihui.name On Thu, Jan 30, 2014 at 10:57 AM, Jeff Newmiller wrote: > This sounds like a classic "you need to write a custom

[R] glmmADMB error

2014-01-30 Thread Sebastián Daza
I ran an example using glmmADMB, and I got this error: data(bacteria,package="MASS") bacteria$present <- as.numeric(bacteria$y)-1 (bfit <- glmmadmb(present ~ trt + I(week > 2), random = ~ 1 | ID, + family = "binomial", data = bacteria)) Error in paste0(symbol1, paste0(paste0(var, collapse = symb

[R] Universal regression program - multiple regressions on set of data

2014-01-30 Thread nooldor
Hi, Here is description of my problem: I have data frame that contains 110 variables (columns) and 595 observations each. Some of the variables will be my Y-dependent variable, some will be Q, X, Z - independent variables I need to estimate robust regression models: Y~X+Z+Q I want to create 4 su

Re: [R] mlogit: message "invalid 'row.names' length" after subsetting data

2014-01-30 Thread craux
Thanks a lot Bill! This works fine! Useful to know regarding subset. Charles -- View this message in context: http://r.789695.n4.nabble.com/mlogit-message-invalid-row-names-length-after-subsetting-data-tp4684465p4684470.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] mlogit: message "invalid 'row.names' length" after subsetting data

2014-01-30 Thread William Dunlap
> dd1 <- subset(dd, dd$Effect=='none') Try using dd1 <- dd[dd$Effect=='none', ] instead of that call to subset(). subset() strips the non-data.frame attributes from its output. It may be easier to use the subset argument to mlogit, as in mlogit(..., subset = Effect=="none"). Bill Dunlap TIBC

[R] mlogit: message "invalid 'row.names' length" after subsetting data

2014-01-30 Thread craux
Hi! I have been working a while with mlogit, estimating successfully a series of models. This time I have a series of discrete choice experiments with different effects to test (var Effect). I use the code below: when I estimate the first model with "dd" it's work fine. Then I perform subset dd in

Re: [R] Coeficiente Gini

2014-01-30 Thread Achim Zeileis
On Thu, 30 Jan 2014, Gastón Acosta c. wrote: He le?do bastante sobre el programa R. No soy un programador. Estoy interesado en obtener un programa R que permita calcular los intervalos de confianza y el error de estimaci?n del coeficiente Gini, ya que el c?lculo de este y su gr?fico, as? como

[R] Coeficiente Gini

2014-01-30 Thread Gastón Acosta c .
He leído bastante sobre el programa R. No soy un programador. Estoy interesado en obtener un programa R que permita calcular los intervalos de confianza y el error de estimación del coeficiente Gini, ya que el cálculo de este y su gráfico, así como la curva de Lorenz, se logran por otros medios

[R] [R-pkgs] SPRINT 1.0.5 release (parallelised R functions for complex analysis and large data sets)

2014-01-30 Thread SPRINT Project
Dear All We have recently released SPRINT v1.0.5. SPRINT provides parallelised versions of existing R functions (e.g. statistical, machine learning, utility) that often exceed computational limits (speed or memory) with large or complex data sets. SPRINT 1.0.5 now: -> runs on Mac OSX mu

Re: [R] Parsing Complex Text in Single Cell

2014-01-30 Thread arun
Another way would be: library(qdap) library(stringr) x <- scan(what="character",)  x1 <- c(x,x) x1 <- paste(x1,collapse=" ")  x2 <- gsub('"',"",bracketXtract(x1,"curly"))  res2 <- as.data.frame(str_trim(do.call(rbind,genXtract(paste0(x2,","),":",","))),stringsAsFactors=FALSE) res2[,1:3] <-

Re: [R] Inconsistent results between first run of Rprof and next runs of Rprof

2014-01-30 Thread Nathan Uyttendaele
Hello again, turns out that I was able to nail down the problem to this : Prior to the profiling, I have several packages and a lot of functions to load, including the C code for the rankCR function in R. If I load first the libraries and next the C functions library(Hmisc) library(copula) libr

Re: [R] Controlling font size on code chunk outputs using Knitr

2014-01-30 Thread Jeff Newmiller
This sounds like a classic "you need to write a custom CSS file" problem... Which is off-topic here, so is homework for you. --- Jeff NewmillerThe . . Go Live... DCN:Basics: #

Re: [R] Controlling font size on code chunk outputs using Knitr

2014-01-30 Thread Jeff Johnson
Hi Yihui, The package I have installed is "knitr". To generate the HTML, I run Knit HTML from within R Studio version .98.490 (there's an icon to initiate it. Here's a simple example: showcode <- FALSE commentchar <- NA You can load this data as 'mydf'... dput(mydf) structure(list(PERSONPROFILE_

[R] print.matrix

2014-01-30 Thread Göran Broström
In the documentation of 'prmatrix' (base) I read, under Details: ‘prmatrix’ is an earlier form of ‘print.matrix’ but 'print' doesn't seem to have a 'matrix' method. And in the 'Examples' section: chm <- matrix(... chm # uses print.matrix() Is this a bug in the documentation? R-3.0.2 on ubun

[R] package mgcv - predict with bam: Error in X[ind, ] : subscript out of bounds

2014-01-30 Thread Katharina May
Dear R-Community, I`m trying to apply the mgcv package to fill gaps in sensor data from different sites (9 sites, 2 sensors per site) and do the filling on a site-wise level. Based on http://r.789695.n4.nabble.com/mgcv-gamm-predict-to-reflect-random-s-effects-td3622738.html my model looks like th

Re: [R] Parsing Complex Text in Single Cell

2014-01-30 Thread Rui Barradas
Hello, Maybe something like the following. x <- scan(what = "character", text = ' {"trial":1,"corr":1,"resp_dur":799,"stim":"â†*�*â†*�*â†*�*â†*�*â†*�*","cond ":"congruent"},{"trial":2,"corr":1,"resp_dur":0,"stim":"xx→xx","cond":" nogo"},{"trial":3,"corr":0,"resp_dur

[R] Regsubsets n

2014-01-30 Thread wild_manul
Hello I am trying to run a model where the number of observations is less than the number of predictors. When I've tried to run the regsubsets on a dummy dataset of random normally distributed numbers it gives me an error whenever I set n to be less than p. For instance, when number of observation