Re: [R] confine a dataframe with a polygon

Roger Bivand wrote: On Tue, 3 Apr 2007, temiz wrote: hello how can I confine a dataframe with a known polygon that is part of whole data ? Could I suggest following up this question on the R-sig-geo list? It is not clear what you want to do - are the rows in the data frame

Re: [R] HPDinterval problem

Hi, Sorry for multiple postings. Please consider the following reproducible example: grp1 - rep(1:24, times = 24, each = 1) grp2 - rep(1:12, times = 2, each = 24) set.seed(1) out - as.integer(abs(rnorm(576))*10) x1 - out [25:552] x2 - out [49:576] mydf - data.frame(cbind(out[1:528], x1,

[R] transition matrices

I am working with transition matrices of sequences of animal behaviours. What I would like to do is parse the original matrices, adjusting row/column order so that the matrix has its main values in blocks surrounding the diagonal. This would cause behaviours involved in functional groupings

Re: [R] REPLACE and REMOVE fns

Tina Robles wrote: I want to use the replace function on R_Jeter_04_post (see attached) so that S=1, D=1, T=1, HR=1, O=0, K=0, E=0, FC=0, and W,IW,and HP are removed, so that i have a simple list of 1's and 0's. There are just 0's and 1's in your attachment ... Uwe Ligges I understand

Re: [R] time zone problems

Hi Marc, This R Help Desk article was very helpful for me: http://cran.r-project.org/doc/Rnews/Rnews_2004-1.pdf Gabor Grothendieck and Thomas Petzoldt And Gabor noted subsequent changes in behavior of POSIXt classes since the aforementioned R News publication: (1)

Re: [R] Testing additive nonparametric model

I'd use mgcv::gam to fit a 2D smooth, rather than loess, and then do an *approximate* likelihood ratio/ F test. For example library(mgcv) ## fake some data set.seed(0) n-400;sig2-4 x0 - runif(n, 0, 1);x2 - runif(n, 0, 1) f - 2 * sin(pi * x0) f - f + 0.2 * x2^11 * (10 * (1 - x2))^6 + 10 *

[R] Problem with R library

When I installed RSRuby (Ruby interface to the R programming language), an error occured : cannot find the R library. I defined Rlibrary : R\library in PATH on Windows but it didn't work with the same error. RSRuby ask R to permit the option 'shared library' but when I install R, it didn't ask me

[R] Need mtrace (malloc debugging) help

Hi all, I am using mtrace malloc debugging utility to trace memory leaks in my application. I followed following steps: 1. Added mtrace() function in the beginning of my main function. 2. exported the MALLOC_TRACE environment variable. $ export MALLOC_TRACE=memory.log 3. Run the program 4.

[R] Silhouette width using K-means

Hello, I am doing clustering and I want to know how can i find Silhouette width using K-means. Just like PAM (code below). (2) Secondly, I have mixed data all sort of variables numeric, categorical, nominals so I first change to all nominal to binary and normlise the data before any clustering.

Re: [R] Problem with R library

bisous wrote: When I installed RSRuby (Ruby interface to the R programming language), an error occured : cannot find the R library. I defined Rlibrary : R\library in PATH on Windows but it didn't work with the same error. RSRuby ask R to permit the option 'shared library' but when I install

Re: [R] gnomeGUI under KDE

I use the Kate editor (FC6, KDE, R 2.4.1) whihc allows me to pipe commands to a console window with a key binding. Works well. Richard On Wed, 2007-04-04 at 00:19 -0400, Luis Lopez Oliveros wrote: Dear all, Recently, I installed the KDE desktop under my Linux system, but Included also the

[R] Info on SPATSTAT window and maps

I wonder if there is any tutorial explaining, step by step, how to convert a (georeferenced) map boundary (from esri shape-file) into a Spatstat window, for performing the analysis of marked point patterns surveyed inside that map. Any help on the subject will be really appreciated. Giuseppe

[R] Annotate a levelplot (using abline) - Difficulty with trellis.

Hi, I am generating a beautiful plot with the 'levelplot' function over my square matrix of data. In order to help visualise the data I would like to draw a diagonal line on the matrix. Because the plot is actually a trellis object, I am having difficulty working out how to do this. I have been

Re: [R] Info on SPATSTAT window and maps

On Wed, 4 Apr 2007, Giuseppe Brundu wrote: I wonder if there is any tutorial explaining, step by step, how to convert a (georeferenced) map boundary (from esri shape-file) into a Spatstat window, for performing the analysis of marked point patterns surveyed inside that map. Any help on the

Re: [R] Need mtrace (malloc debugging) help

On Wed, 4 Apr 2007, $|-|@r@|) wrote: Hi all, I am using mtrace malloc debugging utility to trace memory leaks in my application. I followed following steps: 1. Added mtrace() function in the beginning of my main function. 2. exported the MALLOC_TRACE environment variable. $ export

Re: [R] Reading user input

That is exactly what the code does. On 4/3/07, projection83 [EMAIL PROTECTED] wrote: I dont know if thats what i want (its not working like i want). I would like to have R pause, wait for a user input of a number, then store that number as a variable to use in the next calculation it

Re: [R] Newbie: Simple loops: complex troubles

One of the things about R that you have to learn is vector operations. You try to avoid loops and also 'generating' variables -- this is where a list comes in very handy. To relate to Java, it is similar to 'struct'. Here is a program that does what you want to do; it uses lists and vectorized

[R] to findout maximized log likelihoods by using rlarg.fit (for several r order statistics)

Dear R helpers, I need to find out maximized log likelihoods, parameters estimates and standard errors (in parentheses) of r largest-order statistics model, with different values of r by using the function rlarg.fit. I want to specify required number of order statistics to the model. I attached

Re: [R] Annotate a levelplot (using abline) - Difficulty with trellis.

On Apr 4, 2007, at 7:19 AM, Dan Bolser wrote: Hi, I am generating a beautiful plot with the 'levelplot' function over my square matrix of data. In order to help visualise the data I would like to draw a diagonal line on the matrix. Because the plot is actually a trellis object, I am having

Re: [R] Reading user input

I think the problem the OP may be having is that the code will not work if you put it in a document window in R and then tell R to source the document. At least not with R.app in MacOSX. This is what happened when I did it: source(/tmp/Rtmp0TQfA6/file10d63af1) enter the number of groups:

Re: [R] Info on SPATSTAT window and maps

And what does this have to do with the lead thread subject line Transition Matrices? I guess the only way we can start enforcing thread discipline is to stop responding to threads that hijack others? Any thoughts? Ranjan On Wed, 4 Apr 2007 11:31:33 +0200 Giuseppe Brundu [EMAIL PROTECTED]

Re: [R] transition matrices

Hi, It appears that you are trying to separate the states in the transition matrix such that you have recurrent classes all next to each other. Here is an idea off the top of my head: make equivalence classes from your statespace and then use that to create a block diagonal matrix. I did a

Re: [R] the numimum number of fixed factors lme package can deal with

On 4/3/07, Fang, Yongxiang [EMAIL PROTECTED] wrote: Dear Douglas, Thanks for your help. the error message is: 'Error in eval(expr, envir, enclos) : Object dx.200 not found'; In fact, dx is the design matrix and now in a data frame: dframe I have checked if dx.200 in dframe. dframe$dx.200

Re: [R] the numimum number of fixed factors lme package can dealwith

As Doug noted yesterday, you may have a limit on the memory needed to create the model matrix for the fixed effects. But, aside from that, based on what I see below, it appears you have manually created a model matrix yourself. If that's true, you don't need to do that. You can create a single

[R] GAM/GLM parameters

I think this might be a very basic question, but is there a simple way to characterise the relationships that a gam or lm model have identified? I am trying the create species distribution models based on climate, and want to know whether, for example, higher temperatures (one of the predictor

Re: [R] REPLACE and REMOVE fns

See attached...sorry On 4/4/07, Uwe Ligges [EMAIL PROTECTED] wrote: Tina Robles wrote: I want to use the replace function on R_Jeter_04_post (see attached) so that S=1, D=1, T=1, HR=1, O=0, K=0, E=0, FC=0, and W,IW,and HP are removed, so that i have a simple list of 1's and 0's. There

[R] selection of character

Hi all, I have a simple question: I have a character vector like this [1] 10_1 1_1 11_1 2_1 3_1 4_1 5_1 6_1 7_1 8_1 [11] 9_1 1_3 2_3 3_3 4_3 1_4 2_4 3_4 4_4 5_4 [21] 6_4 7_4 1_5 2_5 10_6 11_6 12_6 1_6 2_6 3_6 [31] 4_6 5_6 6_6 7_6 8_6 9_6 1_7 1_8 2_8 3_8 [41] 10_9

Re: [R] selection of character

one way is the following: x - c(10_1, 1_1, 11_1, 2_1, 3_1, 4_1, 5_1, 6_1, 7_1, 8_1) # sapply(strsplit(x, _), [, 2) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address:

[R] Problems loading package GeoXp

Dear useR's, I wanted to load the package GeoXP in my computer but I couldn't because this warning message appeared: require(GeoXp) Carregando pacotes exigidos: GeoXp Carregando pacotes exigidos: tcltk Loading Tcl/Tk interface ... Erro em fun(...) : couldn't connect to display :0 Além disso:

[R] Power analysis and mixed model

Hi Is there any way to compute power analysis for a mixed model ? Julien Julien Martin Candidat au Doctorat Laboratoire de Marco festa-Bianchet Université de Sherbrooke Sherbrooke, Québec J1K 2R1 Canada 819 821 8000 poste 62059 -- Checked by AVG Free Edition. 12:48

[R] pattern search

Dear R-experts, I'm looking for an easy possibility for pattern search. I have got a string and a special pattern. I would like to know if the pattern appears in the string, if yes where does it appear (position) and the number of appearance. Example: Text = c(If the sun shines, no clouds

Re: [R] Annotate a levelplot (using abline) - Difficulty with trellis.

My question was thus; Given library(lattice) my.m - matrix(seq(1,100,1),nrow=10) levelplot(my.m) How can I add a diagonal line onto the resulting 'color square'? The answer I found was to hack the 'panel.levelplot' function. Here is the diff between the old (panel.levplot) and the new

Re: [R] Power analysis and mixed model

One option is to use lmer in a monte carlo simulation. I just did this last week. Check out the article published in the American Statistician and can be found at http://maven.smith.edu/~nhorton/R/r.pdf. The article is not about power per se, but is about R as a toolbox for mathematical

[R] argmax

Hello, Is there any function that returns the argmax of a vector ? For example I have the vector x=(1,4,15,6,7) max(x)=15 I want a function f such that f(x)=3 because this is where max(x) is. Thank your very much. P.S. : I know it is a basic question but I have not figured it out by

Re: [R] Annotate a levelplot (using abline) - Difficulty with trellis.

Dan Bolser said the following on 4/4/2007 7:52 AM: My question was thus; Given library(lattice) my.m - matrix(seq(1,100,1),nrow=10) levelplot(my.m) How can I add a diagonal line onto the resulting 'color square'? The answer I found was to hack the 'panel.levelplot'

Re: [R] gnomeGUI under KDE

On 4/3/07, Luis Lopez Oliveros [EMAIL PROTECTED] wrote: Dear all, Recently, I installed the KDE desktop under my Linux system, but Included also the Gnome desktop (actually, I installed everything that comes in the openSUSE cd). My question about R is the following. I would like to run

Re: [R] argmax

x - c(1,4,15,6,7) which(x==max(x)) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Wednesday, April 04, 2007 11:14 AM To: r-help@stat.math.ethz.ch Subject: [R] argmax Hello, Is there any function that returns the

Re: [R] argmax

check ?which.max(), e.g., x - c(1,4,15,6,7) which.max(x) 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:

Re: [R] Power analysis and mixed model

-Original Message- From: Jean G. Orelien [mailto:[EMAIL PROTECTED] Sent: Wednesday, April 04, 2007 11:22 AM To: 'Julien' Subject: RE: [R] Power analysis and mixed model Julien, You could write your Linear Mixed Model as a General Linear Multivariate Model (GLMM) and use the power

Re: [R] pattern search

how about: length(gregexpr(the, Text)[[1]]) ? b On Apr 4, 2007, at 10:49 AM, Schmitt, Corinna wrote: Dear R-experts, I'm looking for an easy possibility for pattern search. I have got a string and a special pattern. I would like to know if the pattern appears in the string, if yes where

Re: [R] pattern search

Hi, I'm looking for an easy possibility for pattern search. I have got a string and a special pattern. I would like to know if the pattern appears in the string, if yes where does it appear (position) and the number of appearance. Example: Text = c(If the sun shines, no clouds should

Re: [R] argmax

which.max(x) Il giorno mer, 04/04/2007 alle 15.13 +, [EMAIL PROTECTED] ha scritto: Hello, Is there any function that returns the argmax of a vector ? For example I have the vector x=(1,4,15,6,7) max(x)=15 I want a function f such that f(x)=3 because this is where max(x) is.

Re: [R] transition matrices

Hi Richard, You might want to have a look at Bertin's reorderable matrices, which seem to be solving a similar problem. Googling for bertin reorderable matrix gives a few automatic implementations. Hadley On 4/4/07, Richard Rowe [EMAIL PROTECTED] wrote: I am working with transition matrices

[R] gpclib library

Hi all, I'm tring to download the gpclib library but I don't find it in the old repository install.packages(c(gpclib), dependencies=TRUE, repos=http://cran.cnr.berkeley.edu/;) Avviso in install.packages(c(gpclib), dependencies = TRUE, repos = http://cran.cnr.berkeley.edu/;) : argument

Re: [R] transition matrices

I think what you are looking for is also called seriation or ordination, i.e. find the best enumeration order of a set of objects. Take a look at the following link for more information and on criteria for defining best http://www.lirmm.fr/~caraux/PermutMatrix/EN/Seriation.htm The link is for a

[R] Cannot install iWidgetsRGtk

Dear All I am trying to install iWidgetsRGtk (on a Linux machine), but getting the following error: install.packages(iWidgetsRGtk,repos=http://www.math.csi.cuny.edu/pmg;) Warning in download.packages(unique(pkgs), destdir = tmpd, available = available, : no package 'iWidgetsRGtk' at

[R] (sans objet)

Hello, Is there any function argmax in R ? For example, if I have the vector x=c(1,4,15,7,6), max(x)=15. But I need a function f such that f(x)=3, 3 being the number where 15 is. Thank you very much. _ ues clics pour

Re: [R] selection of character

Another way, with regexps : gsub(.*_, , x) HTH, -- Julien __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented,

Re: [R] gpclib library

It's a _package_ that you are looking for (there is something else called a gpclib library, to do with GNU Pascal). You haven't told us your OS or R version, so we don't know what file you are looking for. But try another mirror if one appears not to work (and using a US mirror with an email

Re: [R] Cannot install iWidgetsRGtk

Hi, iWidgetsRGtk is an obsolete package. Instead, install gWidgetsRGtk from CRAN. That's g not i. Michael On 4/4/07, Paul Smith [EMAIL PROTECTED] wrote: Dear All I am trying to install iWidgetsRGtk (on a Linux machine), but getting the following error:

[R] Finding a single unique item in duplicated vectors

I have a very simple problem and am completely missing the solution. I have two character variables (character ID's from two datasets) Data set 'b' is data set 'a' with one more subject added. How do I find out which is the added subject? I have tried duplicated and unique without much success.

Re: [R] Cannot install iWidgetsRGtk

On 4/4/07, Michael Lawrence [EMAIL PROTECTED] wrote: iWidgetsRGtk is an obsolete package. Instead, install gWidgetsRGtk from CRAN. That's g not i. Thanks, Michael. That is solved now. Paul On 4/4/07, Paul Smith [EMAIL PROTECTED] wrote: Dear All I am trying to install iWidgetsRGtk

Re: [R] Finding a single unique item in duplicated vectors

On Wed, 4 Apr 2007, John Kane wrote: I have a very simple problem and am completely missing the solution. I have two character variables (character ID's from two datasets) Data set 'b' is data set 'a' with one more subject added. How do I find out which is the added subject?

Re: [R] Finding a single unique item in duplicated vectors

--- John Kane [EMAIL PROTECTED] wrote: I have a very simple problem and am completely missing the solution. I have two character variables (character ID's from two datasets) Data set 'b' is data set 'a' with one more subject added. How do I find out which is the added subject? I have

[R] rgl.postscript generating pdf: floating xlab, font size, transparency

Hello, I very much enjoy rgl! But when converting the screen output to pdf using rgl.postscript I have encountered three problems: (1) I don't seem to have control over font sizes and point sizes (of added points) in the pdf output. I tried cex from plot but that does not have an effect.

Re: [R] Finding a single unique item in duplicated vectors

Arrgh! I knew it going to be was blindingly obvious once someone showed it to me. I had even looked at diff() at one point but never found setdiff. ?sets did not seem to find anything. Thanks very much. --- Charles C. Berry [EMAIL PROTECTED] wrote: On Wed, 4 Apr 2007, John Kane wrote: I

Re: [R] (sans objet)

On 4/4/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Hello, Is there any function argmax in R ? For example, if I have the vector x=c(1,4,15,7,6), max(x)=15. But I need a function f such that f(x)=3, 3 being the number where 15 is. which.max __

Re: [R] transition matrices

On Wed, 4 Apr 2007 07:46:37 -0500 Ranjan Maitra [EMAIL PROTECTED] wrote: Hi, It appears that you are trying to separate the states in the transition matrix such that you have recurrent classes all next to each other. Sorry, I meant non-communicating classes. Btw, thinking about this some

[R] R and bash script problem

Hi, I have a problem with R package and bash. I have the package scatterplot3d installed and working. I'm using bash as the command line interpretor. When I open R interactively I can load scatterplot3d without problem and use it. When I create an R script it also works and

Re: [R] Annotate a levelplot (using abline) - Difficulty with trellis.

On 04/04/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Dan Bolser said the following on 4/4/2007 7:52 AM: My question was thus; Given library(lattice) my.m - matrix(seq(1,100,1),nrow=10) levelplot(my.m) How can I add a diagonal line onto the resulting 'color square'?

Re: [R] rgl.postscript generating pdf: floating xlab, font size, transparency

On 4/4/2007 1:17 PM, Kim Milferstedt wrote: Hello, I very much enjoy rgl! But when converting the screen output to pdf using rgl.postscript I have encountered three problems: (1) I don't seem to have control over font sizes and point sizes (of added points) in the pdf output. I tried

[R] compile errors with rgl-0.70.570 on FreeBSD

Thank you for your work on rgl. Reading in another thread about your new (inofficial) version of rgl (see below) I tried it with R version 2.5.0 alpha (2007-03-31 r40986) under FreeBSD 7.0-CURRENT. Unfortunately I got the following error: --- #R CMD INSTALL rgl_0.70.570.tar.gz * Installing

[R] Power analysis and mixed model

Like Harold, I have also used simulation-based methods for power analyses, which I discuss (briefly) in the following article (with R code in an appendix): Atkins, D. C. (2005). Using multilevel models to analyze marital and family treatment data: Basic and advanced issues. Journal of

[R] rcmdr help

Dear All, I use Debian etch and downloaded R from its repo. Although it is running I would like to make the first steps using rcmdr. I went through on wikis and howtos, installed rcmdr but it is not running. Apparently it is properly installed, there is an Rcmdr library in

Re: [R] REPLACE and REMOVE fns

Does this do it? x - readLines(/temp/Giambi_03_post.txt) x.in - c(S, D, T, HR, O, K, E, FL) x.out - c( 1, 1, 1, 1,0, 0, 0, 0) # map the characters. Assume that missing will be deleted x - x.out[match(x, x.in)] x - x[!is.na(x)] x [1] 0 0 0 0 1 0 0 1 0 0 0 0 0 1 1 0 0 0 0

[R] unexpected behavior when creating a list of functions

I wanted to create a list of functions whose output differs depending the value of a variable when the function was created. Generally this did not work. Each function was exactly the same, as in the simple example below: get_data_function - function(v) { function() { print(v) } }

[R] transfer function models

Dear Sam, I´ve found your posting on https://stat.ethz.ch/pipermail/r-help/2005-January/063966.html about transfer funtion in SAS. I´m writing my diploma thesis in Germany right now and can´t find any good documentation which would explain how dynamic regression is realized in SAS.Do you

Re: [R] unexpected behavior when creating a list of functions

'v' was not evaluated when you defined the function; this is 'lazy' evaluation in R; try get_data_function - function(v) { v # cause 'v' to be evaluated function() { print(v) # now it is defined when the function is } } data_functions - sapply(1:10,function(v) get_data_function(v))

Re: [R] unexpected behavior when creating a list of functions

On Wed, 4 Apr 2007, Matthew Suderman wrote: I wanted to create a list of functions whose output differs depending the value of a variable when the function was created. Generally this did not work. Each function was exactly the same, as in the simple example below: get_data_function -

[R] how to image.plot a XY grid file into a lat-lon map

Hi All, I have a netcdf gridded file with LCC projection. I can easily use image.plot to visualize it. However, as the axises are in X,Y, not Lat and Lon, I could not add state or country maps onto it (or lat lon information). I do have a grid2d file that describes the lat and lon for each (X,Y)

[R] Moving plot labels

Hi. I'm sure this is a simple problem, yet I can't seem to get simple help for it. I am simply trying to move my xlab in plot(). I am currently using the following commands: plot(c(-0.25,18),c(0, max(patient10)),type=n, ylab=SD of POST estimator, xlab=Scans \n (a)) But when the plot prints,

Re: [R] rcmdr help

Dear Daniel, I'm not a Debian user, but your description of what you did makes me wonder whether (1) you ever loaded the Rcmdr package, and (2) whether you installed its dependencies. After starting up R, you can load (and start) the R commander via the command library(Rcmdr). Prior to that, you

Re: [R] Newbie: Simple loops: complex troubles

This helped a lot. Thank you so much - im now able to get some basic stuff moving around in R! Petr Klasterecky wrote: Not sure whether this is exactly and everything you want, but at least it may give you some ideas how to proceed. You do not need loops at all: Let's try a simplified

[R] Problem installing rgdal on Mac OS X

Please forgive me for posting this here if it is not appropriate (where should this question be posed?)… I am attempting to install rgdal (from source, since there is no binary) on Mac OS X 10.4.9 on a PowerPC G3 Macintosh. I have R version 2.4.1 (2006-12-18) installed. Following the install

Re: [R] Moving plot labels

Dean Try: plot(c(-0.25,18),c(0, max(patient10)),type=n, ylab=SD of POST estimator, xlab=) title(xlab='Scans\n(a)', line=3) Peter Alspach -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Thursday, 5 April 2007 12:21 p.m.

[R] Package for common rater agreement statistics?

Could someone point me to a package that generates common agreement statistics for classified observations (particularly the F1 measure)? It seems like there would have to be a package in cran that could handle this, but I cant seem to find it. The irr package does alot of interrater reliability

[R] changing existing graphs optimally

hi there, right now i am going about things a round about route with graphs in that if i want to change the scale to log, I have to to a reset with add=F, then program things to be redrawn from scratch - bad coding basically:P in plot(), the first arguement it needs is data - how do you

Re: [R] Problem installing rgdal on Mac OS X

On Wed, 4 Apr 2007, Thomas Adams wrote: Please forgive me for posting this here if it is not appropriate (where should this question be posed?)? The R posting guide asked you to ask the maintainer, but as this is a MacOS X issue, on the R-sig-mac list. It's not appropriate for R-help

[R] TM convert Newsgroup to plain text

Hi, Anyone tried to do stemming for emails (mbox format) using the Text Mining (tm) package in R. According to the documentation the text document collection has to be converted to text and there is no method to convert Newsgroup (email) to plain text. Thanks, Saeed