similar to: constrained splines in GLM

Hello, I have a big dataframe with *NO* na's (9 columns, 293380 rows). # doing memory.limit(size = 1000000000) cor(x) #gives Error in cor(x) : missing observations in cov/cor In addition: Warning message: NAs introduced by coercion #I found the obvious workaround: COR <- matrix(rep(0, 81),9,9) for (i in 1:9) for (j in 1:9) {if (i>j) COR[i,j] <- cor (x[,i],x[,j])} #which works fine,

Hello I am looking for a function to create an Integrated (I) spline basis, somehting similar to the likes of 'bs' and 'ns'. I have come across the funcitons, fda::eval.monfd Values of a Monotone Functional Data Object fda::/.fd FDA internal functions fda::monfn Evaluates a monotone function fda::smooth.monotone Monotone

Hello everyone, Hope it is the good place for this (I discuss the question of the right place below). Most of the time, warnings are more than 1000 characters long and thus are truncated. Most of the time, this generates printouts with unbalanced parentheses. Intelligent text editors which do parentheses highlighting get very confused with this. After too many warnings, they give errors, and

*Paid internship "web application development with R"* The main skills for this proposed paid internship are: R, Shiny, and components used by Shiny (plotly, javascript). Other skills (Debian, Docker, CKAN...) are useful. This internship offer is part of a partnership between the Grenoble Alpes Cybersecurity Institute and the CYBIS project. The *Cybersecurity Institute* aims to

Hi everyone. There's been several threads on baseline hazard in Cox model but I think they were all on cumulative baseline hazard, for instance http://tolstoy.newcastle.edu.au/R/help/01a/0464.html http://tolstoy.newcastle.edu.au/R/help/01a/0436.html "basehaz" in package survival seems to do a cumulative hazard. extract from the basehaz function: sfit <- survfit(fit) H

Hello, coxph does not use any information that are in the dataset between event times (or "death times") , since computation only occurs at event times. For instance, removing observations when there is no event at that time in the whole dataset does not change the results: > set.seed(1) > data <- as.data.frame(cbind(start=c(1:5,1:5,1:4),stop=c(2:6,2:6,2:5),status=c(rep(

I get the following error: Error in backSpline.npolySpline(sp) : spline must be monotone Has anyone had this error before? any ideas on a workaround? > > vols=read.csv(file="C:/Documents and Settings/Hugh/My > Documents/PhD/Swaption vols.csv" + , header=TRUE, sep=",") > X<-ts(vols[,2]) > #X > > > dcOU<-function(x,t,x0,theta,log=FALSE){ +

Hello everyone, I am searching for a covariate selection procedure in a cox model formulated as a counting process. I use intervals, my formula looks like coxph(Surv(start,stop,status)~ x1+x2+...+cluster(id),robust=T) where id is a country code (I study occurence of civil wars from 1962 to 1997). I'd like something not based on p-values, since they have several flaws for this purpose. I turned

Hi, this is a newbie question. Would it be able to convert e.g. annual time series of flow data (or an index series) into quarterly data using interpolation splines by means of an existing R-function? The problem is, that the average value of the computed quarterly values must be the annual value, i.e. the spline should cross the annual values (in a stairs-line plot) in the middle of every annual

I'm sure I'm doing something completely boneheaded here, but I've used this idiom (constructing an interpolation spline and using prediction from a backSpline to find an approximation profile confidence interval) many times before and haven't hit this particular problem: r2 <- c(1.04409027570601, 1.09953936543359, 1.15498845516117, 1.21043754488875,

Hi, Can anybody tell me what the message below means and how to overcome it. Thanks, Neil Warning message: X matrix deemed to be singular; variable 2 in: coxph(Surv(age_at_death, death) ~ project$pluralgp + project$yrborn + ......... >

The following creates a point process version of a sinewave (maybe there's a better way): p<-amp*cos(2*pi*freq*(1:n)/n ) + 0.5 as.numeric(runif(1:n)<p) I get something like this: 1 o oo oo o o ooooo o 0 ooo o o o o o ooo ooo In case it's not obvious, this is a noisy version of 1 ooooooo ooooooo 0 ooooooo

Dear list, As can be seen below, adding a NA to the y values in a call to isoreg results in a segfault. ir4 <- isoreg(1:10, y4 <- c(5, 9, 1:2, 5:8, NA, 8)) Adding missing values to the x values, on the contrary, gives an error, but maybe the error message could be tailored to this particular situation. y <- c(5, 9, 1:2, 5:8, 3, 8) x <- c(1:9, NA) isoreg(x, y) ## error message:

Hi all, I'm looking for smooth monotone estimation packages, preferably using splines. I downloaded the 'cobs' package and intend to use it, but since it offers only quadratic splines based on L1 minimization, I'd like to compare its performance to that of a more 'mainstream' cubic-spline, L2-norm minimizing spline. Preferably a smoothing spline. Does anyone know of such

I believe there is a memory leak in isoreg in the current version of R, as I believe the following shows > gc() used (Mb) gc trigger (Mb) max used (Mb) Ncells 120405 3.3 350000 9.4 350000 9.4 Vcells 78639 0.6 786432 6.0 392463 3.0 > for(k in 1:100) { + + y <- runif(10000) + isoreg(x,y) + } > rm(x) > rm(y) > gc() used (Mb) gc

I have developed a GAM model in order to predict Y using 4 X variables. 2 of these X's are factors, and 1 is a spline. Part of the data looks like: Days WRM variety PWM O_EC 31 75 1 90 234 31 79 1 78 283 31 82 1 92 281 31 84 1 96 213 31 99 2 69 247 31 100 2 77 324 31 104 2 74 259 31 118 2 81 282 31 61 3 58 478 31 98 3 83 429 31 98 3 70 379 31 156 3 87 467 31 78 4 56 283 31 97 4 67 282 31

Full_Name: Petr Klasterecky Version: 1.8.1 OS: Windows XP, Linux Submission from: (NULL) (195.113.27.212) The isoreg() function causes R to crash when called repeatedly. Consider the following simple script: { library(modreg) N <- 10 x <- rnorm(N) print("Original x values:") print(x) for(n in (1:N)){print(y <- isoreg(x[1:n])$yf)} } I am able to run (call) it several

Running R : Copyright 1999, The R Development Core Team Version 0.90.1 (December 15, 1999) on NT 4.0 gives me problems with: win.metafile(file="./x.emf") x <- 1:100/7 plot(x,cos(x),type="n") lines(x,sin(x)) abline(v=0:15,h=-2:2/2,col=gray(0.8)) Only labels and titles on axes are in the file. No axes or lines of any kind. (I look at the file by inserting it in Word,

Once you hae in stalled Guido's NT port of R it is nice and easy to install the add-on packages too. However, it would be nice to have links to the documentation of the functions in all the packages from the "All Installed Functions" html-page, and not have to go through each of the packages. Is there a utility somewhere to do this? Bendix \\\|///

Dear Duncan, thank you very much, that is it! Indeed, there was activated in my internet router the "telekom navigation help". When I am (or the system is) looking for an not existent url, this "help" opens its own website. After deactivation of telekom navigation help, the normal message (server not found) comes. I never gave attention to this "help", although it is