similar to: thicker axis in levelplot

Dear R-Help community, I am trying to overlay a single contour line over a correlation plot using levelplot in the lattice package. These are the two arrays: 1) a correlation plot over Africa - so each grid square is a different colour dependent on correlation - this is in an array: result_cor with dim[465,465] 2) a single contour line from a ***different data source*** - this is from data

Dear R-community, I am looking for some simple advice - I have a matrix (therefore 2 dimensional) of global temperature. Having read R-help I think that when I ask R to image() or levelplot() my matrix will it actually appear upside down - I think I therefore need to use the line: > levelplot(temperature.matrix[,ncol(output.temp):1], ........) to get it looking like it was on the globe due

Dear R-help community, Does anybody know of a stats function in R that tells you which distribution best fits your data? I have tried look through the archives but have only found functions that tell you if it's normal or log etc. specifically - I am looking for a function that tells you (given a timeseries) what the distribution is. Any help/advice will be greatly appreciated, All the

Dear R-help, Which random number generator function would you recommend for simply picking 15 random numbers from the sequence 0-42? I want to use replacement (so that the same number could potentially be picked more than once). I have read the R-help archives and the statistics and computing book on modern Applied statistics with S but the advice seems to be for much form complicated

Dear R-help community, I have looked on the R search site and archives but cannot find mention of a way of interpolating a gaussian distribution of data to a standard 2.5 degree grid. I have two global dataset and I need to correlate - unfortunately one is a 2.5 degree grid dim[longitude=144,latitude=72] and one is gaussian dim[longitude=192,latitude=94]. I would rally appreciate hearing

Dear R-help, I have a loop, which is set to take about 26 hours to run at the rate it's going - this is ridiculous and I really need your help to find a more efficient way of loading up my array gpcc.array: #My data is stored in a table format with all the data in one long column #running though every longitute, for every latitude, for every year. The #original data is sotred as

Hi Jenny If you want a general solution I understand. However I just downloaded the file fine (as far as I can tell) so you are welcome to have a copy. I can email it to you if you want. I do not think your test for NA is valid. i.e if(test != "NA"){ } I think you should use if(is.na(test)){ } Or something similar. J --- John Seers Institute of Food Research Norwich

Dear R-Help, I am trying to find a function that will give me the significance of the correlation of 2 variables (in the same dimension arrays) correcting for serial autocorrelation. How can I view the function cor.test's code? I would like to know a lot more detail about the function than written in the documentation at http://finzi.psych.upenn.edu/R/library/stats/html/cor.test.html to

Dear R-help community I am trying to write an R object (data.out) to a file in order to re-access it later and not have to re-load up the array with data every time. Here is the form of data.out > data.out <- list(lats=seq(88.542, -88.542, length=94), lons=seq(0, 360-1.875, length=192), date=vector(length=nyr*12), data=array(NA, c(nyr*12, 94*192)) ) I tried >

Dear R-Help friends, I am unable to get the latest version of R (2.4.1) to compile on my solaris 10 system - has anybody else experienced this problem and are you able to offer me any advice? I appreciate your time, many thanks, Jenny Barnes Here are my CURRENT specifications: platform sparc-sun-solaris2.10 arch sparc os solaris2.10

Dear R-help community, My main goal of this message is to find a way of skipping a file of a month/year in a loop that does not exist (and making it's output into an data.out array would be NA) and moving onto the next year/month in the loop to carry on filling data.out with real precipitation data. The situation so far: I downloaded 50 years worth of GRIB data files from the NCEP data

R does a great job with the fine details regarding plots. e.g in the following: library(lattice) y <- -4:4/10 xyplot(y~1, las=1) the y axis is labelled with numbers -0.4, -0.2, 0.0, 0.2, 0.4 with the numbers aligned on the decimal point. How do I get the same behaviour in the colorkey of a levelplot? e.g. levelplot(matrix(y,3,3)) the numbers in the colorkey seem left-aligned, and

I'm a new user of levelplot, and are not familiar with the terminology very well. Is the bar alongside the levelplot indicating color or shade call "colorkey"? I have to adjust the size of its label, but couldn't make it so far. This is how I did: levelplot(z~x*y, grid, at=seq(0,1,by=0.1), colorkey=list(labels=list(cex=2))) But I got error message "Error in

I'm using the {lattice} "levelplot" function to make a (more or less) 2-d histogram, and for the most part it's working fine with my data. However, I can't get the color key to do what I need. I can give it labels and custom cutoffs, but my cutoff lines (and hence my labels) aren't evenly spaced, instead they're more-or-less logarithmic, starting at

Dear R-help community, I have 2 different arrays of precipitation data each of the same dimensions of [longitude, latitude, time] dim=[30,32,43], called array1 and array2. I need to correlate them. This is the code I used to get one overall correlation value for the whole of the area of interest: > result <- cor(array1,array2,use="complete.obs") > result This give me a

Hi, I would like to add contour lines to a (trellis/lattice-) levelplot. Sure, there is the "contour=TRUE" argument, but this uses "cuts=..." (which is usually chosen very high for my plots. I guess cuts=99 is the best you can do (?)) for plotting the contour lines. Furthermore, I do not like the numbering of the contour lines this way. Therefore, I tried to add a

I have two questions; Do you know how to take out axes in a levelplot? Me doesn't work axes = FALSE And I would like to fix the values range of colorkey or legend from my rainbow col.regions, dou you know how can I fix this values independent of values database? Thanks and happy new year! [[alternative HTML version deleted]]

Dear R gurus, I have got stuck on how to customise the colorkey generated by levelplot and contourplot, in the lattice package. This best illustrated by an example: library(lattice) levelplot(volcano^20/1e45,at=c(0,0.001,0.01,0.1,1,10)) The reason for raising the volcano dataset to the 20th power is to create a dataset with a large vertical exaggeration - this is similar in nature to the data

Dear R'ers, I've scanned available documentation and the web, but I can't seem to figure out where I've gone wrong in adding numbers to scale the colorkey in levelplot (Lattice package). For example, levelplot(z, contour=T, labels=T, cuts=10, region=T,... colorkey=list(T, space="bottom", at=seq(0,2000,length=20),

Hi, I'm a novice with levelplot and need some assistance! Basically, I want a window which contains 6 levelplots of equal size presented in 3 columns and 2 rows. I've tried to approach it two ways. The first way leads to this question: Is there any way to concatenate levelplots from a factor vertically as opposed to horizontally? I'd like to pair the levelplots by factor.2 on