search for: freeny

...i Say I want to add manually an intercept in the function lm. Even if almost all results will be identical, few stats are different as DF counting will be different as intercept will not be included in "automatic" case, while it will be in "manual" case. See: ###usual lm on freeny fr<-lm(freeny.y~freeny.x) ###manual lm on freeny man<-cbind(1,freeny.x) colnames(man)<-c("const",colnames(freeny.x)) fr_man<-lm(freeny.y~man-1) ###coef are the same cbind(coef(fr), coef(fr_man)) ###but summary output is different (but should be the same!). #Difference comes...

Full thread here: https://github.com/tidyverse/broom/issues/287 Reproducible example: is.object(freeny$y) # [1] TRUE attr(freeny$y, 'class') # [1] "ts" class(freeny$y) # [1] "ts" # ts attribute wiped by model.frame class(model.frame(y ~ ., data = freeny)$y) # [1] "numeric" attr(model.frame(y ~ ., data = freeny)$y, 'class') # NULL # but still: is.object...

Hi My goal is to do a (multiple) regression, just knowing that my Y variables will be the say k first variables of a matrix/data frame. I thought I should do it with eval(parse)) but encounter a strange problem. See: lm(y~.-y, data=freeny) #that's what I want to do in the one equation case #Problem is I don't know name of the variable... only that it is the first one... #so idea is to just take first name a<-colnames(freeny) #and then use eval(parse(text=a[1])) #it works if I replace y on either the left or right side:...

Hello! I would love to be able to include an external variable to a lm call, I mean something: if(TRUE) a<-freeny.x[,4] else a<-NULL lm(freeny.y~freeny.x[,-4] +a) but it does not work with a<-NULL, whereas lm(freeny.y~freeny.x[,-4] +NULL) I don't understand why and did not find an answer in the manuals... do you see it? Any idea? Thanks!!

Dear all, I have discovered an object named 'freeny.x' which seems in-built. Can somebody tell me why it is there and what purpose it is serving? 'ls()' syntax does not seem to catch this object in the environment. Thanks,

This one showed up while looking at one of Ripley's other reports: > data(freeny) > model.frame(y~1,data=freeny,subset=1:10) y 1962.25 8.79236 1962.5 8.79137 1962.75 8.81486 1963 8.81301 1963.25 8.90751 1963.5 8.93673 1963.75 8.96161 1964 8.96044 1964.25 9.00868 1964.5 9.03049 > model.frame(y~1,data=freeny,subset=1:10)$y Warning: Replacement length...

[Copied to R-devel for information] This applies to all versions of R I have: 0.62.2, 0.62.3, 0.63. Great care seems needed with glms with offsets, as many things seem wrong. Consider the following: > data(freeny) > freeny.glm <- glm(y ~ offset(lag.quarterly.revenue) + price.index + income.level + market.potential, data=freeny, subset=1:30) > predict(freeny.glm) Qtr1 Qtr2 Qtr3 Qtr4 1962: NA 0.01040457 0.01073223 0.01233351 1963: 0.01211730 0.02744293 0.0325...

Hi, A recent thread provided a (working) construct for lm: lm(as.matrix(freeny[ix]) ~., freeny[-ix]) Can someone explain what is meant by the formula in that expression, that is, what does "mymatrix~." do? I couldn't find any such example in the lm() or formula() help pages. thanks Carl

...Conversion Rates of Euro Currencies euro.cross (euro) Conversion Rates of Euro Currencies eurodist Distances Between European Cities faithful Old Faithful Geyser Data fdeaths (UKLungDeaths) Monthly Deaths from Lung Diseases in the UK freeny Freeny's Revenue Data freeny.x (freeny) Freeny's Revenue Data freeny.y (freeny) Freeny's Revenue Data infert Infertility after Spontaneous and Induced Abortion iris Edgar Anderson's Iris Data ir...

...exhaustive search). The package leaps can help me to generate all possible variable subsets, but I’ll appreciate your guidance as of how to generate one model for each of those possible subsets. I’m new in R! Thanks in advance for your help. library(leaps) subsets <- regsubsets( y ~ . , data=freeny, nbest=100, nvmax=5, method="exhaustive", really.big=T) subset.models <- summary(subsets)$which now fit one lm for each possible combination of variables. [[alternative HTML version deleted]]

...hether it correspond to an unexpected behavior. When head(a,n) is applied on data.frame or matrix, it returns a data-frame or matrix with first n obs of *each* variable. When applied to a mts object, it returns first n obs of *first* variable only, not of all... The same for tail(). See: head(freeny) ###mts object head(EuStockMarkets) #is equivalent to: head(EuStockMarkets[,1]) I guess it comes from absence of a head method for mts. Does it seem reasonable to have also a head.mts or did I misunderstand something? Thanks

...ong the contour, p is the number of parameters and hh is a vector of hat values. It is clear that in the case of a normal linear model, the cl.h values will depend on whether the model was fitted by lm or glm. The following example illustrates this: par(mfrow = c(2, 2)) plot(lm(y ~ ., data = freeny), 4:5) plot(glm(y ~ ., data = freeny), 4:5) For the lm fit we have crit = (cl.h^2 * hh)/(p * (1 - hh)) where cl.h is on the scale of the standardized residuals. This is the Cook's distance as defined e.g. in the Cook & Weisberg reference on ?cooks.distance. For the glm fit we have cr...

...rary/base/data/esoph.R action: create D: file: /usr/local/lib/R/library/base/data/euro.R action: create D: file: /usr/local/lib/R/library/base/data/eurodist.R action: create D: file: /usr/local/lib/R/library/base/data/faithful.R action: create D: file: /usr/local/lib/R/library/base/data/freeny.R action: create D: file: /usr/local/lib/R/library/base/data/infert.R action: create D: file: /usr/local/lib/R/library/base/data/iris.R action: create D: file: /usr/local/lib/R/library/base/data/iris3.R action: create D: file: /usr/local/lib/R/library/base/data/islands.R action: create...

...text html latex example formatc text html latex example formula text html latex example fourfoldplot text html latex example frame text html latex freeny text html latex example ftable text html latex example ftable.formula text html latex example function text html latex example gc...