search for: pcurve

...0 -0.083877 0.000090 0.025965 -0.083877 0.002574 and the code as follow: pp <- scan("bmn007.data", quiet= TRUE) x <- matrix(pp, nc=2, byrow=TRUE) fit <- principal.curve(x, plot = TRUE) points(fit,col="red") By now, I got a right result. But when i changed to use pcurve with matrix x as "pcurve(x)", an error was thrown as following: Estimating starting configuration using : CA Error in h %*% diag(sqrt(d)) : non-conformable arguments How to convert a matrix to the format could be accepted by "pcurve"? Any help appreciated! Regards - Sun

Hey, R-listers I am a new user of R and just found the package of PCURVE which can estimate principal curve for arbitrary dimensional data set. Now I have some 2-Dimensional data set X, which is stored as an Nx2 matrix in data.txt file and looks as following: -1.5551 2.4183 1.0051 1.0102 0.90644 0.82163 1.3248 1.7551 -1.3626 1.8568 -1.3554 1.8371 0.039396 0.0015521 -0.9...

A student of mine wants to use R to do some nonmetric multidimensional scaling. According to the R FAQ, there's a package called pcurve that computes multidimensional scaling solutions, but I was not able to locate it the contrib page (I am a Windows user with R version 1.4.1). Can anyone tell me whether it is possible to do nonmetric multidimensional scaling with R, and if so, how? John -----------------------------------------...

...te diamonds dse e1071 ellipse event.chart exactRankTests fastICA fdim fields foreign fracdiff gafit gee geoR gld gregmisc gss ineq leaps lgtdl lmtest locfit logspline lokern lpridge maptree maxstat mclust mda meanscore mgcv mlbench muhaz multiv mvnmle mvtnorm netCDF nlme nlrq norm odesolve oz panel pcurve permax pinktoe pixmap polymars polynom princurve pspline quadprog quantreg rmeta rpart rpvm scatterplot3d sem sgeostat sm sma sn sna splancs sptests spweights strucchange survival tensor tree tripack tseries twostage vegan wavethresh wle xgobi xtable zmatrix -------------- next part --------...

Can I get some advice on what is required to claim a GPL for a package? A package I maintain (pcurve) uses some S code from statlib. I have an email from the author of the S code saying that he would be happy to issue a GPL for the code. Is the email from the author enough, or do I need something more formal? Chris -- Chris Walsh CRC for Freshwater Ecology Water Studies Centre PO Box 23 M...

...ry provides Kruskal's method for nonmetric MDS, sammon in the MASS library provides Sammon's nonlinear mapping method for nonmetric MDS, and cmdscale in the mva library provides metric MDS. A function xgvis in the xgobi package was mentioned but I have not acquired this package as yet. The pcurve package is mentioned in R-FAQ as an MDS procedure, but it is not on the Windows contributed page. Thank you, John -------------------------------------------------------------------- John Miyamoto, Dept. of Psychology, Box 351525 University of Washington, Seattle, WA 98195-1525 Phone 206-543-080...

...and 8.0. I also compiled the contribute packages It The following packages failed, cause I believe I haven't installed all neccessary includes and libaries: checking for package ipred ...failed! checking for package cramer ...failed! checking for package hdf5 ...failed! checking for package pcurve ...failed! checking for package Bhat ...failed! checking for package SASmixed ...failed! checking for package netCDF ...failed! checking for package pear ...failed! checking for package Rmpi ...failed! checking for package tripack ...failed! checking for package spdep ...failed! checking for packag...

...a reproducible example (including your version of R and OS, indicate add-on packages you use, etc.). For your request above, there is nothing much that we could do. After I had found out that this is about a function in the package princurve (do note that there is a second CRAN package pcurve for the same task), all I can say is that R> library(princurve) R> example(principal.curve) R> fit <- principal.curve(x, plot = TRUE, smooth="periodic.lowess", iter= 0) gives some warnings but does not throw an error as you indicated above. This is on Linux with...

