[R] convering upper triangular matrix into vector

Hi I have two symmetrical distance matrices and want to compute the correlation coefficient between them (after turning them into vectors). Is there a way of selecting only the upper triangular part of each matrix, then convert this into a vector so I can compute the correlation? Many Thanks

Re: [R] convering upper triangular matrix into vector

[EMAIL PROTECTED] writes: Hi I have two symmetrical distance matrices and want to compute the correlation coefficient between them (after turning them into vectors). Is there a way of selecting only the upper triangular part of each matrix, then convert this into a vector so I can

Re: [R] convering upper triangular matrix into vector

look at: '?upper.tri()' Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://www.med.kuleuven.be/biostat/

[R] arima.sim bug?

Hi, I am using the arima.sim function to generate some AR time series. However, the function does not seem to produce exactly the same time series when I specify the innov parameter. For example r - rnorm(300) x - arima.sim(300, model=list(order=c(1,0,0),ar=c(.96)), innov=r, n.start=10) y -

[R] modeling language for optimization problems

Does anyone know whether R has its own modeling language for optimization problems (like SIMPLE in NuOPT for S-plus)? Paolo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] arima.sim bug?

Kemp S E (Comp) wrote: Hi, I am using the arima.sim function to generate some AR time series. However, the function does not seem to produce exactly the same time series when I specify the innov parameter. For example r - rnorm(300) x - arima.sim(300,

Re: [R] plot question when type = b and pch is a vector

Dear tokas, How about: x - seq(0.01, 10, length = 20) xx - x[7] x[7] - NA plot(c(0, 10), c(-20, 20), type = n, xlab = x, ylab = expression(2 * alpha * log(x))) for(i in 1:4){ lines(x, 2 * i * log(x), lty = 1) text(xx, 2 * i * log(xx), i) } I hope this helps, John

Re: [R] How to get to the varable in a list

Lisa Wang wrote: Hello, I have a list lis as the following: $1 x1 x2 4 3 1 $2 x1 x2 3 3 2 5 3 2 $3 x1 x2 2 3 3 6 3 3 How do I get the x1 varible? for example ss1 I can only guess you mean something like lis[[1]][[x1]] Uwe Ligges

Re: [R] User error (was arima.sim bug?)

On Sun, 2 Oct 2005, Kemp S E (Comp) wrote: Hi, I am using the arima.sim function to generate some AR time series. However, the function does not seem to produce exactly the same time series when I specify the innov parameter. For example r - rnorm(300) x - arima.sim(300,

Re: [R] modeling language for optimization problems

On Sun, 2 Oct 2005, Paolo Cavatore wrote: Does anyone know whether R has its own modeling language for optimization problems (like SIMPLE in NuOPT for S-plus)? No. Note that SIMPLE is the language of NUOPT, not of S-PLUS. There is an (extra-cost) interface module S+NUOPT, but it is an

[R] Size of jpegs/pngs

Dear all I have trouble with setting the size for jpegs and pngs. I need to save a dendrogram of 1000 words into a jpeg or png file. On one of my computers, the following works just fine: bb-agnes(aa, method=ward) jpeg(C:/Temp/test.txt, width=17000, height=2000) plot(bb) dev.off() On my main

Re: [R] p-level in packages mgcv and gam

Denis Chabot said the following on 2005-09-29 21:55: OK, I think I understand better but still have two points to clarify. The first one is about the number of df. I think those who replied on this objected to the way I chose df, not the fact that I would run a model with 7.4 df per se.

Re: [R] Memory management on Windows (was Size of jpegs/pngs)

I think this an issue about the amount of graphics memory. You are asking for an image of about 17*2*3 = 102Mb, and you need more than that. From the help page: Windows imposes limits on the size of bitmaps: these are not documented in the SDK and may depend on the version of

[R] rate instead of scale in ?ks.test

I am not sure whether I'm doing something wrong or there is a bug in the documentation of ks.test. Following the posting guide, as I'm not sure, I haven't found this in the bug tracker, and the R FAQ says that stats is an Add-on package of R, I think this is the place to send it. ?ks.test

Re: [R] p-level in packages mgcv and gam

Thank you very much Henric, now I see! Denis Le 05-10-02 à 11:47, Henric Nilsson a écrit : This test concerns only the non-linear part of the term s(Number, 3). In order to simultaneously test both the linear and non-linear part, as mgcv::summary.gam does, you'd kyp1.1 - gam(Kyphosis ~

[R] generalized linear model for multinomial data?

Dear All: Does R have the package as in SAS's generalized logits model for nominal response data? I have searched but cannot find the windows package. Many thanks, HS __ R-help@stat.math.ethz.ch mailing list

[R] question for some tests

Hi! I come from Poland, and my English is not as good as it should be, so I have to apologize for my mistakes. I would like to ask you for a package, where I can find a test, which is named The Test Of Series (in Poland). In econometric models we use this test to check if the residuals are

Re: [R] question for some tests

On Sun, 2 Oct 2005, Krzysztof Sroka wrote: Hi! I come from Poland, and my English is not as good as it should be, so I have to apologize for my mistakes. I would like to ask you for a package, where I can find a test, which is named The Test Of Series (in Poland). In econometric models we

Re: [R] Access to particular predict value

If you would still like help from this list, please provide a very simple reproducible example to help what you tried that didn't quite work. There are several functions that use inverse distance weighting (IDW) and kriging for geostatistical computations. spencer graves

Re: [R] generalized linear model for multinomial data?

Dear Hongyu, See multinom() in the nnet package, associated with Venables and Ripley's Modern Applied Statistics with S. John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] R-code for binormla distribution

What are you trying to do that requires binormal probabilities other than pmvnorm? spencer graves Nabil Channouf wrote: Dear users, does any one have a code (S or R) to compute the binormal distribution (or the upper its quadrant area) other than the pmvnorm. Thanks

Re: [R] generalized linear model for multinomial data?

Dear Professor: Thanks! I just found it. It is bundled within the VR package. Hongyu - Original Message - From: John Fox [EMAIL PROTECTED] To: 'Hongyu Sun' [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Sent: Sunday, October 02, 2005 3:27 PM Subject: RE: [R] generalized linear model

[R] Selecting columns of a table (not consecutive)

Hi everyone, I have a table with 10 columns and many rows. I want to select the 3rd,4th and 10th column. If I wanted to select the 3rd and 4th it would be no problem. q3[ ,(3:4)] but how do I select the 10th along with it? Thanks for your help, George p.s Please respond if you know,

Re: [R] Selecting columns of a table (not consecutive)

George Papayiannis wrote: Hi everyone, I have a table with 10 columns and many rows. I want to select the 3rd,4th and 10th column. If I wanted to select the 3rd and 4th it would be no problem. q3[ ,(3:4)] but how do I select the 10th along with it? Thanks for your help, George p.s Please

[R] What is Mandel's (Fitting) Test?

Hello everyone, A little background first: I have collected psychophysical data from 12 participants, and each participant's data is represented as a scatter plot (Percieved roughness versus Physical roughness). I would like to know whether, on average, this data is best fit by a linear function

Re: [R] R-code for binormla distribution

Why do you say we do not have the package that contains pmvnorm? It is in library(mvtnorm), which you can get via CRAN (the Comprehensive R Archive Network). If you are using Windows or a standard Linux of MacOS, install.packages('mvtnorm') should install it on your computer (unless

Re: [R] Local install of a contributed package under Linux (was: R-code for binormla distribution)

There probably is a way under Linux to install a contributed package in your local directory. I don't know how, but if you read the documentation on download.packages and install.packages, you might be able to do it without bothering your Info Tech department. I know there is a

Re: [R] Local install of a contributed package under Linux (was: R-code for binormla distribution)

Try R CMD INSTALL -l lib pkgs The help file is in the utils package. I'm sure this is documented in the manual too. Jon -- Jonathan Baron, Professor of Psychology, University of Pennsylvania Home page: http://www.sas.upenn.edu/~baron __

[R] grob questions

If I run the following example from: http://www.stat.auckland.ac.nz/~paul/grid/doc/grobs.pdf grid.newpage() pushViewport(viewport(w = 0.5, h = 0.5)) myplot - gTree(name = myplot, children = gList(rectGrob(name = box, + gp = gpar(col = grey)), xaxisGrob(name = xaxis))) grid.draw(myplot)