similar to: Re storing a UPDATES on a data.frame

Dear List-members, Hopefully someone will help through my confusion: In order to get the same coefficients as we get from the following ## require (MASS) summary ( lm(Gas ~ Insul/Temp - 1, data = whiteside) ) ...................... we need to do the following (if we use model.matrix to specify the model) ## summary ( lm(Gas ~ model.matrix(~ Insul/Temp - 1) - 1, data = whiteside) )

R folks, As an R novice, I struggle with the mystery of subsetting. Textbook and online examples of this seem quite straightforward yet I cannot get my mind around it. For practice, I'm using the code in MASS Ch. 6, "whiteside data" to analyze a different data set with similar variables and structure. Here is my data frame: ###subset one of three cases for the variable

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

Hi, How would one go about determining if the slope terms from an analysis of covariance model are different from eachother? Based on the example from MASS: library(MASS) # parallel slope model l.para <- lm(Temp ~ Gas + Insul, data=whiteside) # multiple slope model l.mult <- lm(Temp ~ Insul/Gas -1, data=whiteside) # compare nested models: anova(l.para, l.mult) Analysis of Variance

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

I think this site is always failing. It has an embedded video. http://www.sportsnet.ca/basketball/nba/free-association-podcast-much-toronto-raptors-fear-miami-heats-dwayne-wade-hassan-whiteside-joe-johnson/ I tried looking at the output of strace firefox > fferrors.txt 2>&1 grep ENOENT fferrors.txt hope that helps Note: I also have ffmpeg with nux's repo installed On Mon, May

R-helpers; I tried to get a superscripted 3 in the following xyplot example but failed: >data(whiteside) >xyplot(Gas ~ Temp | Insul, whiteside, panel = function(x, y, ...) { panel.xyplot(x, y, ...) panel.lmline(x, y, ...) }, xlab = "Average external temperature (deg. C)", ylab = paste(paste("Gas consumption (1000", expression(ft^3),")"), aspect

Page 140 of MASS uses the function xyplot. But I don't find it in R. Is there a package that I should load to use xyplot. Or there is a function with a different name in R that does the same thing as xyplot in S. xyplot(Gas ~ Temp | Insul, whiteside, panel = function(x, y, ...) { panel.xyplot(x, y, ...) panel.lmline(x, y, ...) }, xlab = "Average external temperature (deg. C)", ylab

We had some custom additions to our site.py file for a third party application. 'yum update' to 4.6 overwrote the file with no backup or warning. Not hard to repair, but it did have me worried there for a few minutes when the application failed to start. Is this a python issue, an upstream issue or a CentOS issue? Can something be done about it going forward? Jeff

I'm running a series of simple bivariate linear regressions on grouped data. I want to test the slopes to see if they are parallel. I normally use analysis of covariance to do so, looking at interaction between the covariate and the factor to make this determination. VR3 pp.149 - 154 has a very nice example of an ANOCOVA, ending with a discussion of this very operation. My question has

Oh, Boy. I am in trouble now.... I just installed Centos on a USB drive on my corp notebook. To not TOUCH my corp drive. I spent time with the drive partitioner to make sure that nothing was done to the internal hard drive... Well I missed something and I overwrote the encrypted bootloader on the hard drive. Now what? Can I rescue things? I am leaving for my flight in a couple hours to

Is there anyway to unset the extensions.conf definition of writeprotect=yes while in the CLI interface (or by other mechanism) to enable the dialplan save command? I accidentally overwrote my extensions.conf but still have a running copy of asterisk with the old dial plan running in memory. while it would not be difficult for me to rebuild what I lost - it would be easier if I could just save

Hello, Is there a way to sync the Unix password database with Samba's database? I know that Samba is able to update the Unix database if a user changes their password via smbpasswd, but what if a user changes their Unix password with passwd? I looked into libpam-smbpass, but it didn't seem to work; the documentation was rather sparse, so I just overwrote my /etc/pam.d/samba with the

Hi, I'm implementing storeRegToStackSlot() and, in order to store some specific registers (floating point regs and address regs) I've to copy them to more standard regs and copy these last ones to the slot. I tried to generate instructions that use physical registers, but by doing that I overwrote registers already assigned by the register allocator. Is it possible to use virtual registers

This situation arose in a computer practical this morning. A student overwrote "c" with a function. This messes up lots of things, but it is difficult to remove the function "c" from the global environment because it is used as a function in both "rm" and "remove". [The solution was to reassign "c" as a number, and then remove it.] Would it be

In "Modern Applied Statistics with S-Plus" 3rd ed., footnote on page 153 regarding a model lm(Gas~Insul/(Temp+I(Temp^2))-1,whiteside), I read "Notice that when the quadratic terms are present, first degree coefficients mean 'the slope of the curve at temperature zero', so a non-significant value does not mean that the linear term is not needed.

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(). I would go even further and say that the examples that are in lm() now should stay at the top. Because people may be used to

I have developed 3 small applications using rails (nothing production worthy, just some small stuff "prototype" style). And I often end up changing my database scheme in the middle of development. I also use "script/generate scaffold" (this is somehow relevant to my question). Whenever I change my database schema I end up having to regenerate the scaffold. In the past this

Thanks for the discussion. I do feel quite strongly that the variables should always be a part of a data frame. Then functions such as summary() and pairs() can operate on them all simultaneously.... regression is only one part of the analysis. And what if there are lots of variables? Have them all scattered about the workspace? One of them could be easily overwritten. The generic predict() will

Hello list, I'm new to this list, so please forgive my ignorance. I have searched R-help for some hints into what might be my problem, but I truly have no idea where to go from here. I have an object of approximately 15,000 rows and 2 columns. There are many duplicates in the first column, all with different corresponding values in the second column. For example (2 is duplicated):