Re: [R] Interpreting summary.lm for a 2 factor anova

2016-12-03 Thread Ashim Kapoor
Dear Sir, Alright. Best Regards, Ashim On Sun, Dec 4, 2016 at 10:59 AM, Richard M. Heiberger wrote: > As Petr Pikal mentioned, the difficulty in interpretation is entirely due > to the set of contrasts you chose.The default treatment contrasts are > not orthogonal and are

Re: [R] Interpreting summary.lm for a 2 factor anova

2016-12-03 Thread Ashim Kapoor
On Sun, Dec 4, 2016 at 10:03 AM, Ashim Kapoor wrote: > Dear Sir, > > Many thanks for the explanation. Prior to your email (with some help from > a friend of mine) I was able to figure this one out. If we look at the > model : - > > y = intercept + B1.woolB + B2. tensionM +

Re: [R] Interpreting summary.lm for a 2 factor anova

2016-12-03 Thread Ashim Kapoor
Dear Sir, Many thanks for the explanation. Prior to your email (with some help from a friend of mine) I was able to figure this one out. If we look at the model : - y = intercept + B1.woolB + B2. tensionM + B3.tensionH + B4. woolB.TensionM + B5.woolB.TensionH + error Here woolB, tensionM,

[R] tcltk: use of .Tcl.callback

2016-12-03 Thread Cleber N.Borges
Dear, How to properly use the function: .Tcl.callback( ) to trigger a single R function with different values? Below is my stupid solution ... The expected behavior is this but I would like to know how to do it right. Thanks in advanced cleber ##

Re: [R] error serialize (foreach)

2016-12-03 Thread Doran, Harold
As a follow up to this, I have been able to generate a toy example of reproducible code that generates the same problem. Below is just a sample to represent the issue, but my data and subsequent functions acting on the data are much more involved. I no longer have the error, but, the loop

[R] tcltk: use of .Tcl.callback

2016-12-03 Thread Cleber N.Borges
Dear, How to properly use the function: .Tcl.callback( ) to trigger a single R function with different values? Below is my stupid solution ... The expected behavior is this but I would like to know how to do it right. Thanks in advanced cleber ##

Re: [R] Bootstrap using ARIMA model

2016-12-03 Thread David Winsemius
> On Dec 1, 2016, at 1:58 PM, Ashwini Patil wrote: > > Hi David, > > here is my code including what i did for the tsboot: > rm(list = ls()) > library(boot) > library(tseries) > library(TTR) > library(quantmod) > library(scales) > library(forecast) > library(zoo) >

[R] error serialize (foreach)

2016-12-03 Thread Doran, Harold
I have a portion of a foreach loop that I cannot run as parallel but works fine when serialized. Below is a representation of the problem as in this instance I cannot provide reproducible data to generate the same error, the actual data I am working with are confidential. Within each foreach

Re: [R] How to setup a multiplicative dummy function in R

2016-12-03 Thread Jeff Newmiller
I think you need an offset term, or maybe I just don't understand your question. A sample data set, particularly if you can show us how your equation could be used to generate the sample data, would be helpful. -- Sent from my phone. Please excuse my brevity. On December 1, 2016 7:22:37 AM

Re: [R] rstan error: C:/Rtools/mingw_64/bin/g++: not found

2016-12-03 Thread Uwe Ligges
On 02.12.2016 12:23, S Ellison wrote: Apologies for posting a possibly package-specific question, but I'm not sure whether this is an R or rstan ussue. Running rstan under R 3.1.1 in windows 10 I get the well-known error "Compilation ERROR, function(s)/method(s) not created!

Re: [R] data

2016-12-03 Thread jim holtman
This should be reasonably efficient with 'dplyr': > library(dplyr) > input <- read.csv(text = "state,city,x + 1,12,100 + 1,12,100 + 1,12,200 + 1,13,200 + 1,13,100 + 1,13,100 + 1,14,200 + 2,21,200 + 2,21,200 + 2,21,100 + 2,23,100 + 2,23,200 + 2,34,200 + 2,34,100 + 2,35,100") > > result <- input

[R] data

2016-12-03 Thread Val
Hi all, I am trying to read and summarize a big data frame( >10M records) Here is the sample of my data state,city,x 1,12,100 1,12,100 1,12,200 1,13,200 1,13,100 1,13,100 1,14,200 2,21,200 2,21,200 2,21,100 2,23,100 2,23,200 2,34,200 2,34,100 2,35,100 I want get the total count by state, and

Re: [ESS] how to edit openbugs model file?

2016-12-03 Thread Sparapani, Rodney
Hi Paul: Use BUGS mode which is part of ESS… (require ‘ess-bugs-d) Then name your file model.bug Rodney __ ESS-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/ess-help

Re: [R] Interpreting summary.lm for a 2 factor anova

2016-12-03 Thread Fox, John
Dear Ashim, Sorry to chime in late, and my apologies if someone has already pointed this out, but here's the relationship between the cell means and the model coefficients, using the row-basis of the model matrix: -- snip > means <-

Re: [ESS] how to edit openbugs model file?

2016-12-03 Thread Martin Maechler
Also, M-x R-mode but that would need to be done again every time... If you don't want to rename to model.R Add a line ## -*- R -*- at the top (or close enough to the very top; this tells emacs to use R-mode ; there are several (more flexible, more to write) versions of such so-called "file