[R] procrustes results affected by order of sites in input file

2010-12-01 Thread Christine Dolph
Dear All, I am using a Procrustes analysis to compare two NMDS ordinations for the same set of sites. One ordination is based on fish data, the other is based on invertebrate data. Ordinations were derived using metaMDS() from the {vegan} library as follows: fish.mds<-metaMDS(fish.data, distance=

Re: [R] Font family not found in Windows font database

2010-12-01 Thread Mark Ebbert
That did it. Thanks! On Dec 1, 2010, at 12:38 PM, Paul Murrell wrote: > Hi > > On 2/12/2010 4:55 a.m., Mark Ebbert wrote: >> Dear R Gurus, >> >> I have a fairly simple problem, but I haven't been able to find the >> answer on 'the google' or in the r-help archives. >> >> I am generating plots

Re: [R] parametric estimators for species richness in R

2010-12-01 Thread Dennis Murphy
Hi: When the question is: 'Are there any functions in R for ?', package sos is a good place to start: library(sos) # install if necessary findFn('MMeans') Well actually, MMMeans and MMruns came up dry, but Michaelis-Menten is a commonly occurring model in the biosciences: findFn('Michaelis-

Re: [R] draw categorical histogram

2010-12-01 Thread phoebe kong
This would not draw a histogram. On Wed, Dec 1, 2010 at 11:54 AM, Andrew Miles wrote: > Try: > > plot (myCatVariable) > > Andrew Miles > Department of Sociology > Duke University > > > On Dec 1, 2010, at 2:51 PM, phoebe kong wrote: > > Hi, >> >> Can someone tell me how to draw a histogram for th

Re: [R] default arguments and '...' in a function

2010-12-01 Thread Sarah Goslee
Antonio, You need to compare the names of list(...) with the arguments you wish to check. Here's one way to do so (note that I replaced c with d, because c is a function): f <- function(a, ...) { argnames <- names(list(...)) # check whether b is an argument if(!("b" %in% a

Re: [R] draw categorical histogram

2010-12-01 Thread Jorge Ivan Velez
Hi Phoebe, Try x <- c(12, 33, 56, 67, 15, 66) names(x) <- c('Richard','Minnie','Albert','Helen','Joe','Kingston') barplot(x, las = 1, space = 0) HTH, Jorge On Wed, Dec 1, 2010 at 2:51 PM, phoebe kong <> wrote: > Hi, > > Can someone tell me how to draw a histogram for the following summary? >

Re: [R] missing values

2010-12-01 Thread Iasonas Lamprianou
thank you, I'll have a good look and come back to you if necessary Dr. Iasonas Lamprianou Assistant Professor (Educational Research and Evaluation) Department of Education Sciences European University-Cyprus P.O. Box 22006 1516 Nicosia Cyprus Tel.: +357-22-713178 Fax: +357-22-590539 Honor

Re: [R] map() and pdf clipping

2010-12-01 Thread Paul Murrell
Hi On 1/12/2010 4:44 a.m., Peter Ehlers wrote: On 2010-11-30 04:56, Ben Tupper wrote: Hi, I thought it might help if I posted the resulting images. This is the pdf file where the map polygons are not clipped to the plotting boundary. http://dl.dropbox.com/u/8433654/test-map.pdf And this is

Re: [R] draw categorical histogram

2010-12-01 Thread Andrew Miles
Phoebe, In addition to the barplot method below, you really can use plot() to draw what it sounds like you are looking for IF you have a categorical variable. To illustrate, try running the following code: x=sample(c("Richard", "Minnie", "Albert", "Helen", "Joe", "Kingston"), 50, replace

[R] TukeyHSD with Type III SS

2010-12-01 Thread Ali S
Is there a way to run TukeyHSD with Type III instead of Type I SS? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-proje

[R] read a matrix to find a value

2010-12-01 Thread Barroso, Judit
Dear all of you, I would like if someone can help me to know how I can read a matrix looking for positions in the matrix that are not null. I have a matrix with most of the "cells" or "positions" null, I am trying to know the invasion power of a plant (number 1 in the matrix when t==1). I would

Re: [R] read a matrix to find a value

2010-12-01 Thread Phil Spector
Judit - Does this example provide any insight? pop <- matrix(0,7,7)# Population matrix when t==0 pop[sample(1:length(pop),1)] <-1 #Population matrix when t==1 which(pop!=0,arr.ind=TRUE) row col [1,] 3 4 - Phil Spector

[R] Perl "cut" equivalent in R

2010-12-01 Thread Jonathan Flowers
Does anyone know of a command in R that is equivalent to the =cut function in Perl? Thanks. Jonathan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the post

Re: [R] Perl "cut" equivalent in R

2010-12-01 Thread Sarah Goslee
It would help if you told us what you wanted this function to do, and provided an example. Not everyone speaks Perl. Sarah On Wed, Dec 1, 2010 at 6:10 PM, Jonathan Flowers wrote: > Does anyone know of a command in R that is equivalent to the =cut function > in Perl? > > Thanks. > > Jonathan > -

Re: [R] Perl "cut" equivalent in R

2010-12-01 Thread Jonathan Flowers
Thanks Sarah, you are right. In perl, "cut" serves the function of eliminating a block of code from being executed in a script. When "=cut" is placed above and below the code that you do not wish to execute then the interpreter will skip over the code. There are lots of ways to solve the proble

[R] Multivariate time series - Poisson with delayed lags

2010-12-01 Thread Ophedia
Hi all, How can a multivariate Poisson time series be modeled? Aspects of glm, forecast, dse and dynlm seem relevant but not quite complete--but hopefully what I am missing is how to assemble them effectively. What I am looking to do is model my dependent variable y_t as a Poisson family functi

Re: [R] log-normal Centile

2010-12-01 Thread Peter Ehlers
On 2010-12-01 07:09, Robert Quinn wrote: I am having problems trying to get R to graph data input that is log-normal on the horizontal (x) axis. The data is log (base 10), and I am more interested in viewing the tails of the distribution. The closest I can get with this is log on the vertical (

Re: [R] Perl "cut" equivalent in R

2010-12-01 Thread David Winsemius
On Dec 1, 2010, at 6:59 PM, Jonathan Flowers wrote: Thanks Sarah, you are right. In perl, "cut" serves the function of eliminating a block of code from being executed in a script. When "=cut" is placed above and below the code that you do not wish to execute then the interpreter will ski

Re: [R] Perl "cut" equivalent in R

2010-12-01 Thread Peter Ehlers
On 2010-12-01 15:59, Jonathan Flowers wrote: Thanks Sarah, you are right. In perl, "cut" serves the function of eliminating a block of code from being executed in a script. When "=cut" is placed above and below the code that you do not wish to execute then the interpreter will skip over the co

Re: [R] default arguments and '...' in a function

2010-12-01 Thread Gabor Grothendieck
On Wed, Dec 1, 2010 at 2:34 PM, wrote: > Dear R-users, > > I'm trying to work out a way to set default values for arguments in a > function which includes the optional argument '...'. > In practice, I have a 'plot' method for a function which specifies different > types of plots. Every differen

Re: [R] default arguments and '...' in a function

2010-12-01 Thread Bert Gunter
Look at the code for plot.default to see how it handles this issue usling locally defined versions of functions with (possibly) locally set parameter values. Cheers, Bert On Wed, Dec 1, 2010 at 6:39 PM, Gabor Grothendieck wrote: > On Wed, Dec 1, 2010 at 2:34 PM,   wrote: >> Dear R-users, >> >>

Re: [R] Problem in reading Excel spreadsheets

2010-12-01 Thread Stephen Liu
Hi Ivan, Thanks for your further advice. RODBC works without problem. Following steps are correct. > library(RODBC) > data=odbcConnectExcel(file.choose()) > sqlTables(data) TABLE_CAT TABLE_SCHEM 1 C:\\Users\\satimiswin764\\Documents\\copy_of_RDemoDe

[R] Downloading quote data from yahoo finance

2010-12-01 Thread Deb Midya
Hi R users,   Thanks in advance.   I am using R 2.12.0 on Windows XP.   May I request you to assist me in the following please.   1. I am getting error while downloading quote data from yahoo finance.   The example code is below (taken from tseries help):   library(tseries)   con <- url("http://quo

Re: [R] Perl "cut" equivalent in R

2010-12-01 Thread Petr Savicky
On Wed, Dec 01, 2010 at 07:12:37PM -0500, David Winsemius wrote: ... > One of the workarounds... the one I remember anyway... has been to > insert: > > if(FALSE) { > stuff your don't want executed > } > Switching a block of code off/on with editing a single character may be done usin

<    1   2