Re: [R] How to measure mode (central tendancy)

2007-08-30 Thread Julian Burgos
mode()... of course! Raymond Balise wrote: What is the name of the function to give me the mode (central tendancy) of a numeric variable that can be negative? [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] Calculating diameters of cirkels in a picture.

2007-08-23 Thread Julian Burgos
Hi Bart, If you only have 36 circles, the fastest way would be to use some image processing software and measure the circles by hand. One option is to use ImageJ, which you can download here http://rsb.info.nih.gov/ij/ Julian Bart Joosen wrote: Hi, Maybe this is more a programming

[R] weights in GAMs (package mgcv)

2007-08-14 Thread Julian Burgos
Dear list, I’m using the ‘mgcv’ package to fit some GAMs. Some of my covariates are derived quantities and have an associated standard error, so I would like to incorporate this uncertainty into the GAM estimation process. Ideally, during the estimation process less importance would be given

Re: [R] labelling plots with ancillary data in data.frame

2007-08-14 Thread Julian Burgos
Hi Wesley, You can use the text() function to add text to an existing plot. See ?text. Julian Wesley Roberts wrote: Hi All, I am busy using R to do some regression modelling and have been using plot(x,y,) to visualise my variables. I would now like to label my points using data

Re: [R] Spatial sampling problem

2007-08-07 Thread Julian Burgos
Hello SK, I suggest you try the RandomFields package. It has specific functions for simulating random fields with an array of different algorithms and distributions. The GaussRF function in particular normal generates random fields. Julian Julian M. Burgos Fisheries Acoustics Research Lab

Re: [R] SOS

2007-07-20 Thread Julian Burgos
Multiply by 100? Add R=R*100 Fabrice McShort wrote: Dear all, I am a new user of R. I would like to know how to get fund's returns in percentage (%). For example, I use: R - ts(read.xls(FundData), frequency = 12, start = c(1996, 1)) Whith this program, the returns are like 0.0152699.

Re: [R] Graphical Parameters

2007-07-20 Thread Julian Burgos
Hi Amna, To have more control on the grid of your plot use the grid function. Remove the tck argument from your plot call and add, as following: plot(Year,Dir,ylab=Annual Maximum Daily Rainfall (mm),cex.lab=1.3,type=o,lab=c(20,20,25),mgp=c(2.5,1,0),asp=1) grid() See ?grid to see how to

Re: [R] R and Copula

2007-07-19 Thread Julian Burgos
Hi Meeryana, It seems you have not loaded the package. To use a package that you already installed, do: library(copula) I recommend you review the R documentation. There are several good references and tutorial on the CRAN site: http://cran.r-project.org/manuals.html

[R] gamm (package mgcv) with large datasets

2007-07-18 Thread Julian Burgos
Dear list, I am interested in fitting a Generalized Additive Mixed Model with spatially correlated errors to a large, spatially indexed, data set (~4000 observations). My initial analysis was a Generalized Additive Model that included a two dimensional smooth term to model spatially

Re: [R] automatically generating and accessing data frames of varying dimensions

2007-07-12 Thread Julian Burgos
Replace assign(nam, 1:auto.length[i]) with assign(nam,data.frame(matrix(1:auto.length[i], ncol=3))) Verify if the matrix (and hence the data frame) is filled the way you intended. Otherwise, use the argument byrow=T. Julian M. Burgos Fisheries Acoustics Research Lab School of Aquatic and

[R] mgcv: lowest estimated degrees of freedom

2007-06-21 Thread Julian Burgos
is if the estimated degrees of freedom (EDF) for the term are close to their lower limit. What would be the minimum number of EDF's for a) Univariate thin plate regression splines(TPRS) with shrinkage, i.e. s(...,bs=ts) b) Bivariate tensor products of TPRS with shrinkage? Thanks for any help, Julian

[R] Recursive decreasing sequences

2006-10-20 Thread Julian Burgos
Hello fellow R's, I'm sure there must be an easy way to do this. But after digging in the documentation and thinking about it for a while I couldn't figure it out. I need to get a decreasing recursive vector in. I mean something like this: if starting at 2, and ending at 6, the vector

[R] Running executable files from R

2006-06-23 Thread Julian Burgos
Hello fellow R's, I apologize if this question was answer elsewhere. I have an executable file that I need to run from R. Basically I want to use R to create the input files this executable requires, then run it, and finally use R again to analyze the output files. Because I have to do this

[R] Variograms and large distances

2005-11-08 Thread Julian Burgos
Hello R list, I need to compute empirical variograms using data from a large geographic area (~10^6 km2). Although I could not find a specific reference, I assume that both geoR and gstat calculate distances among data points assuming points are on a flat surface (using the Pythagorean