similar to: generating contrast names

This is my first post to this list so I suppose a quick intro is in order. I've been using SPLUS 2000 and R1.6.2 for just a couple of days, and love S already. I'm reading MASS and also John Fox's book - both have been very useful. My background in stat software was mainly SPSS (which I've never much liked - thanks heavens I've found S!), and Perl is my tool of choice for

After a year my R programming style is still very "C like". I am still writing a lot of "for loops" and finding it difficult to recognize where, in place of loops, I could just do the same with one line of code, using "sapply", "lapply", or the like. On-line examples for such high level function do not help me. Even if, sooner or later, I am getting my R

Hi, I would like to use the sem package to perform a path analysis (no latent variables) with a mixture of 2 nominal exogenous, 1 continuous exogenous, and 4 continuous endogenous variables. I seek advice as to how to calculate the appropriate covariance matrix for use with the sem package. I have read through the polycor package, and am confused as to the use of "numeric" for

Is there any logical reason why glm prints out the labels of factor levels after variable names when baseline contrasts (contr.treatment) are used but the codes for the levels when mean contrasts (contr.sum) are used? Jim -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I have been trying to learn the various "apply" functions but am still learning their appropriate use. I appreciate any help the R community can offer me. Sorry for the length of this post. Background: I have data on my hard drive organized in the following manner: The data pertains to many different "samples" of data. (e.g. sample 001, sample, 002, sample 003, etc.) Each

I would like to plot a survival curves of a group with different categories after running a Cox model with frailty() random effects terms. I just could display a survival plot of the covariable?s mean. Here an example: library(survival) fit<-coxph(Surv(time, status) ~ sex+ frailty(litter, dist='gamma', method='em'), rats) summary(fit ) suf<-survfit(fit) plot(suf,

Dear R list members, I'm encountering a problem with expand.model.frame(): Suppose that I define the following simple function (meant just to illustrate the problem): > fun <- function(model){ + expand.model.frame(model, all.vars(formula(model))) + } > and I have the following model, created with an explicit data argument: > mod Call:

Hi, I've got a Powerware Prestige 9 1kVA UPS that seems to not be working with nut. It's not in the current "supported hardware" list, but the old bcmxcp page lists it as being a UPS that speaks that protocol. I have the following in my ups.conf [powerware] driver = bcmxcp port = /dev/ttyS0 desc = "Powerware Prestige 9 EXT" baud_rate = 2400 And when I start the

I am getting an error message from scatterplot: > library(car) > scatterplot(Prestige$income~Prestige$type) Error in Summary.factor(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, : range not meaningful for factors In addition: Warning message: In Ops.factor(x[floor(d)], x[ceiling(d)]) : + not meaningful for factors > The command does output the kind of graph that I want (boxplots).

Hi, Does anyone have suggestions hooking Zyxel Prestige 2002 to Asterisk? I have tested Zyxel Prestige with both supported codecs. Call with G.711 sounds very choppy and cracking. Almost can't understand a word. Today I installed G.729 support into Asterisk but unbearable voice quality remains. It's a little bit better though. I have tested that Zyxel ATA with some commercial SIP

Does anybody have any experience with the ZyXEL Prestige 2000W? I am having problems with the line tear down when I call another extension. If nobody picks up at the other end when I hangup the 2000W, the other extension continues to ring. Is there any way to hangup a SIP call if there is no more traffic? Asterisk seems to think that there is still a connection open. This is pretty annoying

Dear List, I don't get scatter3d to color the sheres according to the '|' argument. library(car) scatter3d(prestige ~ income + education|type, data=Prestige) The spheres on my screen are all colored the same and they are not conditional on Prestige$type. On the other hand: Fit3d and Ellipse3d are colored according to the group argument. rgl_0.92.879 car_2.0-12 R version 2.15.0

Hello Rui, Thanks for your helpful suggestions. Just for illustration, let's use the well known Duncan dataset of prestige vs education + income that is contained in the "car" package. Suppose I wish to use boot function to bootstrap a linear regression of prestige ~ education + income and use the following script: duncan.function <- function(data, indices) {data =

Dear R-help list members, I'm encountering problems using expand.model.frame. To illustrate, consider the function > test <- function(model){ + expand.model.frame(model, "income") + } > The data frame Prestige (from the car library) has several variables, including prestige, income and education. I've attached this data frame >

Hi all, I've more-or-less got NUT 2.0.4 working with my Powerware Prestige 1000VA UPS (serial), but unfortunately NUT is flagging LOW_BATTERY continuously, which is a tad embarrassing when power fails ;-) The numbers from upsc blah at localhost seem reasonable: ambient.temperature: 24.5 battery.charge: 99 battery.runtime: 3999 battery.voltage: 69.2 driver.name: bcmxcp

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, Anyone have any experience with either of these, I 'd appreciate some feedback? Plus it seems pretty easy to steal a connection with this. Zyxel Prestige 2000W WiSIP thanks, - -- Steve "They that would give up essential liberty for temporary safety deserve neither liberty nor safety." Benjamin

If you see a wj0011 version of firmware for Zyxel Prestige 2000W floating around (I found it in a German forum), KEEP AWAY. It completely trashed the wireless networking in my phone. -- ========================================== Rod Bacon Empowered Communications Ground Floor, 102 York St. South Melbourne Victoria, Australia. 3205 Phone: +613 99401600 Fax: +613 99401650 FWD: 512237

Hello, does NUT support Powerware Prestige EXT 1000VA / PowerWorks A30 models in smart mode? I don't see them mentioned in the compatibility list (only Powerware 3115 and 5119 are there, listed as supported only by the "genericups" driver). Is the protocol known to be compatible with any other, supported UPS types? Thanks, Marek

Dear Trevor, I've noticed a discrepancy in the degrees of freedom reported by gam() from the gam package in R vs. gam() in S-PLUS. The nonparametric df differ by 1; otherwise (except for things that depend upon the df), the output is the same: --------- snip ------------ *** From R (gam version 0.93): > mod.gam <- gam(prestige ~ lo(income, span=.6), data=Prestige) >

Hello, Much clearer now, thanks. It's a matter of changing the function boot calls to return the predicted values at the point of interess, education = 50, income = 75. I have changed the way the function uses the indices a bit, the result is the same, it's just the way I usually do it. pred.duncan.function <- function(data, indices) { mod <- lm(prestige ~ education +