similar to: Visualising Vectors

I am wanting to plot a series of wind vectors onto a contoured area map for a series of weather stations (eg arrows showing wind speed/direction for a particular time snapshot), can someone please advise me how best to approach this? My desired end point is to be able to link a time series of such data together so that I will in effect have a "movie" displaying the evolution of these

Quick question - I am having problems creating an orographic image (similar to volcano example). I have created a map matrix with 3 columns and over 2million rows. I have created the matrix as follows: map<-read.table("map.dat",header=TRUE) long.grid<-sort(unique(map$long) lat.grid<-sort(unique(map$lat) map.matrix<-matrix(map$height,nrow=length(lat.grid),byrow=TRUE) The

Probably a stupid question, but I don't seem to be able to find the answer I'm looking for from any of the R literature. Basically I have a matrix with several thousand rows and 20 columns(weather stations) of wind direction data. I am wanting to extract a matrix which contains data for all columns conditional on column 20 having a value of _either_ less than 45 or greater than 315. (ie I

Dear all: Does anyone have any suggestions on how to make a spatially explicit landscape with spatial autocorrelation in R? In other words, a landscape where all cells have a spatial reference, and the environment values that are closer in space are more similar (positive spatial autocorrelation). Thank you, Laura

Hello, I have daily wind speed data and need to fit seasonal ARIMA model, problem is that my period is 365. But when I use arima(...) function, with period 365, I?m getting error message: ?Error in makeARIMA(trarma[[1]], trarma[[2]], Delta, kappa) : maximum supported lag is 350?. Can someone help me with this problem? Thank you Sincerely yours, Laura Saltyte

Hi, this may sound like a very stupid question and perhaps it is, so I apologise in advance if anyone feels bored by it. I have done some principal component analysis with prcomp and what I would like to do is to visualise the effect of a principal component, i.e. show a graph of my data and then show how that graph would change if I go, let's say, one standard deviation, in either

Hi: Summary: I am trying to determine the 90th percentile of ambulance response times for groups of data. Background: A fire chief would like to look at emergency response times at the 90th percentile for 1 kilometer grids in Cape Coral, Florida. I have mapped out ambulance response times on a GIS map. Then I superimpose a regularly-spaced grid over the response times and spatially join the

Dear R users, I have used the following function (in blue) aiming to find the linear regression between MOE and XLA and nesting my data by Species. I have obtained the following results (in green). model4<-lme(MOE~XLA, random = ~ XLA|Species, method="ML")summary(model4) Linear mixed-effects model fit by maximum likelihood Data: NULL         AIC     BIC   logLik  -1.040187 8.78533

Hello, Is there a simple way of stacking/merging two dataframes in R? I want to stack them piece-wise, not simply add one whole dataframe to the bottom of the other. I want to create as follows: x.frame: aX1 bX1 cX1 ... zX1 aX2 bX2 cX2 ... zX2 ... ... ... ... ... aX99 bX99 cX99 ... zX99 y.frame: aY1 bY1 cY1 ... zY1 aY2 bY2 cY2 ... zY2 ... ... ... ... ... aY99 bY99 cY99 ...

Hi all, I've been trying to figure out how to superimpose a histogram on a biplot that shows the relative contribution of each axis. I have been using the NIPALS function (http://biomserv.univ-lyon1.fr/~dray/ files/softwares/nipals.R) to run principal component analyses. Here is a toy example. source("http://biomserv.univ-lyon1.fr/~dray/files/softwares/nipals.R")

Hello, I'm trying to superimpose a line plot onto a histogram but I'm not having any luck. I've attached the dataset. What I did was: > hist(data,freq=F) Now I'm trying to superimpose the following points with a line connecting them onto the histogram: x y 100 0.535665393824959 200 0.212744329736556 300 0.0844933242968584 400 0.0335572838043417 500

Dear R Help Community, I have a question and an answer (based on reading this forum and online research), but I though I should share both since probably there's a much better way to go about my solution. My question is specifically about how to best visualise multiple response contingency tables. What I mean by 'multiple response' is that the total number of responses per row of a

Hi, Just a quick query - if I'm creating a function to produce a number of histograms per page of output (one per column from a matrix), how can I pass the column name of the matrix into the title (or indeed to form part of the x-axis label)? TIA, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax:

Ok, I have reinstalled R-1.9.0 and this appears to have fixed the problems I was having with png(). However, I have a further question regarding png() Is it possible to pass a par() argument to the png() command? I am wanting to produce 4 plots per object, which I normally acheive on an X window by par(mfrow=c(1,4)). I have tried calling a new plot and setting par in this way but this has no

Hello, I requested help a couple of weeks ago creating a dipole field in R but receieved no responses. Eventually I opted to create a 3d sinusoidal plot and concatenate this with its inverse as a means for a "next best" situation. It seems that this isn't sufficient for my needs and I'm really after creating a continuous 3d gaussian mesh with a "positive" and

Dear R users, I am in trouble using the function solve(matrix) in a loop. When it happens that I have a singular matrix, I get an error message saying so and the loop stops. My question is if there in any way to keep on going in the loop. Thank you in advance, Laura *---------------------------------------------------------* Laura Antolini, PhD Department of

Dear R People: I used the "Surv" function to produce the following object: >a <- Surv(1:4,2:5,c(0,1,1,0)) a [1] (1,2+] (2,3 ] (3,4 ] (4,5+] >str(a) Error in "[.Surv"(object, 1:ile) : subscript out of bounds > Why does str(a) give an error, please? Or did I do something wrong? Thanks in advance. R Version 1.9.1 Windows Sincerely, Laura Holt mailto: lauraholt_983

??????? Original Message ??????? On Thursday, April 11, 2019 12:55 AM, John Fawcett via dovecot <dovecot at dovecot.org> wrote: > On 11/04/2019 00:51, Laura Smith via dovecot wrote: > > > ??????? Original Message ??????? > > On Wednesday, April 10, 2019 11:48 PM, John Fawcett via dovecot dovecot at dovecot.org wrote: > > > > > On 11/04/2019 00:18, Laura

Dear list, I just tried to use the function cluster.stat in the package fpc. I just have a couple of questions about the syntax: cluster.stats(d,clustering,alt.clustering=NULL, silhouette=TRUE,G2=FALSE,G3=FALSE) 1) the distance object (d) is an object obtained by the function dist() on my own original matrix? 2) clustering is the clusters vector as result of one of the many clustering methods?

Hello there, I'm trying to plot vectors with p<0.1 in a NMDS ordination plot using p.max. Below the scripts I'm using. I guess I'm missing something! could you please give me a hand? species<-metaMDS(species_matrix)ef<-envfit(species,environmentaldata_file,permu=999,na.rm=TRUE)efplot(species, dis="sites")plot(ef,p.max=0.1) Error in plot.envfit(ef, p.max = 0.1) :