similar to: R-beta: Assigning column names in a data frame

I am trying to show some techniques to my graduate regression class. The textbook mentioned using bootstrap samples of regression coefficients for assessing variability. I decided to show them reasonably effective ways of doing the resampling. The following is a function I wrote to create bootstrap samples of coefficients from a fitted linear regression model. bsCoefSample <- ##

Greetings, In my opinion the documentation or behaviour of qr.solve, qr.coef, qr.resid, and qr.fitted is not easily comprehensible and unfortunate. We all know that a linear system Ax=b can have 0, one or infinitely many solutions. To treat all these cases uniformly we can rephrase the problem as x = argmin_u||Au-b||,

When regressing by month, how do I get the coefficients out into a new data set? I'm looking for [ month, a, b, c ] from the Pastor-Stambaugh model I'm using which is: r[i+1] = a + b * r[i] + c * v[i] + e the model I'm using wants to create a new dataseries based on the coefficient in each month. I'm doing a simple linear regression on DataSet, and >

R developers, I have a first attempt to make an aov function. Eventually I want to build in Error() structure, but first I am trying to get this presentable for balanced data with only a single stratum, just using residual error. I am following R. M. Heiberger's Computation for the Analysis of Designed Experiments, Wiley (1989) I a using a wrapper (aov.bal) to call the

> If I were Mathsoft I would be less than pleased > at the development of R and would try to stop it if I could. I would not assume that Mathsoft must have such a negative view towards R. The reasons are roughly summarized by the analogy that a smaller part of a big pie is often better than a bigger part of a small pie. Mathsoft's limited success in addressing the student market has

There may be a bug in the anova.glm function. deathstar[32] R R : Copyright 2004, The R Foundation for Statistical Computing Version 2.0.1 (2004-11-15), ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project

This is an Automatic Reply ********************************** Thank you for e-mailing Mathsoft Technical Support. We offer support for Mathcad 2001i or higher, Axum 6 or higher and the StudyWorks product line. Cumulative Update 1 for Mathcad 11 Now Available! http://www.mathcad.com/download/ Mathcad 2001i Service Release 2 was posted June 2002. New Support Policy: With the release of Mathcad

Hi list, I tried to use nls to do some nonlinear least square fitting on my data with 340 observations and 10 variables, but as I called nls() function, I got this error message: Error in qr.qty(QR, resid) : 'qr' and 'y' must have the same number of rows Then I traced back a little bit into nls() function, the error seemed to happen when calling nlsiter() internal function, but

Martyn Plummer writes: > How important is it to avoid being sued, or less facetiously, what is > the legal status of R? If I were Mathsoft I would be less than pleased > at the development of R and would try to stop it if I could. I have > been wondering for some time if this is possible. But I am not a lawyer > and the issue seems very unclear to me. To me too (Robert spent a

Hi, I installed R-3.1.1 on AIX-5.3 and my installation hanged up while installing Matrix Package. so i killed the gmake ; gmake install the software nad tired installing Matrix package manually - [ncmr0202][/gpfs1/home/shivali/gang/R-3.1.1/bin/package]> /gpfs1/home/shivali/gang/R-3.1.1/bin/R CMD INSTALL Matrix the package compiled successfully but while loading Matrix package the

Dear R-experts! I have y(x) data, dim(y) == dim(x) == c(2000, 2) I'd like to fit them with nls: fit.result <- nls ( y ~ f(x, p1, p2, p3), start = list(p1 = ... , p2 = .. , p3 = ..) ) Actually I want to fit y[,1] ~ x[,1] and y[,2] ~ x[,2] *simulaneously*, with the same parameters set {p1, p2, p3}. I tried to feed R tha above formula, R errors with:

I am trying to figure out why 'biglm' can handle large data set... According to the R document - "biglm creates a linear model object that uses only p^2 memory for p variables. It can be updated with more data using update. This allows linear regression on data sets larger than memory." After reading the source code below? I still could not figure out how 'update'

I am setting up a test domain to learn bind, and how it,s working. The master domain is bbb.ttt.domene 10.11.12.13 the child domain ftp.bbb.ttt.domene 10.11.55.77 The thing I don't understand is, why don't the query report additional info (ip address)? Her is two queries(the first the one that troubles me), cut of named.conf, test of zone file, the zone file and cut from log file.

Full_Name: Christian Brechbuehler Version: 2.4.1 Patched (2007-03-25 r40917) OS: Linux 2.6.15-27-adm64-xeon; Ubuntu 6.06.1 LTS Submission from: (NULL) (24.61.47.236) Splus and R have different ideas about what qr.coef(qr()) should return, which is fine... but I believe that R has a bug in that it is not internally consistent, and another separate bug in the documentation. In particular, on

Hi, I have a problem with fitting a model: I made a dataframe with this data: a <- 1:3 b <- 1:3 c <- c(3, 2, 3, 2, 1, 2, 3, 2, 3) df <- expand.grid(a,b) df$result <- c names(df) <- c("A","B", "result") Although I can make a graph of the data: require(lattice) wireframe(result~A*B, data=df) I can't get a model to

I suggest adding a 'pivot' argument to qr.R, to obtain columns in the same order as the original x, so that a <- qr(x) qr.Q(a) %*% qr.R(a, pivot=TRUE) returns x. -------------------------------------------------- # File src/library/base/R/qr.R qr.R <- function(qr, complete = FALSE, pivot = FALSE) { # Args: # qr: a QR decomposition, produced by qr() # complete:

Hi all, Regular and avid readers of this column will know that Don Driscoll and I have recently posted two messages requesting assistance concerning an apparent failure of "se.contrast" to produce an se for a contrast. So far, an ominous silence rings in our ears, but read on Gentle Reader, and see if even the machinations of "debug" doesn't stimulate you to respond with a

I think there are two bugs in aov() that shows up when the right hand side of `formula' contains both `-1' and an Error() term, e.g., aov(y ~ a + b - 1 + Error(c), ...). Without `-1' or `Error()' there is no problem. I've included and example, and the source of aov() with suggested fixes below. The first bug (labeled BUG 1 below) creates an extra, empty stratum inside

Here is a modified version of qraux.Rd, an edited version of R-2.14.0/src/library/base/man/qraux.Rd This gives some details and an example for the case of pivoting. In this case, it is not true that X = QR; rather X[, pivot] = QR. It may save some other people bugs and time to have this information. Tim Hesterberg -------------------------------------------------- % File

Dear all, I am running a code using bootstraps for estimating standard errors but the mac crashes. When I use small number of bootstraps (100) it works fine but if I increase that number it crashes. Thanks in advance dimitris The code, the error and my mac characteristics are the following: ##############code##################### qr.1<- rq(y~factor(year)+factor(state)+x1+I(x^2)+I(x^3),