similar to: problem locfit

Hi, i'm trying to apply StepAIC with coxph...but i have the same error: stepAIC(fitBMT) Start: AIC=327.77 Surv(TEMPO,morto==1) ˜ VOD + SESSO + ETA + ........ Error in dropterm.default(fit,scope$drop, scale=scale,trace=max(0, : number of rows in use has changed: remove missing values? anybody know this error?? Thanks. Michele [[alternative HTML version deleted]]

Dear R-help, I am trying to estimate a Cox model with nested effects basing on the minimization of the overall AIC; I have two frailties terms, both gamma distributed. There is a error message (theta2 argument misses) and I don?t understand why. I would like to know what I have wrong. Thank you very much for your time. fitM7 <- coxph(Surv(lifespan,censured) ~ south + frailty(id,

R-help - I'm trying to create axis breaks similar to this : http://www.r-bloggers.com/wp-content/uploads/2010/08/bar-chart-natural-axis-split1.png . Is there a way to do this in R? Here's my code thus far: structure(list(condition = structure(c(2L, 1L, 3L), .Label = c("con", "exp", "unedit"), class = "factor"), trial.avg = c(4.04583333333333,

I have to study censured datas concernig the occurence of infection at the point of insertion of catheter in patients with renal disease. Catheter may be removed for ather reasons than infection, in this case, the observation is censored. Each has exactly 2 observations. The question are the pronostic factors of infection (the other variables are; age, sex, type of renal disease...). Could you

Hey Xen devel team, I’ve never mailed a developer list before, but I’m worried I may have stumbled across a bad Xen/Debian packaging bug. My confidence is somewhere above “just mail -user” and below “open a bug in xen-system-amd64”. Apologies if this is the wrong usage of the list. I’ve updated to Buster, and my Xen kernel no longer boots. I believe it’s because Xen only offers versions 4.11 and

Hi there, this is my first posting to this list, so let me quickly introduce myself. I'm Alex and currently working on a new version of the rsync package for the eisfair Linux distribution?. I have some problems understanding the behaviour of the 'auth users' option in the rsyncd.conf file when running rsync in daemon mode. I set up a module and a secrets file. This is the behaviour

Hi, I have attached a word document to explain the problem i am having creating a for-loop in R with conditions to create a frequency table. I am new to R so any help would be greatly appreciated. Thanks Jones

Hello :-) I want to create a boot/shutdown script including the usual restart, start, status and stop options (on Slackware64 13.1). Netsearching for prior art, I did not find any comprehensive examples. The normal technique seems to be to kill the process recorded in the PID file (let's call it the daemon) but this does not terminate all the rsync server processes. Experimenting and

Full_Name: Somkiat Apipattanavis Version: 2.1.1 OS: Windows Submission from: (NULL) (128.138.44.123) Bug found in predict.locfit for density estimation # Example of bug found in prdict.locfit (Locfit) library('locfit') # generate data y =c(4281,2497,4346,5588,5593,3474,4291,2542,5195,4056, 3114,2864,4904,7625,3377,4001,4999,7191,8062,5668) x1=c( 0.258729, 1.460156, 0.192323,

Apologies for the coming to this late... 1. By now I hope Somkiat has realized that R-bugs is not the place to report problems in contributed packages. Please direct such reports to the package maintainer. 2. This is really user error. predict() expect the newdata to be a data frame containing variables with the same names as those used in the fitting process. E.g., you fitted the model with

Hi, I am running a local logistic regression using locfit. Now, I want to choose the bandwidth using cross-validation. I don't know if there is an additional command to do so or if I can do it in the locfit. I would appreciate any help about this matter. Thank you. Regards, [[alternative HTML version deleted]]

I'm getting the following warning from locfit: Warning: procv: parameters out of bounds Does anyone know what this means? (I haven't found anything about it in the locfit documentation or the list's archives). thanks, Barnet Wagman

Dear All, I am trying to plot a smoothed hazard funciton using the LOCFIT library. Suppose my locfit model is "loc.fit", my plot command looks like this: plot(loc.rw, lty="solid", pch=1, xlab="Age at Risk", ylab="Hazard Rate", ylim=c(0, .5), lwd=2) What's interesting is that whichever number I used in the "pch=" option, I still get a line

Dear memberships, I'm trying to estimate the following multivariate local regression model using the "locfit" package: BMI=m1(RCC)+m2(WCC) where (m1) and (m2) are unknown smooth functions. My problem is that once I get the regression done I cannot get the fitted values of each of this smooth functions (m1) and (m2). My program is the following: library(locfit)

Dear, I am doing some statistics analysis using R software. I am currently trying to implement your library locfit. I am using Win NT and R1.1.1 version obtained on the following Webmaster : http://cm.bell-labs.com/cm/ms/departments/sia/project/locfit/index.html When trying to run it, the following message appears: ********** Error in .C("guessnv", nvm = integer(1), ncm =

The version of locfit on the web site mentioned apparently has been revised by Prof. Loader, and is newer than the CRAN version that I have been maintaining. If Prof. Loader is OK with it, I will take a look and see if I can get the new version into CRAN-conforming form and upload to CRAN. Meanwhile, make sure you're using the package from Prof. Loader's web page, instead of the one on

Hello all! I have recently started using the LOCFIT package, together with Clive Loader's book. I need to implement some method for automatic (plug-in) bandwidth selection in a multivariate kernel regression. From the book, and the LOCFIT documentation, it is not clear whether this is possible. As far as I can see, the only time the various automated procedures, e.g., SJPI, are used is in

hi, after m=locfit(y~x,..........., family=binomial) plot(m,band="local") gives me a plot of locfit's result with a confidence interval around it. i would like to get the actual values that are being used to plot the lines in this figure. i tried using predict, but the standard error it returns when i supply the se=T argument, appears not to be the same as the CI produced by

I just copy parts of my correspondence with the author of Locfit: >>Dear Mr <...>, >>I have tried to apply locfit to a three dimensional data set for density >>estimation. >>The algorithm seems to work well but when I want to plot the result, the >>routine "levelplot" is missing. >>This is not the case in 2-dimensional fitting. >>I=B4d be

Hello!!! Thank you all for your suggestion on my last question about locfit. Now I have another question: I would like to change the evalution structure ev in locfit() to a vector/matrix. I created a vector and a matrix: > vec <- mat.or.vec(101, 1) > vec <- c(0:100) > mat <- matrix(data = vec, nrow = 101, ncol = 1) and I tried to change the ev parameter: > fit <-