similar to: Baffled by drop1

My apologies to the list for sending this without adequate research. I have found my answer; please ignore! Thanks. I've been experimenting with drop1 for my biostatistics class, to obtain the so-called Type III sums of squares. I am fully aware of the deficiencies of this method, however I feel that the students should be familiar with it. What I find baffling is that when applied to a fully

I have a very large dataframe and wish to extract a subset of rows. I have a two column matrix listing the starting and ending indices of one subset on each row. My idea is to create a vector of indices that could be applied to the dataframe and I have a solution using a for loop (below). But surely there is some more elegant way to do this! I looked thorough the archives without

In S-Plus, there is a function ms() [minimum sums] that estimates parameters in a multivariate model and returns, along with the estimates, a loglikelihood score. It's very handy, but I have not found an R version yet. Does anyone know if such a thing exists? Thanks. Cheers, George ================================================================== George W. Gilchrist

I need to extract the degrees of freedom and log likelihoods from a series of mixed models computed using lmer/lme4. If I ask for logLik (lmer.object), I get something like > logLik(lmer.object) 'log Lik.' -177.1000 (df=10) Can I easily get that df from there (or elsewhere) into an object? Thank you for any ideas. George

Does anyone have an idea of how one might plot male and female symbols on a graph using R? Thanks! .................................................................. George W. Gilchrist Email #1: gwgilc at wm.edu Department of Biology, Box 8795 Email #2: kitesci at cox.net College of William & Mary Phone: (757) 221-7751 Williamsburg, VA

I am trying to do a series of xyplots plots, using a "for" loop to substitute the appropriate variables for each plot. The basic command works fine by itself and produces a nice plot: > i<-3 > trellis.device(theme="col.whitebg") > xyplot(as.formula(paste(tmp00[2*i], "~ ", tmp00[(2*i)-1], +

For several years, I have been using Splus to analyze an ongoing series of datasets that have a bimodal distribution. I have used the following functions, in particular the ms() function, to estimate the parameters: two means, two standard deviations, and one proportion. Here is the code I've been using in S: btmp.bi <- function(vec, p, m1, m2, sd1, sd2) {

In S-Plus, I can obtain polynomial contrasts for an ordered factor with contr.poly(). The function also exists in R, however is limited to factors where the levels are equally spaced. In S-Plus, one can obtain the contrasts for a set of numeric values representing unequally spaced ordered factors. Has anyone implemented this in R? I see that the S-Plus function calls another function (poly.raw())

A recent question in r-help made me realize that I should add a drop1 method for coxph and survreg. The default does not handle strata() or cluster() properly. However, for coxph the right options for the "test" argument would be likelihood-ratio, score, and Wald; not chisq and F. All of them reference a chi-square distribution. My thought is use these arguments, and add an

On Wednesday 07 May 2003 22:21, Michael Collette wrote: > Scenario: > FreeBSD box running IPFW acting as a gateway to private network. The > private network is made up of entirely routeable IP addresses. External > users running Win2k and XP on DSL connections with dynamic IPs. > > Goal: > To have the FreeBSD gateway securely authenticate and encrypt the traffic > between

1.6.2.13, sip.conf: [155] type=friend context=longdistance callerid="Admin" <155> secret=test host=dynamic dtmfmode=rfc2833 allow=all defaultuser=155-trust ............ On aastra: Basic SIP Authentication Settings Screen Name Phone Number 155 Caller ID 155 Authentication Name 155-trust Password test But: WARNING[1737]:

Hello. I am converting my Java app to Rails and I''ve hit something that has me a bit baffled. I''ve been fiddling with this for over an hour and trying many different combinations so at this point I''ve lost perspective. It''s probably something simple but I can''t see it anymore. I have a component that looks like this: <% for advertiser_category in

I'd like to fit and explore a collection of hierarchical loglinear models that might range from the independence model, ~ 1 + 2 + 3 + 4 to the saturated model, ~ 1 * 2 * 3 * 4 I can use add1 starting with a baseline model or drop1 starting with the saturated model, but I can't see how to get the model formulas or terms in each model as a *list* that I can work with further. Consider

the data is : >table.8.3<-data.frame(expand.grid( marijuana=factor(c("Yes","No"),levels=c("No","Yes")), cigarette=factor(c("Yes","No"),levels=c("No","Yes")), alcohol=factor(c("Yes","No"),levels=c("No","Yes"))), count=c(911,538,44,456,3,43,2,279))

Dear List, recently tried to reproduce the results of some custom model selection function after updating R, which unfortunately failed. However, I ultimately found the issue to be that testing with pchisq() in drop1() seems to have changed. In the below example, earlier versions (e.g. R 2.4.1) produce a missing P-value for the variable x, while newer versions (e.g. R 2.7.1) produce 0 (2.2e-16).

This is probably simple, but I just can't see it... I want to calculate the R^2s for a series of linear models where each term is dropped in turn. I can get the RSS from drop1(), and the r.squared from summary() for a given model, but don't know how to use the result of drop1() to get the r.squared for each model with one term dropped. Working example: library(vcd) # for

If you can help me in any way, it would be much appreciated. e3po at e3po-linux:/media/UDF Volume$ wine Startup.exe fixme:powrprof:DllMain (0x7e8c0000, 1, 0x1) not fully implemented fixme:ole:CoInitializeSecurity ((nil),0,(nil),(nil),1,3,(nil),0,(nil)) - stub! fixme:msxml:bsc_QueryInterface interface {6d5140c1-7436-11ce-8034-00aa006009fa} not implemented fixme:win:WIN_CreateWindowEx Parent is

Yes, it's homework . . . delete now if desired . . . but I think it is an interesting problem. Running R 2.11.1, LaTeX on WinXP, via Sweave. A drop1() object from a glm() produces, as part of its output, a string that looks like this: <none> The trouble I run into is that running latex() on a drop1() object from glm() produces a string that looks like this in the generated .tex

Hi, If I used drop1 in a weighted lm fit, it seems to ignore the weights in the AIC calculation of the dropped terms, see the example below. Can this be right? Yan -------------------- library(car) > unweighted.model <- lm(trSex ~ (river+length +depth)^2- length:depth, dno2) > Anova(unweighted.model) Anova Table (Type II tests) Response: trSex Sum Sq Df F value

Our backup mail server (which I have just re-configured) tries to contact the primary mail server, and fails. My log shows repeatedly "connection refused": May 15 22:21:41 mx2 sm-mta-rx[8674]: q4FIhPij007483: makeconnection (mail.greenspot.fi. [83.143.217.182]) failed: Connection refused by mail.greenspot.fi. May 15 22:21:41 mx2 sm-mta-rx[8674]: q4FIhPij007483: to=<myuser at