similar to: Smooth monotone estimation on R

Hi all, I think this question could be quite trivial, but I can?t find out the solution... How can you compute the statistic "mode" of a sample, in case it exists (as mode() returns the mode of an object)? I tried help.search("mode") but I couldn't find a clue... Any help would be much appreciated. Regards, Aurora

Hi all, I want to create a rather standard overlaid qqnorm plot on a single variable, with different subgroups of the same dataset plotted using different colors/symbols/etc. (I don't want side-by-side, rather different-colored curves on the same graph) I managed to do it rather tediously using "lapply" and "split", and wondered whether there is any single-command

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){ +

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,

Hello, i have a little problem with R and i hope you can help me. I want to use splines to estimate a function but i want to force the interpolation to be monotone. Is this possible with R ? Thank you, Rémi. --------------------------------- [[alternative HTML version deleted]]

Dear R users, I have repeated measurements on individuals. I want to estimate the time-varying effect of a factor variable X (taking three levels), e.g. a model in the spirit of Hastie and Tibshirani (1993). I am considering using the package "mgvc" which implements generalized additive models, especially the function gamm, which estimates generalized additive mixed models, and thus,

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

Dear R-Users, I have problems with the cobs library. When doing the cobs example, I get the folling error message: example(cobs) cobs> x <- seq(-1, 3, , 150) cobs> y <- (f.true <- pnorm(2 * x)) + rnorm(150)/10 cobs> con <- rbind(c(1, min(x), 0), c(-1, max(x), 1), c(0, 0, 0.5)) cobs> Rbs <- cobs(x, y, constraint = "increase", pointwise = con)

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, I have been trying to take the derivative of a quadratic B-spline obtained by using the COBS library. What I would like to do is similar to what one can do by using fit<-smooth.spline(cdf) xx<-seq(-10,10,.1) predict(fit, xx, deriv = 1) The goal is to fit the spline to data that is approximating a cumulative distribution function (e.g. in my example, cdf is a 2-column matrix with x

Hi all I us Linux for some two years now, most of the time i us ubuntu ultimate edition. since version 1.2 of wine I could also run Solidworks 2009. The installation was made by the instruction on the AppDB: http://appdb.winehq.org/objectManager.php?sClass=version&iId=18101 Now I an trying to us Linux Mint 9 ( ubuntu 9.10 ) I installed wine 1.2 and tried to install Solidworks like I

Hi All, Did anyone came across a 411 XML service I can feed to the "service" button with XML? Some other feed I can manipulate to XML query? Assaf Benharoosh -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20041004/dbf552ac/attachment.htm

For some reason I can't get the Dial command 'C' flag to work. The calls are recorded in the CDR with the 'C' on. Does anyone have an idea? extensions.conf: exten => 114,1,Dial(SIP/114,,C) It shows in the lastapp: cdr: | 2004-10-05 13:16:02 | "112" | | 114 | intern | SIP/112-3fb6 | SIP/114-0e7a | Dial | SIP/114||c | 6 |

Hi all, Does anyone know how to get the channel ID on the other side of the call? For example: When SIP/50 calls SIP/21, and the call is answered by SIP/21 I get: SIP/21-6735 answered SIP/50-b456 ${CHANNEL} will show me SIP/50-b456. Is there a parameter or a workaround to get the SIP/21-6735 part? Thanks. Assaf Benharoosh -------------- next part -------------- An HTML attachment was

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

I've noticed some problems with solidworks when using Enterprise Web Hosting (http://www.voipreview.org/enterprise-web-hosting) with my Small Business Web Hosting (http://www.voipreview.org/small-business-web-hosting)... After a couple of trying to redownload and messing around with it a couple of times, I can get it to work.

I'm sure there must be various peak-finding algorithms out there. Not knowing of any, I have written one myself*, but I thought I'd ask to see what's out there. Basically, I have a 2-dimensional data set and I want to identify local peaks in the data, while ignoring "trivial" peaks. My naive algorithm first identifies every peak and valley (point of inflection change

Hello, We would like to apply the smooth monotone function to our data which correspond to a non-linear function. We follow the example posted on the web, but in our case it did not apply. We always get a straight line in response. Which parameters we should change. ind.basis = create.bspline.basis(c(min(time),max(time)),nbasis=38,norder=4) Wfdob =

Hello List. I am having a very interesting problem with my samba server printer setup. Some background , the samba server is Suse Enterprise 9 (x86_64 platform) , we are running samba-3.0.26a-0.7 , with a connection to an openldap and sharing the printers via cupsd . File sharing is working properly , so it's smooth sailing on that end , but the printing is acting up . some users are able

I am having an issue where a feature works fine but the corresponding test fails, and I can''t understand why. I created a mock application to demonstrate the problem. https://github.com/assafshomer/clear_button All it does is show an input form with a clear button. The clearing happens by redirecting back to the same form. It does clear it when you run the app, but the test for it