similar to: Using mixed models to analyze Longitudinal intervention

Hi, When running z <- plot(fit, age=NA, cholesterol=NA, perim=boundaries, method='image') Legend(z, fun=plogis, at=qlogis(c(.01,.05,.1,.2,.3,.4,.5)), zlab='Probability') And after pointing the cursor to the plot() screen in R, I obtain the following message: Using function "locator(2)" to place opposite corners of image.legend Error in

dear£º in the example£¨nomogram£©£¬I don't understand the meanings of the program which have been marked by red line.And how to compile the program(L <- .4*(sex=='male') + .045*(age-50) + (log(cholesterol - 10)-5.2)*(-2*(sex=='female') + 2*(sex=='male'))). n <- 1000 # define sample size set.seed(17) # so can reproduce the results age <- rnorm(n, 50, 10)

Dear All, using the example from the help of summary.rms library(rms) n <- 1000 # define sample size set.seed(17) # so can reproduce the results age <- rnorm(n, 50, 10) blood.pressure <- rnorm(n, 120, 15) cholesterol <- rnorm(n, 200, 25) sex <- factor(sample(c('female','male'), n,TRUE)) label(age) <- 'Age'

Hi, just making some experiments with design library i get an error if i want plot(fit) - show below from onlineHelp !? ..perhaps is here another mask problem?, but label from xtable which was my first problem is now off ! Thanks for advance & regards, Christian $ n <- 1000 # define sample size $ set.seed(17) # so can reproduce the results $ age <- rnorm(n, 50, 10)

Hello. One of my teaching assistants was experimenting and encountered unexpected behaviour with the lrtest function in the rms package. It appears that when you have a pair of non-nested models that employ an RCS, the error checking for non-nested models appears not to work. Here is a reproducible example. > library(rms) Loading required package: Hmisc Loading required package: lattice

Hello, I'm currently using auto.arima to verify the order of my arima model. I would like to use the model to conduct an intervention analysis. The problem is that, when I include a step function in auto.arima, by including a binary variable in "xreg", the arima order that auto.arima gives is different from when I don't include it. From my understanding, the

Dear everyone, I need R codes to fit an Intervention Time Series Model for two separate policies for a single observed time series. Please, I have prior knowledge of how to fit the impact parameter(omega) but I no not the codes for the decay parameter(delta). Can someone help me? [[alternative HTML version deleted]]

Hello, One of the requirements for our distribution that came down from management (doesn't it always? :) is to have the boot menu on the CD completely optional, but very easy to activate. I'm attaching a patch that does just that. By adding the following lines to the configuration file: MENU ACTIVATETIMEOUT 30 MENU ACTIVATEMESSAGE Press any key for Advanced Menu... The

Hi, Sorry for asking it again. I just thought I asked the question the wrong way. Too much unnecessary details. And therefore, got no replies. I'm configuring a server from, say, a script (automatically). I need it to be done without user saying, "Yes, it's okay to import the key." I was told to run the following command: yum makecache -y --disablerepo=*

Am 12.03.2017 um 16:42 schrieb Yuri Kanivetsky <yuri.kanivetsky at gmail.com>: > > Sorry for asking it again. I just thought I asked the question the > wrong way. Too much unnecessary details. And therefore, got no > replies. > > I'm configuring a server from, say, a script (automatically). I need > it to be done without user saying, "Yes, it's okay to

Am 12.03.2017 um 18:30 schrieb Yuri Kanivetsky: >> what about importing the key beforehand > >> rpm --import keyfile > > That's what I do. But when I try to install a package from that > (passenger) repository manually, yum wants my confirmation: > > # curl --fail -sSLo /etc/yum.repos.d/passenger.repo >

> It looks like the passenger RPM has a script importing the GPG key. Check using "rpm -q --scripts passenger". >From what I can see, it doesn't: # rpm -q --scripts passenger postinstall scriptlet (using /bin/sh): if /usr/sbin/sestatus | grep 'SELinux status' | grep -q enabled; then actual_selinux_variants=; for selinux_policy_file in /etc/selinux/*/policy;

...Check out the full typescript of what happens when installing passenger, please: https://gist.github.com/x-yuri/1dc92db44f89253679ab44f6c3de125c Regards, Yuri

On 13/03/2017 04:38, Yuri Kanivetsky wrote: > ...Check out the full typescript of what happens when installing > passenger, please: > https://gist.github.com/x-yuri/1dc92db44f89253679ab44f6c3de125c > > Regards, > Yuri In my kickstart scripts I call yum with yum -t -y -e 0 This just works for me. -t may be what you are looking for. Tris

Check to see what port the second printer was assigned to... > -----Original Message----- > From: Karl E. Jorgensen [mailto:samba@e-jorgensen.freeserve.co.uk] > Sent: Monday, May 20, 2002 12:41 PM > To: Samba Mailing List > Subject: [Samba] Printer: User Intervention Required - Use Printer > Offline > > > For some reason, I cannot get samba (or Win98, depending on

Hi Team, I am running ARIMAX with TSA package. my code is fit2 <- arimax(yseries, order = c(1,0,1),xtransf = data.frame(X1var),transfer=list(c(1,0))) my question is 1st Q.--> If I need to take difference of X1var then what should i do?. What i am doing like submitting R code as X1vard <- diff(X1var) and then i am including in the xtransf. Same time if i need to take difference of

Dear all, I am wondering how/if it is possible to implement the general methodology of Box and Tiao: "Intervention analysis with applications to economic and environmental problems" (JASA, 1975, pages 70-79) in R? This question has been posted before but without a positive response (at the time): tolstoy.newcastle.edu.au/R/help/06/03/22934.html Perhaps there have been some new

I have a plot produced by function bplot (package = rms) that is really a lattice plot (class="trellis"). It is similar to this plot produced by a very minor modification of the first example on the bplot help page: requiere(rms) n <- 1000 # define sample size set.seed(17) # so can reproduce the results age <- rnorm(n, 50, 10) blood.pressure <- rnorm(n, 120,

Dear R-helpers, I'm having a problem in using plot.design in Design Library. Tho following example code produce the error: > n <- 1000 # define sample size > set.seed(17) # so can reproduce the results > age <- rnorm(n, 50, 10) > blood.pressure <- rnorm(n, 120, 15) > cholesterol <- rnorm(n, 200, 25) > sex <-

Dear Rlist, How can the Legend of the plot.Design() function can be visualized? Following the documentation in R, I did the following (see below), only the 'Legend' function doesn't visualize the legend of the plot (method='image') of the lrmfit. I tried to change par( margin setting) but this didn’t solve it. How can this be solved? Thanks a lot, Jan