similar to: Creating a map .

summary: I can display a lon-lat map on a lattice::layerplot, and I can display a Lambert conformal conic (LCC) map on a spam::image, but I can't display an LCC map on a lattice::layerplot. Example follows. What am I doing wrong? details: I've been using `lattice` (via `rasterVis`) successfully to display global atmospheric data, which works well enough (though I am definitely intrigued

I have nc files for global soil moisture,here is one file https://echange-fichiers.inra.fr/get?k=f9DDllPKdUKs5ZNQwfq from the metadata ,the projection is cylindrical and the resolution is 25 km(it is based on authalic sphere based on International 1924 ellipsoid).As I want to compare with other data, I have to make them identical. - my other data are in WGS84 with 0.25*0.25 degree resolution

Hola que tal alguien de ustedes me podría auxiliar, lo que pasa es que estoy tratando de hacer un programa en R utilizando la siguiente tabla: x<1 km x>1km Estabilidad a c d f c d f A 213 440.8 1.94 9.27 459.70 2.094 -9.6 B 156 106.6 1.15 3.30 108.20 1.098 2.0 C 104 61.0 0.91 0.00 61.00 0.911 0.0 D 68 33.2 0.73 -1.70 44.50 0.516 -13.0 E 51 22.8 0.68 -1.30 55.40 0.305

Dear all, I am using ggplot with geom_tile to print as an image a matrix  I have. My matrix is a squared one of 512*512 cells.  The code that does that is written below > print(v + geom_tile(aes(fill=dB))+ opts(axis.text.x=theme_text(size=20),axis.text.y=theme_text(size=20), axis.title.x=theme_text(size=25) , axis.title.y=theme_text(size=25), legend.title=theme_text(size=25,hjust=-0.4) ,

Dear all, I am writing as I would really need your help on the problem with gamm4. I have tried to find a solution online but I wasn't very successful. I am running a gamm4 model with an interaction between two variable using the tensor term, t2. I have a group variable (super end group) with six factors; I would like to run the model to see the how the interaction term varies across the

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 meters and not correlated beyond that. I have been having problems

Dear all, I have asked yesterday of how I can read a simple bitmap file in R cran. I was suggest to use either readGDAL or raster for loading my bitmap a. I have done it with readGDAL like     store<-readGDAL(fname='lena256.bmp')     and it works,... but it converts my matrix-like notion of a bitmap to a large vector b. Raster also returns a class that I can not understand. So 1,

Hi Maria If you have problems just start with a small model with predictions and then plot with xyplot the same applies to xyplot Try library(gamm4) spring <- dget(file = "G:/1/example.txt") str(spring) 'data.frame': 11744 obs. of 11 variables: $ WATERBODY_ID : Factor w/ 1994 levels "GB102021072830",..: 1 1 2 2 2 3 3 3 4 4 ... $ SITE_ID

Hello everyone, Could you please help me find out if R supports matrixes inside matrixes? This is what I would like to do I have an area map of humidity per km. I would like at every cell to keep also information about the height of this area, the current temperature etc. Is something like that supported? I would like to thank you in advance for your help Best Regards Alex

SUMMARY: Specific problem: I'm regridding biomass-burning emissions from a global/unprojected inventory to a regional projection (LCC over North America). I need to have boundaries for Canada, Mexico, and US (including US states), but also Caribbean and Atlantic nations (notably the Bahamas). I would also like to add Canadian provinces and Mexican states. How to put these together? General

Hi everyone, I wonder if anyone knows the reason why the outputs of the same reprojection in r and arcgis are different?. The magnitude of the change can be up to 40 km in the poles. Basically, I have a database of points equally separated by one degree over the globe. In ARCGIS, I am projecting the data in GCS-WGS-1984 and then reprojected it to Berhmann to ensure equal area

hi all, We are using squid 25s1 with kernel 2.4.19 and iproute2(+julian''s Pathes) with the following acl.. acl short_path dst 128.0.0.0/8 tcp_outgoing_address myIp2nd short_path we are linked to 2 isp--one having satelite & the other OFC. We want the above mentioned network to go thru OFC(ispB) as the sibling resides there. But when i use squidclient mgr:server_list

hi to all, frustated for this error, to day i buy a 1 GB memory slot for my laptop now it have 1,28GB instead the old 512, but i've the same error :-( damn!damn!....how can i do? repeat for a little area (about 20X20 km and res=20m) it work fine! have you any suggestion? is ther a method for look if this error depend from my ram or other....? thanks foe any suggestion! i need your help.

Hello everyone, This is Kangmin. I am trying to produce shape and scale of my wind data. My data is based on wind speed frequency with 1km/hr increment. data is described below. Windspeed (km/h) Frequency 1 351 2 147 3 317 4 378 5 527 6 667 7 865 8 970 9 987 10 907 11 905 12 642 13 1000 14 983 15 847 16 842 17 757 18 698 19 632 20 626 21 599 22 529 23 325 24 391

Hi, I am trying to do a statistical analysis of a SPSS data set, which has 7 different parameters, 4 of which are Yes/No. When I try to select in one of the parameters all the rows with a Yes it always gives me the same error : object 'Yes' not found although: > suburban [1] No No Yes No No Yes No Yes No No No No No No No No No No [19] No No Yes No Yes Yes

I run thunderbird (both 3.x and 5.x) against dovecot 2.x and frequently see delays "loading message". I admit to having a very large mailbox and running multiple imap clients. I've turned on thunderbird debugging, and at least one pattern I see in the log during the delays (maybe a minute long) is the failed protocol message above. These are no the only delays, but are the ones

Hello, I finally run my measurement model in sem - successfully. Now, I am trying to print out the path diagram that is based on the results - but for some reason it's not working. Below is my script - but the problem is probably in my very last line: # ANALYSIS OF ANXIETY, DEPRESSION, AND FEAR - LISREL P.31 library(sem) # Creating the ANXIETY, DEPRESSION, AND FEAR intercorrelation matrix

Hello everyone. I would like to simulate a small map area and store information regarding it. I am trying to find out what might be the best way to do that and I need some feedback, as you might be more experienced. Usually a map can be described by a k * m dimensions matrix where every cell can store one value(eg. number) that expresses something for that area. In my map I would like to

hi, i'm trying to perform a clustering on a big dataframe the code is this: print("load required R packages") require(spgrass6) require(cluster) gmeta6 <- gmeta6() print("read in our 7 raster files from GRASS") x <- readFLOAT6sp(c ("er","crosc","longc","slope","profc","minic","maxic"))

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24),