similar to: Smooth monotone estimation on R

I would like to create a functional data object where both domain and range are matrices. My erroneous attempts are below. Can you suggest corrections? library(fda) domain<-matrix(c(1,1.1,1.2,2,2.1,2.1,3,3.1,3.2,4,4.1,4.2,5,5.1,5.2), nrow = 5, ncol=3, byrow=TRUE) range<-sin(domain)+matrix(rnorm(15,sd=0.1),nrow=5, ncol=3) myfd<-Data2fd(argvals=domain, y=range) #Error in

Full_Name: David D Degras Version: 2.8.0 OS: Mac OS X Submission from: (NULL) (128.135.239.11) I have recently installed the version 2.8.0 of R along with package fda (v 2.0.2) and its dependencies (including package splines v. 2.8.0). Here are my problems: 1) The package splines should feature functions such a predict.bs, predict.bSpline and such and it does not! I can make calls to bs, ns,

Hello , i am working in the functional time series using the multivariate time series data(hourly time series data). Sir? i am using FAR model more than one order for which no statistical package is available in R, so for this i convert my data into functional form and obtained the functional principle component and from those FPCA i extract their corresponding? FPCscores. Know i use the VAR model

Hello , i am working in the functional time series using themultivariate time series data(hourly time series data). Sir? i am usingFAR model more than one order for which no statistical package is available inR, so for this i convert my data into functional form and obtained thefunctional principle component and from those FPCA i extract theircorresponding? FPCscores. Know i use the VAR model on

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

Good morning,??Please help me to code this code in R. I working in the multivariate time series data, know my objective is that to one year forecast of the hourly time series data, using first five as a training set and the remaining one year as validation. For this??I transform the the data into functional data through Fourier basis functional, apply functional principle components as dimensional

Good morning dear administrators, Please help me to code this code in R. I working in the multivariate time series data, know my objective is that to one year forecast of the hourly time series data, using first five as a training set and the remaining one year as validation. For this??I transform the the data into functional data through Fourier basis functional, apply functional principle

Hello, I am using library fda and I can not run a lot of functions because I receive the error: Error in bsplineS(evalarg, breaks, norder, nderiv) : couldn't find function "spline.des" do you know how I can fix that? Thnaks. Liliana

Dear all, I'm trying to use the Pspline add-on package to fit a quintic spline (norder =3), but I keep running into a Singularity error. > traj.spl <- smooth.Pspline(time, x, norder=3 ) Error in smooth.Pspline(time, x, norder = 3) : Singularity error in solving equations > Playing around with the other parameters produces an "unused arguments" error: > traj.spl

Hi everyone, I have a question. Now I am reading the resource code of the package "ssfcov". The resource code is as following. I cannot find the resource code of the function "myss2d" anywhere in the package. Can anyone give me a hint how to find it in the package. Thanks a lot!!bv > ssfcov function (time, x, subject, nbasis = 5, centered = FALSE, noDiag = TRUE) {

To whom it may concern, I don't know whether this is really a bug with the Pspline package or only a problem with my installation. Things work fine in Linux but not in Mac OS X (Darwin). Both system run the latest public versions of R and Pspline. predict.smooth.Pspline produces only NaN instead of predicted values when norder>2: > library (Pspline) > tt <- seq

Dear all, I need to compute tensor product of B-spline defined over equi-spaced break-points. I wrote my own program (it works in a 2-dimensional setting) library(splines) # set the break-points Knots = seq(-1,1,length=10) # number of splines M = (length(Knots)-4)^2 # short cut to splineDesign function bspline = function(x) splineDesign(Knots,x,outer.ok = T) # bivariate tensor product of

Hey all, I seem to be having trouble fitting a spline to a large set of data using PSpline. It seems to work fine for a data set of size n=4476, but not for anything larger (say, n=4477). For example: THIS WORKS: ----------------------------- random = array(0,c(4476,2)) random[,1] = runif(4476,0,1) random[,2] = runif(4476,0,1) random = random[order(random[,1]),] plot(random[,1],random[,2])

Hi All, I want to fit a piecewise spline of degree 1, i.e. a spline consisting of a straight line over each piece. I downloaded the R package pspline, then I have following questions: 1. in the program, the degree of the spline is specified by 2*norder-1. Why do they adopt such scheme that we can only fit a spline with odd degree? 2. norder cannot be set to 1. Is there any specific reason

Hi, All: What might cause "R CMD check" to report, "could not find function" for a function that has long been in the 'fda' package? Both Jim Ramsay in Ottawa, Canada, and I in San Jose, CA, get this same error. I replicated it with a fresh, anonymous checkout from R-Forge (svn checkout svn://svn.r-forge.r-project.org/svnroot/fda). With this, I did

I have a simple question that irritatingly I haven't been able to figure out on my own. It seems that some functions from the "Pspline" package are successfully installed while others are not. The code with which I'm working is more complicated, but the following highlights my problem. If I run the following code > tt <- seq (0,1,length=20) > xt <- tt^3 > fit

I'm a Beginner on R-language and I' have tried to use some packages but the following message have appeared: "Error in .Fortran("pspline", as.integer(n), as.integer(nvar), as.integer(norder), : C/Fortran function name not in load table" I'm using a Windows version and I'd like to know if it is a installation problem or problems with paths? Thanks,

On 15 February 2017 at 00:57, Daniel P. Berrange <berrange@redhat.com> wrote: > What is the actual error you're getting during startup. # virsh -d0 start instance-0000037c start: domain(optdata): instance-0000037c start: found option <domain>: instance-0000037c start: <domain> trying as domain NAME error: Failed to start domain instance-0000037c error: monitor socket did

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

On 02/15/2017 03:43 AM, Blair Bethwaite wrote: > On 15 February 2017 at 00:57, Daniel P. Berrange <berrange@redhat.com> wrote: >> What is the actual error you're getting during startup. > > # virsh -d0 start instance-0000037c > start: domain(optdata): instance-0000037c > start: found option <domain>: instance-0000037c > start: <domain> trying as