search for: ol

Dear all, When fitting an "ols.model", the confidence interval at 95% doesn't cover the plotted data points because it is very narrow. Does this mean that the model is 'overfitted' or is there a specific amount of serial correlation in the residuals? Which R functions can be used to evaluate (diagnostics) maj...

Thanks to Gabor for setting me right. My code is as follows. I found it useful for learning optim(), and you might find it similarly useful. I will be most grateful if you can guide me on how to do this better. Should one be using optim() or stats4::mle? set.seed(101) # For replicability # Setup problem X <- cbind(1, runif...

I'm running nut 2.0.1 with an MGE pulsar evolution 3000 with the mge-utalk (serial) driver on NetBSD 3.1 . I had a very strange occurance in my basement datacenter this early AM. the following was logged by upslog: 20080219 023452 100.0 123.0 049.0 [OL CHRG] NA 60.00 20080219 023522 100.0 122.8 049.0 [OL CHRG] NA 60.00 20080219 023552 098.0...

Hi! I need to run ols regressions with Huber-White sandwich estimators and the correponding standard errors, without an intercept. What I'm trying to do is create an ols object and then use the robcov() function, on the order of: f <- ols(depvar ~ ind1 + ind2, x=TRUE) robcov(f) However, when I go f <- o...

...h the driver output and the last few sections of truss output leading to the hang. Both outputs end at the CTRL+C i pressed when i was forced to end the processes. Hope this is helpful. Please bear with me if i didn't run it with the right options etc - as I mentioned, I'm a little new to Solaris. Regards, Richard On Sun, Apr 5, 2015 at 1:09 AM Richard Flint <richard.flint at gmail.com> wrote: > Thank you for the rapid response. I will try and investigate getting > answers to some of your points but I'm a little new to Solaris so I'll need > some time. Glancing...

Hi everyone, I'm having some difficulty getting "simple effects" for the ols() function in the rms package. The example below illustrates my difficulty -- I'll be grateful for any help. #make up some data exD <- structure(list(Gender = structure(c(1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L), .Label = c("F", "M"), class = "factor...

Thank you for the rapid response. I will try and investigate getting answers to some of your points but I'm a little new to Solaris so I'll need some time. Glancing at the configure output, it looks like it built against v0.1.7 of libusb (yes i think that is derived from the one you mention), checking for libusb version via pkg-config... 0.1.7 found checking for libusb cflags... checking for libusb ldflags... -lusb che...

Firefox used to render <ol type="a"> <li>.................</li> as a. ............... Now it does not. The HTML 4.1 spec at http://www.w3.org/TR/html401/struct/lists.html#type-values says use CSS instead. These <ol type="a"> <ol type="A"> <ol type="i...

I am trying to extract significance levels out of a robcov+ols call. For background: I am analysing data where multiple measurements(2 per topic) were taken from individuals(36) on their emotional reaction (dependent variable) to various topics (3 topics). Because I have several emotions and a rotation to do on the topics, I'd like to have the results pu...

Peter, Here are the transitions that I have tested with your latest changes (on the APC/Back-UPS ES 650): 1) OL ==> OB DISCHRG followed by OB DISCHRG ==> OL (no CHRG inbetween) 2) OL CHRG ==> OL 3) OL CHRG ==> OL CHRG LB (by manipulating battery.charge.low) 4) OL CHRG LB ==> OL CHRG (i.e. transition from below) 5) OB DISCHRG ==> FSD OB DISCHRG LB followed by proper shutdown In other words,...

.... I build a full model consisting of five predictors with the use of penalization by rms package (lrm, pentrace etc) because of events per variable issue. Then, I tried to approximate the full model by step-down technique predicting L from all of the componet variables using ordinary least squares (ols in rms package) as the followings. I would like to know whether I am doing right or not. > library(rms) > plogit <- predict(full.model) > full.ols <- ols(plogit ~ stenosis+x1+x2+ClinicalScore+procedure, sigma=1) > fastbw(full.ols, aics=1e10) Deleted Chi-Sq d.f. P Re...

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...

...the error generated in R seems to imply it cannot locate a file. I've checked my paths and they are correct. Also, my data is loading correctly. Many thanks, Joe ######################################################### # R code # Runs Bayesian Ordered Logit by calling WinBUGS from R # ologit2.txt: WinBUGS commands library(R2WinBUGS) setwd("c:/docume~1/admini~1/mydocu~1/r_tuto~1") load("oldat.Rdata") # R data file containing data frame ol.dat # with vars: q02, bf23f, bf23b, bf22, bf34a, bf34.1, bf34.2 q02 <- ol.dat$q02 bf23f &...

How can I extract the adjusted R^2 value from an ols object (using rms package)? library(rms) x <- rnorm(10) y <- x + rnorm(10) ols1 <- ols(y ~ x) Typing "ols1" displays adjusted R^2 among other things, but how can I assign it to a variable? I tried str(ols1) but couldn't see where to go from there. Thanks, Mark Seeto

Hi, I'm new to this list and I have a question concerning using the genericups driver (type 22) with the Cyberpower 485 in serial mode. When running the driver in debug mode: genericups -a myupsname_genericups22 -DDDDDD I observe the following statuses as the mains power dies, switching over to battery power, and then letting it discharge. Initially, with mains power on, status = OL. When mains power goes off, and we run on battery, status = OB. When battery level drops sufficiently the UPS starts beeping, and status = OB OL....

Dear all, I have tried to merge two cross sections (crossA and crossB) as follows, both crossA and crossB include the same variables, e.g., ID,y,x1,x2: Panel <- merge(crossA,crossB,by="ID") I have, thereafter, tried to run the following pooled OLS: OLS <- lm(Panel$y~Panel$x1+Panel$x2) and have received the following error message "Error in model.fra...

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...

I cleaned up my old benchmarking code and added checks for missing data to compare various ways of finding OLS regression coefficients. I thought I would share this for others. the long and short of it is that I would recommend ols.crossprod = function (y, x) { x <- as.matrix(x)...

Hi folks! Does anybody know how to use the OLS function in fMultivar to do dynamic regression? I've tried specifying lags in OLS using a data series created in fSeries and it doesn't seem to work. I've done dynamic regression using dyn$lm and I was wondering how to accomplish the same th...

I might be sorry for asking this question :-) I have two equations and I tried to estimate them individually with "lm" and "gls", and then in a system (using systemfit) with "OLS", "WLS" and "SUR". Quite surprisingly (for myself at least) the results are identical to the last digit. Could someone (please!) give a hint as to what am I doing wrong? Thanks, mihai --------------------------------- [[alternative HTML version deleted]]