similar to: Spatial sampling problem

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. But, I would like to have 1.52%. Please advise me about the function. Thanks! Fabrice

Dear list, Hello! I have a question about how to print a label in the plot. I am using the following code: <pdf("mel4_chr3_11cancer_cghFLasso.pdf", height=6, width=5);plot(Disease.FL, index=i, type="Single",main="Plot of Labels");dev.off(); But "Plot of Labels" has not been printed. Any suggestions? Thanks a lot! Allen

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 should be 2 3 4 5 6 3 4 5 6 4 5 6 5 6 6 An easy way would be to do this x <-

Not really a bug (I eventually discovered), but ... When running GaussRF()/DoSimulateRF() with n>1, the default pch="!" prints backspaces to the terminal, which causes problems (LaTeX barfs) if you're running this within an Sweave() job. Setting RFparameters(pch="") fixes the problem. Could put this in the documentation, but I don't know where you could

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 Theorem). Because the location of my data is large and located near the pole, assuming that latitude and

I have asked the same question on stackoverflow but did not get a satisfying answer. I am trying to simulate a lognormal spatial random field but I need the simulated value in a certain range. So I need some easy to use functions to generate a truncated Gaussian field to start with. To be specific, I need a function like GaussRF from the RandomFields package or grf from the geoR package to

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 to observations whose covariates have high standard errors. The gam() function in the ?mgcv? package

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]]

Hi, I'm using RF to simulate a correlated variable with GaussRF set.seed=1 GaussRF(sim.kfinegrid, grid=F, model="exponential", param=c(0,0.5,0,0.2)) However when I simulate again using the same random seed I get different results. > set.seed=1 > summary(GaussRF(sim.kfinegrid, grid=F, model="exponential", param=c(0,0.5,0,0.2))) Min. 1st Qu. Median Mean 3rd

I don't understand what is happening! I have a (large) object sim1, an matrix list with dim c(101,101) where each element is an 3*3 matrix. I am subsetting that with a matrix coo, of dim c(100,2), of unique indices, but the resulting object has length 99, not 100 as expected. Code reproducing the problem follows: library(RandomFields) set.seed(123) sim0 <- GaussRF(x=seq(0, 100, by=1),

Dear list, Could someone show me where can I find the detailed formula on how to calculate the two way anova with unbalanced design? Say, if I have 2*2 design with 10,20,30,40 samples in each of the 2*2 cells. Most of the places I've googled only show how to calculate using software such as R, but not clear the detailed formula for calculating this. Thanks, Jack [[alternative HTML version

Hi. I am plotting graphs for values ranging between -1 and 10, for example: (1,2,1,1,6,7,-1,-1,5,-1) I am trying to plot the graphs so that the points with value of -1 will be in one specific color, and the rest of the points will be in one different specific color. I would be grateful for any idea of how to do that in two colors. Thanks for any hint. -- View this message in context:

Hi, Please see the email exchanges below. I am having trouble generating output that is large enough for our needs, specifically when using the GaussRF function. However, when I wrote Dr. Schlather (the author of the GaussRF function), he indicated that there is also a limit imposed by R itself. Is this something that we can overcome? Thank you very much for any assistance you may provde.

Hi, I'm working with a numeric matrix with 55columns and 581012 rows and I'm having problems in allocation of memory using some of the functions in R: for example lda, rda (library MASS), princomp(package mva) and mvnorm.etest (energy package). I've read tips to use less memory in help(read.table) and managed to use some of this functions, but haven't been able to work with

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 stored in the data.frame used for the regression analysis. For example each of my data points is made up of a field measured forest volume value and a remotely sensed vegetation estimate (NDVI). Each point is an enumeration

Hi! After unziping the CRAN file e1071.zip on my Win95 library directory, I start R and type >library(e1071) After that, help commands of the e1071 functions work (i.e., help(read.pnm), BUT some functions themselves are not available: > read.pnm Error: Object "read.pnm" not found I've looked at the file library/R/e1071 with an ascii editor and it seems like the functions

I'm interested in writing a function that constructs a new plot on the current graphics device if no plot exists there yet, but adds lines to the existing plot if a plot is already there. How can I do this? It seems to me that the exists() function might be co-opted to do this, but it's not obvious how. Many thanks, Roger -- Roger Levy Email: rlevy at ucsd.edu

Hi I am working on R 2.5.0 on window. I am trying to provide a Hyper-link to the user as a result, I have tried using gtkLinkButton to exercise the facility, however, i am not able to perform the required task, i.e. when I clicked on the LinkButton actually nothing happened. I have gone through the documentation for the same and found that GtkLinkButtonUriFunc is a function which is require to

Good Night I made different test to check normality and multinormality in my dataset, but I don´t know which test is better. To verify univariate normality I checked: shapiro.test, cvm.test, ad.test, lillie.test, sf.test or jaque.bera.test and To verify multivariate normal distribution I use mardia, mvShapiro.Test, mvsf, mshapiro.test, mvnorm.e. I have a dataset with almost 1000 data and 9

Hello. I would like to estimate the correlation coefficient from two samples with Bootstrapping using the R-function sample(). The problem is, that I have to sample pairwise. For example if I have got two time series and I draw from the first series the value from 1912 I need the value from 1912 from the second sample, too. Example: Imagine that a and b are two time series with returns for