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 therefore the most di

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 + B3.tensionH + B4. woolB

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, tensio

[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 ## libr

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 runn

[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 ## libr

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) > library(TSA) > security<-"NFLX

[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 lo

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 P

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! C:/Rtools/mingw_6

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: [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 <- with