similar to: How to use a matrix in pcurve?

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.99622 0.99245 0.2314

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

Hi All: Package TSERIES is stated to depend on MVA. However, there is no MVA package to be found under the list of package sources. Best wishes, ANDREW tseries: Package for time series analysis Package for time series analysis with emphasis on non-linear and non-stationary modelling Version: 0.7-6 Depends: ts, mva, quadprog Date: 2001-08-27 Author: Compiled by Adrian

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

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

Dears, I've been learning biplot (Gabriel, 1971) and I found the function 'biplot', inside of the package 'stats', useful but, a bit limited. So, I'm thinking to start a colaborative package to enhance this methods to other multivariate methods. In this way, I would like to start it, making public a new function (biplot.pca, still in development, but running) that make

Those of you who monitor the daily package checks 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

I sent a query to R-Help about the availability of nonmetric multidimensional scaling (MDS) algorithms in R. I would like to thank Tony Rossini, Jonathan Baron, Sundar Dorai-Raj, and Brian Ripley for helpful replies. The gist of the replies is that isoMDS in the MASS library provides Kruskal's method for nonmetric MDS, sammon in the MASS library provides Sammon's nonlinear mapping method

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, during today I build R 1.5.0-packages on a SuSE 7.3 they normaly should also run on a 7.2 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

Moved to R-help from R-SIG-gR. On Thu, 21 Oct 2004 17:44:27 +0800 XP Sun wrote: > hi, everyone, > > this is the first time for me to use R, > would you like to tell me how to use this command? thank you in > advance. > > fit <- principal.curve(y, plot = TRUE, smooth="periodic.lowess", > iter=0) > > when i uesd the parameter > >

Dear all R-users, I wanted to fit a Garch(1,1) model to a dataset by: >garch1 = garch(na.omit(dat)) But I got a warning message while executing, which is: >Warning message: >NaNs produced in: sqrt(pred$e) The garch parameters that I got are: > garch1 Call: garch(x = na.omit(dat)) Coefficient(s): a0 a1 b1 1.212e-04 1.001e+00 1.111e-14 Can any one

Dear list, I have been checking my package ('analogue') using R2.13-0-alpha (details of exact svn version appended below) and the R CMD check procedure is generating an error rebuilding a vignette in the package, which raises a NOTE in the check. The log printed to screen during check shows: .... * checking re-building of vignettes ... NOTE Error in re-building vignettes: ... Error in

Hi all, what indicates the parameter theta in the summary of a random effect panel model estimated with the plm function? example: data("Produc", package = "plm") zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, model="random", data = Produc, index = c("state","year")) summary(zz) Effects: var std.dev

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

Dear all, I've been learning biplot (Gabriel, 1971) and some days ago I sent for this list a procedural function with invitation for a collaborative package. Jari Oksanen made some suggestions and I agree with all. So, I reworked the function under the object-oriented programming (OOP/S3). I think it is now a good frame for more resources. Below it is the function and a small script to

I need to remove certain interactions and keep only the one between the second level of the factor and the continuous var t2 bin4 <- glm(resp2~ t*t2+c5.vrm,data=dfa,family="quasibinomial") > summary(bin4) Call: glm(formula = resp2 ~ t * t2 + c5.vrm, family = "quasibinomial", data = dfa) Deviance Residuals: Min 1Q Median 3Q Max -6.5464

I am working with an existing schema and thus limited in how much I can bend my existing column names. I have a table Users, and a table ForumUsers, for which I have a "has_one" relationship - a user has_one forumuser. When I look at my user records I wish to include some information from forumuser: @users = User.find :all, :include => :forum_user, :order =>

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

Hi there. Have this code, and i need to find put why it fails. >>Controller ------------ def update @user = User.find(@params[:id]) if @user.update_attributes(@params[:user]) flash[:notice] = ''The user was successfully updated.'' redirect_to :action => ''show'', :id => @user else flash[:notice] = "Error: The user