similar to: cluster analysis on extreme event

I am a student doing my MSc Research Methods, i am working on my thesis research on analysing and modelling of crop failure risks due to drought in selected districts in Malawi. The analysis and modelling will focus on two crop stages of development: just after planting and flowering stages. I have rainfall, temperature, Relative humidity, wind speed, coordinate location data for the past 30 years

Dear All, Let: dp: depth of the?river tp: temperature with respect to depth We can?have a simple scatter plot, between depth as y-axis and temperature as x-axis,?by using a?plot function as shown?below. #####################? dp <- c(1,4,3,2,5,7,9,8,9,2) tp <- 1:10 plot(tp,dp, type= 'l') ##################### ? Could someone advice me on?the way to?plot the same pair of

Dear All, Let a, b and c are?three matrices with same no. of column but different no. of row. a <- matrix(1, 1, 2) b <- matrix(2, 2, 2) c <- matrix(3, 3, 2) Could someone help me to combine these matrices together as a single matrix? Thank you Fir

Dear All, I was looking for the color index in?image function,?such as from topo.colors(n) and etc. but still never found it. For instance, from the help menu. ########################################### # Volcano data visualized as matrix. Need to transpose and flip # matrix horizontally. image(t(volcano)[ncol(volcano):1,]) # A prettier display of the volcano x <- 10*(1:nrow(volcano)) y

Dear All, I have a series of data in which the first column consist of a combination of date and time, for instance 17 April 2008 at 4.01pm, such data is recorded as: 4/17/2008 16:01 I'd like to seperate it into four different columns which consist of Day, Month,Year and Time, respectively. Could someone please advice me on this mater? Thank you, Fir

Dear All, Could?someone please suggest me the way to calculate the optimal threshold in POT method via?any available ?packages in R? Thanks, Fir

Dear All,   I am working on a time series of hourly river flow measurements from 2000 - 2003  and have been trying to compute the extremum index from the original series as well as a new series from  sampling with replacement . The extremum Index produced from the original data series looks  fine but sadly, a strange value of index is given from the sampling. The flow measurements are in range of

Dear All, I've been?searching for?appropriate codes to compute the rate of change and the curvature?of ?nonparametric regression model whish was denoted by a smooth function?but?unfortunately?don't manage to?do?it. I presume that such characteristics from a smooth curve can be determined by the first and second derivative operators. The following are the example of fitting a

Dear All, Could anyone?give me a hand?to suggest few packages in R to running Kalman prediction and filtration ? Thanks Fir

Dear R helpers, I need to use KS and AD test for Generalized Pareto and Generalized extreme value. E.g. if I need to use KS for Weibull, I have teh syntax ks.test(x.wei,"pweibull", shape=2,scale=1) Similarly, for AD I use ad.test(x, distr.fun, ...) My problem is fir given data, I have estimated the parameters of GPD and GEV using lmom. But I am not able to find out the distribution

Dear All, Let? 499 piece-wise lines were buit up by 500 pair of observations, via R code below. x <- 1:500 y <- rnorm(500) plot(x, y, type = 'b') I was trying to compute all the slopes?for the ?lines which were connected between two adjacent points. For instance,?slopes?of lines between first and second points, second and third points, and so on, until?between??points 499 and 500,

Dear All, I'm sorry if my question does not suit with this R group. I have recently installed?R software with version 2.9.2, but i found the program took almost 1 minute as soon as it was opened, before it can be used. However,?the previous version 2.9.1 only take few seconds?after the menu bar was clicked. This circumstance has caused me to wait for couple of minutes as several R windows

Dear All, I have a question about package "ismev", its function "gpd.fit", and interpretation of the results. I used the package ismev to do an extreme value analysis on a fire dataset. Two variables are used in the analysis. The focal variable is acreage burned per fire, ranging from 1 to 5000 acres per fire. In total, there are 69,980 observations. The date covers

Dear All, ? Let: dp: depth of the river tp: temperature with respect to depth These pair of observations are in 3 different groups i.e: Obs. 1,3,5,7?from the first group Obs. 2,4 and?10?from second group Obs 6,8 and 9?from third group. We can have a simple scatter plot, between depth as y-axis and temperature as x-axis, with each pairs are denoted by a red dot, by using a plot function shown

Hi,   I have an original data frame with 8 columns of variables, which are stored in 'data1' frame.   data1 <- read.csv("E:\\PHD GLASGOW UNIVERSITY\\Data\\R\\Colin\\Cailness21.csv") attach(data1) names(data1)   [1] "Date"       "d"          "m"          "y"          "Time"      [6] "Depth"      "Temp"      

Dear All, Let have 10 pair of observations, as shown below. ###################### x <- 1:10 y <- c(1,3,2,4,5,10,13,15,19,22) plot(x,y) ###################### Two fitted? models, with ranges of [1,5] and [5,10],?can be easily fitted separately by lm function as shown below: ####################### mod1 <- lm(y[1:5] ~ x[1:5]) mod2 <- lm(y[5:10] ~ x[5:10]) #######################

Hi, I tried to use tapply function to find the mean of the data?in each group as the following command, but the result are NA, as there are several missing values in each group. tapply(data,group,mean) Could someone please advice me the way to ?ignore the missing data in order for the fucntion to run successfully? Thanks Fir

I have a time series of rainfall in a dataframe. I would like to be able to aggregate this using a sliding window approach- i.e. a new 24 hourly total is calculated for each hours rainfall. I'm struggling to understand how this might be achieved - currently I've tried looping a sum function to re-calculate a new total at every stage of the loop. for (inp[[9]] in

Dear All,   Could someone please give some advice the way to do linear modelling via best subset regression in R? I'd really appreciate for your kindness.   Thanks, Kagba [[alternative HTML version deleted]]

Dear All, I'm looking an appropriate way in R to compute/estimate? points of intersection between a line and a curve and will really appreciate for any suggestion or ideas? Thank you, Fir