similar to: ols Error : missing value where TRUE/FALSE needed

Dear R help, I'm having some trouble with model formulas for the ols function in the rms package. I want to have two variables represented as restricted cubic splines, and also include an interaction as a product of linear terms, but I get an error message. library(rms) d <- data.frame(x1 = rnorm(50), x2 = rnorm(50), y = rnorm(50)) ols(y ~ rcs(x1,3) + rcs(x2,3) + x1*x2, data=d) Error in

I posted this yesterday to r-help and Ben Bolker suggested reposting it here... Dickison, Daniel <ddickison <at> carnegielearning.com> writes: > > The documentation for agrep says it uses the Levenshtein edit distance, > but it seems to get this wrong in certain cases when there is a > combination of deletions and substitutions. For example: > > >

Hello, I have a couple of questions about the ols function in Frank Harrell's rms package. Is there any way to specify variables by their column number in the data frame rather than by the variable name? For example, library(rms) x1 <- rnorm(100, 0, 1) x2 <- rnorm(100, 0, 1) x3 <- rnorm(100, 0, 1) y <- x2 + x3 + rnorm(100, 0, 5) d <- data.frame(x1, x2, x3, y) rm(x1, x2, x3,

Hi, I wonder if anyone can answer the following or point me in the direction of how to obtain answers to the questions. Below is Output from R and further down are the questions raised and explanation of the study. Output from R: glm(formula = CB95TO00 ~ URB + INC, family = poisson) Deviance Residuals: Min 1Q Median 3Q Max -1.2272 -1.1290 0.2709 0.4272 2.1376

Dear All, I'm trying to use the ols function in the Design library (version 2.1.1) of R to estimate parameters of a linear model, and then use the contrast function in the same library to test various contrasts. As a simple example, suppose I have three factors: feature (3 levels), group (2 levels), and patient (3 levels). Patient is coded as a non-unique identifier and is

Dear R users, I tried to fit a cox proportional hazard model to get estimation of stratified survival probability. my R code is as follows: cph(Surv(time.sur, status.sur)~ strat(colon[,13])+colon[,18] +colon[,20]+colon[,9], surv=TRUE) Error in if (!length(fname) || !any(fname == zname)) { : missing value where TRUE/FALSE needed Here colon[,13] is the one that I want to stratify and the

hi all - i have a function that needs to call glm() with a weights argument that includes a variable whose name comes from the caller. so instead of: fit <- glm(formula, poisson(), data, weights = 1-z, ...), i do something like this: fit <- do.call("glm", list(formula=formula, family=poisson(), data=data, weights = call("-", 1, as.name(zname)), ...))

Hello, Is is possible to do "shorewall reject 1.1.1.1 tcp 25" ? So I can dynamically blacklist offensive smtp senders, but only have shorewall reject certain types (smtp) traffic from them? Thanks, Alex Martin http://www.rettc.com

Hi, A error message arose while I was trying to fit a ordinal model with lrm() I am using R 2.8 with Design package. Here is a small set of mydata: RC RS Sex CovA CovB CovC CovD CovE 2 1 0 1 1 0 -0.005575280 2 2 1 0 1 0 1 -0.001959580 2 3 0 0 0 1 0 -0.004725880 2 0 0 0 1 0 0 -0.005504850 2 2 1 1 0 0 0 -0.003880170 1 2 1 0 0 1 0 -0.006074230 2 2 1 0 0 1 1 -0.003963920 2 2 1 0 0 1 0

Hello, I have problem running WinBUGS from R. The following example works in WinBUGS but it does not work in R through package R2WinBUGS. Does anyone know what the problem is? x <- c(0.2, 1.1, 1, 2.2, 2.5, 2.9, 2.9, 3.6, 3.8, 0.6, 1, 2, 2.4, 2.6, 2.8, 3.2, 3.9, 3.5) y <- c(0.5, 1.3, 0.1, 0.7, -0.4, 0.5, -0.9, -0.3, -0.3, 0.6, 0.4, 0.9, -0.1, -0.4, -0.5, -0.2, 0.3, -1.5) eco <- c(1, 3,

Hello all; I have been writing this LLVM pass that extracts loops in a function and then duplicates the function so the loop is executed twice. The pass executes fine but the produced bit code does not behave correctly as it should (both the original and duplicated function in which the loop is extracted is called, but from the output it seems only one of them is executed.) Here is the source:

While it was a good move from Rails part to escape ERB <%= %> tags by default, it doesn''t seem to happen to Sprockets as well. The strange bit is that according to Sprockets documentation, it would be just a matter of naming your template as .jst.eco to enable Eco: https://github.com/sstephenson/sprockets#javascript-templating-with-ejs-and-eco Then, extracted from Eco

From: saint-filth@hotmail.com To: saint-filth@hotmail.com Subject: RE: Date: Wed, 9 Jun 2010 09:53:20 +0000 OK! sorry thats my fault, here the translations of the stata commands 1st step is to get the mean values of the variables, well that doesnt need explanation i guess, 2nd step is to estimate the model on panel data estimation method which is:

Hello R users: I would like to produce a world map with countries colored according to whether they fall into one of 7 eco-climatic zones. For simplicity, each country is allocated to exactly 1 eco-climatic zone. For this purpose I have looked at the map and mapdata packages, which contain world maps composed of polygons (1 for each country, it seems). Each of the polygons can be referred to by a

Dear testers, 5.0.0-rc3 was just tagged. This is a release candidate in the real sense: if nothing bad comes up in testing, this is what the release is going to look like. Please build, test and upload binaries to the sftp (use the /data/testers-uploads/ directory) and let me know what issues remain. I know we're a little bit behind schedule, but hopefully we can get to 'final'

All, I've a question on A*k hardware. I'm running 1.2.18 on a Dell DC051 (Intel(R) Celeron(R) CPU 2.80GHz) with 512mb RAM. I'm supporting 60 users (Cisco 7940s each + Xlite PCs). Call loads are low, max of about 10 simultaneous SIP/IAX calls. CPU for A*k rarely goes above 2% as I can tell. Its IP only, no E1/T1 cards present. However, I get complaints of bad voice quality,

On Wed, 10 Apr 2019, Bartosz wrote: > BTW,does it have passive cooling? In the sense that there is no fan, The Eaton Ellipse ECO 1600 has passive cooling. It gets warm but not hot. > Are you happy with your Eaton Elipse ECO 1600? Does it work well with ondelay? "ECO" means that it does not have "buck and boost" or a fancy control panel. I replace my batteries

On 02.10.2014 03:42, Charles Lepple wrote: > On Oct 1, 2014, at 5:07 AM, Tarvo.Tiits at automaatika.ee wrote: > >> Why is ups.status empty? > If I understand the commit log for SVN r3590, then that revision was supposed to set the status to OL if the numeric state was 9 for ECO mode. It does not appear to do that: > >

Dear all, I need to use R for one estimation, and i have readily available stata command, but i need also the R version of the same command. the estimation in stata is as following: 1. Compute mean values of relevant variables . sum inno lnE lnM Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------

Thanks for help! Are you happy with your Eaton Elipse ECO 1600? Does it work well with ondelay? I will need some UPS that supports the ondelay parameter. Bartosz śr., 10 kwi 2019 o 13:06 Roger Price <roger at rogerprice.org> napisał(a): > On Wed, 10 Apr 2019, Bartosz wrote: > > > I don't see ups.delay.start anywhere. Here is the output of upsrw: > > > > upsrw