search for: isoreg

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

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 times, let say 6 o...

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) + } >...

Before re-inventing the wheel I would like to ask: does anyone know about an optimizer in R which can reliably identify which value of X (Xopt) leads to Y (Yopt) closest to Ytarget in Y <- MonotoneStepFun(X) optionally with the restriction that Yopt <= Ytarget (at least if any Y <= Ytarget, otherwise any Yopt > Ytarget would be the preferred answer) If none is known, I will write

Full_Name: Benjamin Tyner Version: 2.6.0 (43063) OS: WinXP Submission from: (NULL) (171.161.224.10) At the end of the examples for isoreg, there is cat("R^2 =", formatC(sum(residuals(ir4)^2) / (9*var(y4)), digits=2),"\n") I think this should be cat("R^2 =", formatC(1 - sum(residuals(ir4)^2) / (9*var(y4)), digits=2),"\n" Thanks Ben platform i386-pc-mingw32 arch...

Hello, If I want to see how, say, apply function is written, how would I be able to do that? Just typing "apply" at the prompt does not work. Thank you for help! Sergey

...plot.HoltWinters* plot.POSIXct* [4] plot.POSIXlt* plot.TukeyHSD plot.acf* [7] plot.data.frame* plot.decomposed.ts* plot.default [10] plot.dendrogram* plot.density plot.ecdf [13] plot.factor* plot.formula* plot.hclust* [16] plot.histogram* plot.isoreg* plot.lm [19] plot.medpolish* plot.mlm plot.ppr* [22] plot.prcomp* plot.princomp* plot.profile.nls* [25] plot.spec plot.spec.coherency plot.spec.phase [28] plot.stepfun plot.stl* plot.table* [31] plot.ts plot.tskernel* No...

Hi, I am trying to use R to mimic what I did in SAS. proc transreg data=x ; model identity(GSI)=monotone(group1); output out=d2 pprefix=M; run; Accroding to SAS documentation, the MONOTONE transfomation algorithm comes from (Kruskal 1964, secondary approach to ties). I have tried ace. it does provide some kind of monotone transformation, but it is not what I expected. Here is how

Hi, I'm trying to find an implementation of monotonic regression in R and I haven't been able to find anything that's really related to this. isoMDS in the MASS package uses monotonic regression, however, I was wondering if there is any standalone function for monotonic regression? Basically what I'm trying to do is implement monotonic regression where I can see not just the

...portsrc.f -o portsrc.o cc -xtarget=generic64 -G -L/mounts/devel/SUNWspro/lib/amd64 -L/usr/sfw/lib/amd64 -L/mounts/devel/GNU/repoz/readline43/lib/amd64 -o stats.so init.o kmeans.o ansari.o bandwidths.o chisqsim.o d2x2xk.o fexact.o kendall.o ks.o line.o smooth.o prho.o swilk.o ksmooth.o loessc.o isoreg.o Srunmed.o Trunmed.o dblcen.o distance.o hclust-utils.o nls.o HoltWinters.o PPsum.o arima.o burg.o filter.o mAR.o pacf.o starma.o port.o family.o bsplvd.o bvalue.o bvalus.o loessf.o ppr.o qsbart.o sbart.o sgram.o sinerp.o sslvrg.o stxwx.o hclust.o kmns.o eureka.o stl.o portsrc.o -xlic_lib=s...

...ing the latest version of R so that I can avoid these errors? configure make ... ... ... f90: CODE: 0 WORDS, DATA: 0 WORDS gcc -G -L/usr/local/lib -o stats.so init.o kmeans.o ansari.o bandwidths.o chisq sim.o d2x2xk.o fexact.o kendall.o ks.o line.o smooth.o prho.o swilk.o ksmooth .o loessc.o isoreg.o Srunmed.o Trunmed.o dblcen.o distance.o hclust-utils.o nl s.o HoltWinters.o PPsum.o arima.o burg.o filter.o mAR.o pacf.o starma.o port.o family.o sbart.o bsplvd.o bvalue.o bvalus.o loessf.o ppr.o qsbart.o sgram.o si nerp.o sslvrg.o stxwx.o hclust.o kmns.o eureka.o stl.o portsrc.o -L../../...

...a.frame" "plot.default" "plot.density" [7] "plot.design" "plot.factor" "plot.formula" [10] "plot.function" "plot.hclust" "plot.histogram" [13] "plot.isoreg" "plot.lm" "plot.mlm" [16] "plot.mts" "plot.new" "plot.ppr" [19] "plot.spec" "plot.spec.coherency" "plot.spec.phase" [22] "plot.tab...

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

Hi, from the man page ?methods, I expected to be able to build pairs (class,package) for a given S3 method, e.g., print, using attr(methods(print), 'info'). However all the methods, except the ones defined in base or S4 methods, get the 'from' value "registered S3method for print", instead of the actual package name (see below for the first rows). Is this normal

...xample - stats.so currently the command to make stats.so is: gcc -Wl,-bnoquiet -std=gnu99 -shared -Wl,-brtl -Wl,-G -Wl,-bexpall -Wl,-bnoentry -o stats.so init.o kmeans.o ansari.o bandwidth s.o chisqsim.o d2x2xk.o fexact.o kendall.o ks.o line.o smooth.o prho.o swilk.o ksmooth.o loessc.o monoSpl.o isoreg.o Srunmed.o d blcen.o distance.o hclust-utils.o nls.o rWishart.o HoltWinters.o PPsum.o arima.o burg.o filter.o mAR.o pacf.o starma.o port.o fa mily.o sbart.o approx.o loglin.o lowess.o massdist.o splines.o lm.o complete_cases.o cov.o deriv.o fft.o fourier.o model.o optim .o optimize.o integrate.o r...

On 2015-10-15 15:02, Prof Brian Ripley wrote: > On 15/10/2015 13:32, Michael Felt wrote: >> Hi. >> >> Just wanted to let you know I am getting close to packaging R for AIX in rephrase - would like to be active in keeping R binaries current for AIX. My interest in not in R per se (rather a colleague who has a project that uses R, so I hope to assist him, and others like him).

...lot.dendrogram* plot.density [10] plot.ecdf plot.factor* plot.formula* [13] plot.function plot.gam plot.gls* [16] plot.hclust* plot.histogram* plot.HoltWinters* [19] plot.intervals.lmList* plot.isoreg* plot.lm [22] plot.lme* plot.lmList* plot.medpolish* [25] plot.mlm plot.nffGroupedData* plot.nfnGroupedData* [28] plot.nls* plot.nmGroupedData* plot.pdMat* [31] plot.ppr* plot.prcomp*...

With the recent changes in R 1.9.x, it is now impossible to properly call plot on a formula() where the formula is provided via a named first argument. On today's R-1.9.1-alpha: > x <- 1:10 > y <- rnorm(x,0.25) > > plot(x~y) > > plot(x=x~y) Error in terms.formula(formula, data = data) : argument is not a valid model > > plot(formula=x~y) Error in

...veral variables that I need concurrently in the same model. I looked at package mgcv but I do not know if/how I can use it in GLM (not GAM) : I could not manage to adapt the mono.con(mgcv) example to GLM. The help for package fda is not complete. Not sure that backSpline(splines) does what I need. isoreg (modreg) seems to do univariate regressions. Some of my covariates are linear. Three covariates (x1,x2 and x3) must be transformed in a decreasing and convex way like this: |o |o | o | o | o | o | ooooo |----------------- Currently, I use exp(-x1/alpha1)+exp(-x2/alpha...

...ail end of make follows. Does anyone have any experience with that? Thanks a lot in advance for any suggestions. gcc -std=gnu99 -shared -L/usr/local/lib -o stats.so init.o kmeans.o ansari.o bandwidths.o chisqsim.o d2x2xk.o fexact.o kendall.o ks.o line.o smooth.o prho.o swilk.o ksmooth.o loessc.o isoreg.o Srunmed.o Trunmed.o dblcen.o distance.o hclust-utils.o nls.o HoltWinters.o PPsum.o arima.o burg.o filter.o mAR.o pacf.o starma.o port.o family.o sbart.o approx.o loglin.o lowess.o massdist.o splines.o bsplvd.o bvalue.o bvalus.o loessf.o ppr.o qsbart.o sgram.o sinerp.o sslvrg.o stxwx.o hclust.o km...