hello when I enter: ftable(xtabs(cnt~jeo+slp+pcurv+asp+hey),data=dt) it gives a frequncy table. But some cells are empty. like this: 19 2 72 0 3 68 0 4 99 0 5 110 0 6 97 0 7 74 0 8 46 0

...hecks will know that some packages generate ERROR reports spasmodically. This is normally due to segfaults, and in almost all cases has been traced (eventually) to memory access errors using the tool valgrind (http://valgrind.kde.org). Current examples include RCurl, RandomFields, geoR, kza and pcurve. This week I have run all the examples in all the CRAN packages under valgrind (which took about 120 hours, not all running to completion), and have sent reports to quite a few package maintainers. A good selection of the reports will be available at http://www.stats.ox.ac.uk/pub/bdr/Rvg for th...

Hi, I'm wondering which function would allow fitting this type of data: tmp=rnorm(2000) X.1 = 5+tmp Y.1 = 5+ (5*tmp+rnorm(2000)) tmp=rnorm(100) X.2 = 9+tmp Y.2 = 40+ (1.5*tmp+rnorm(100)) X.3 = 7+ 0.5*runif(500) Y.3 = 15+20*runif(500) X = c(X.1,X.2,X.3) Y = c(Y.1,Y.2,Y.3) plot(X,Y) The problem with loess is that distances for the "goodness of

Sorry, for disturbing the list again. > Also I got one suggestion of using ORDPACK at http://www.netlib.org/. It's ODRPACK at http://www.netlib.org/, not ORDPACK. Best, -- Etsushi Kato ekato at ees.hokudai.ac.jp -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

...uot;random"), smoother = smoothSpline, complexity, vary = FALSE, maxComp, finalCV = FALSE, axis = 1, rank = FALSE, stretch = 2, maxit = 10, trace = FALSE, thresh = 0.001, plotit = FALSE, fitFUN = c("princurve", "pcurve"), latent = FALSE, ...) Docs: function(X, method = c("ca", "pca", "random"), smoother = smoothSpline, complexity, vary = FALSE, maxComp, finalCV = FALSE, axis = 1, rank = FALSE, stretch = 2, max...

...nter=center, scale=scale) if(is.null(rownames(x))) rownames = 1:nrow(x) else rownames = rownames(x) if(is.null(colnames(x))) colnames = paste('V', 1:ncol(x), sep='') else colnames = colnames(x) s = svd(x) s2 = diag(sqrt(s$d[1:n.values])) #s2 = diag(s$d[1:n.values]) pca of pcurve is like this!? switch(match.arg(weight), equal = { g = s$u[,1:n.values] %*% s2 h = s2 %*% t(s$v[,1:n.values]) hl = t(h) }, samples = { g = s$u[,1:n.values] %*% s2 h = t(s$v[,1:n.values]) hl = t(h) }, variables = { g = s$u[,1:n.val...

...tch","orientlib","ouch","outliers","oz","pamr","pan","panel","papply","partitions","partsm","party","pastecs","pbatR","PBSmapping","pcaPP","pcurve","pear","permax","permtest","perturb","pgam","pgirmess","pheno","phpSerialize","PHYLOGR","pinktoe","pixmap","PK","PKfit","plotrix","pls","p...

...t; 3) stop('Please, check the parameters: lambda.ini and lambda.end!') x = as.matrix(x) x = scale(x, center=center, scale=scale) svdx = svd(x) s2 = diag(sqrt(svdx$d[lambda.ini:lambda.end])) # 'prcomp.default' of 'stats' package (and 'pca' of 'pcurve') is like the below! #s2 = diag(svdx$d[lambda.ini:lambda.end]) switch(match.arg(weight), equal = { g = svdx$u[,lambda.ini:lambda.end] %*% s2 h = s2 %*% t(svdx$v[,lambda.ini:lambda.end]) hl = t(h) }, objects = { g = svdx$u[,lambda.ini:lambda.e...