similar to: How to include an externally defined NULL value in lm

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

Hi 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)

Hello I'm only user of R and have many little knowledge in programming but I permit to send you some whishes/suggestions for R. which.min like which(), which.min() should also include an argument arr.ind. Note that one can have it with which(a==min(a), arr.ind=TRUE) but if there is a reason to build a special function which.min, why not add also this nice argument? lag() If one wants to

Hi I have a few problems with tail/head when applied to multiple time series. I'm not sure as whether I did not understand the function or whether 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,

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

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:

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 not a

can you recommend a good manual for R that starts with a data set and gives demonstrations on what can be done using R? I downloadedR Langauage definition and An introduction to R but haven't found them overly useful. I'd really like to be able to follow some tutorials using a dataset or many datasets. The datasets I have available on R are Data sets in package 'datasets':

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,

Hi, I’d like to fit one GLM model for each possible combination of inputs (i.e. 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,

[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

Bug PR#9316 noted an inconsistency between the Cook's distance contours on plot.lm(x, which = 5) and the values given by cooks.distance(x) -- as shown in plot.lm(x, which = 4) -- for glms: http://bugs.r-project.org/cgi-bin/R/Analyses-fixed?id=9316;user=guest;selectid=9316 The suggested fix was to modify the contour levels by a dispersion factor, implemented as follows: dispersion <-

Dear all, since it could be of interest for the list, here I post a link to compiled kernel for XEN 3.1 on a DELL PowerEdge 1950 compiled with gcc (GCC) 4.1.2 20061115 (prerelease) (Debian 4.1.1-21) on a Debian4 distribution. This is provided AS IS AND WITHOUT ANY SUPPORT AND/OR ANY RESPONSIBILITY :-) The first link (7 MB) pack all the files you should have in /boot. The second link (12 MB)

Hello, First, thanks for this great piece of software. I'm new to dovecot and i've just converted my old courier configuration on pop/imap server side, everithing is working just fine on this side and speed gain expected was found (not enough for our busy servers :). Setup is done with 3 dovecot virtual servers load balanced with a Hardware load balancer (this part is a big point of the

Hello, I''ve this strange problem. When Xen kernel start and it try to start eth devices, network doesn''t work. It seems to be a problem with NICs firmware. Anyone had this problem? Do you know where I can download firmware upgrade? (on DELL''s site I''ve just searched...) Thanks, Carlo No virus found in this outgoing message. Checked by AVG Free Edition.

Hi We are about to publish a book, which contains figures made with R plots. An important detail that we did not take into account is that the book will not be printed in 4 colors (cmyk mode), but only 2 (black +"spotcolor"). The spotcolor we use is part of the big Pantone family. The problem is that both pdf() and postscript() offer either rgb or cmyk, but no spotcolors such as

hello I'm running debian testing, winth samba as P.D.C, with passdb backend = ldapsam ii samba 3.0.23c-1 a LanManager-like file and printer server fo ii slapd 2.3.25-1 OpenLDAP server (slapd) ii smbldap-tools 0.9.2-3 Scripts to manage Unix and Samba accounts st ii libnss-ldap 251-5 NSS module for using LDAP as a naming servic no libpam_ldap

Hello Willy, Here is the dump from one of my bitcode files: 0000000 42 43 c0 de 21 0c 00 00 25 05 00 00 0b 82 20 00 As expected, 0x42 (= B), 0x43 (= C), xc0 and 0xde are in correct order. In your case, the first byte is read as 37 (= 0x25). I wonder why? When you check the bytes yourself, you get expected results. When the same bytes are read from Stream object, you get a different result (maybe

Hello I looked on the archives but did not find answer for that... Running R CMD check for a package, i get an error: Error in parse(n = -1, file = file) : unexpected symbol at 3341: } But how can I find which file is guilty? What is this 3342 referring to? Finally the solution I found is to source() every file until I find the file where the problem is... Is there some better way to do it?

> As far as I know, ObjectFile is either COFF, ELF or MachO format, > statically compiled. But the pass manager for applying optimisation > passes works only on llvm-ir. I would not need to apply optimisations on the parts of the program in the object file. These are meant to be optimised beforehand. The only parts I want to run passes on are specific function that I know