[R] contour plot

Hi, Could someone give some ideas on plotting a contour by using geoR package, please? Thank you Kagba [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r

[R] Differentiation on mathematical equation

Dear All, I was trying to compute the first and second order differentiation on a number of mathematical equations, but never found any command in R to do this operation. At present, i have almost 3000 different functions, consist of polynomial, harmonic and several other functions. Could some

[R] Color index in image function

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. ima

[R] Color index in image function

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. ima

Re: [R] Color index in image function

alues for intance, what are these colors represent for? Does the first color, "#4C00" represent the lower/medium/higher temperature? I do hope you could advice me on this matter? Thank you Fir - Original Message From: FMH To: Steve Lianoglou Sent: Sunday, September 6, 2

Re: [R] Color index in image function

color index in another table beside the image, started from green followed by green-yellow, yellow, yellow-blue and blue? Could someone please advice on this matter? Cheers Fir - Original Message From: Bernardo Rangel Tura To: FMH Sent: Monday, September 7, 2009 11:13:12 AM

[R] Derivative of nonparametric curve

Dear All, I'm looking for a way on computing the derivative of first and second order of a smoothing curve produced by a nonprametric regression. For instance, if we run the R script below, a smooth nonparametric regression curve is produced. provide.data(trawl) Zone92   <- (Year == 0 & Zone =

[R] Data in Array

Dear All, I have some data which were stored in few  matrices with different orders. Let have three different matrices a, b and c, which have the same number of column but different number of row. a <- matrix(1, nrow = 5, ncol = 1) b <- matrix(2, nrow = 10, ncol = 1) c <- matrix(3, nrow = 15,

Re: [R] Data in Array

Yes, but what actually i want to have is an array that might store these different matrices. - Original Message From: Schalk Heunis To: FMH Cc: r-help@r-project.org Sent: Tuesday, September 8, 2009 11:16:29 AM Subject: Re: [R] Data in Array have you tried rbind? On Tue, Sep 8

[R] gridlines in contour plot

Dear All, Someone suggesting me to use the filled.contour function to plot the image together with the color index, and an example from the help menu is show below.   # require(grDevices) # for colours filled.contour

Re: [R] Derivative of nonparametric curve

Thank you - Original Message From: spencerg To: "Liaw, Andy" Cc: Rolf Turner ; FMH ; r-help@r-project.org Sent: Wednesday, September 9, 2009 3:08:43 PM Subject: Re: [R] Derivative of nonparametric curve     This may be overkill for your application, but you might be int

Re: [R] Color index in image function

Thank you From: Henrique Dallazuanna Cc: r-help@r-project.org Sent: Sunday, September 6, 2009 4:30:14 PM Subject: Re: [R] Color index in image function Try this to see which colors the codes are: Color <- function(color){     z <- matrix(1:length(color), nco

Re: [R] Data in Array

Thank you - Original Message From: jim holtman To: FMH Cc: Schalk Heunis ; r-help@r-project.org Sent: Tuesday, September 8, 2009 3:45:09 PM Subject: Re: [R] Data in Array Not sure what you mean by 'store', but you can use a list: > a <- matrix(1, nrow = 5, ncol = 1)

[R] How to combine matrices?

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 __

[R] How to sort the elements in a matrix?

Dear All, Let ab is the combination of matrices a and b. a <- matrix(c(1,2,3,4), nrow = 2, ncol = 2) b <- matrix(c(4,3,1,2), nrow = 2, ncol = 2) ab <- rbind(a, b)   >From matrix ab, could someone give some advice on the way to sort all elements >in column 2, by the sequence of the elements in c

[R] Problem in graph plotting

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') #

Re: [R] Problem in graph plotting

write the coding? Thank you Fir   - Original Message From: jim holtman To: FMH Cc: r-help@r-project.org Sent: Wednesday, September 23, 2009 1:04:21 PM Subject: Re: [R] Problem in graph plotting try this: plot(tp,dp, type= 'l',ylim=rev(range(dp))) On Wed, Sep 23, 2009 at 7

[R] How to extract one of four plots in a linear regression model

Dear All, A linear regression model could be fitted by using lm function and the plot function can be used to check the assumption of the model. The help menu shows few instances on suitable coding for fitting such a linear model. In addition, four different plots could be extracted simultaneo

[R] How to add a title to represent four different plot in lm function

Dear All, A linear regression model could be fitted by using lm function and the plot function can be used to check the assumptions of the model. The example is as followed. require(graphics) ## Annette Dobson (1990) "An Introduction to Generalized Linear Models". ## Page 9: Plant Weight Data.

[R] Changepoints estimation in a data series

0, 60 >From the above series, i reckon there is more than one changepoint and  >presuming there is a package in R which might enable the estimation on such >changepoints. Could someone please advice me on this matter by using R? Cheers, FMH _

[R] Problem with GroupedData

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

[R] Inflection point on a curve

Dear All, I have been trying to find an inflection point from a nonparametric model, for instance on  a series of Economics, Financial, Environmental data, which was fitted via sm.regression package, but have difficulty to find the derivative of the model. The following are the data and the cod

[R] Derivative of a smooth function

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 th

[R] Third and fourth order of derivative in smooth.lf function from locfit package

Dear All, Could someone please advice me the way to define the derivative of the local polynomial regression in third and fourth order from the smooth.lf function? Thank you Fir [[alternative HTML version deleted]] __ R-help@r-project

[R] Color of the plot which correspond to the group of the observations

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 a

[R] How to speed up R with version 2.9.2?

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

Re: [R] How to speed up R with version 2.9.2?

Thank you for your answer. I'm using Win XP with 2GB RAM in memory. Cheers Fir - Original Message From: stephen sefick To: FMH Cc: r-help@r-project.org Sent: Fri, October 2, 2009 4:38:10 PM Subject: Re: [R] How to speed up R with version 2.9.2? You're fine, but please d

Re: [R] How to speed up R with version 2.9.2?

Yes, i noticed there are few Windows start up programs, shown by the EasyCleaner software. Cheers   - Original Message From: Gabor Grothendieck To: FMH Cc: stephen sefick ; Sent: Fri, October 2, 2009 4:50:52 PM Subject: Re: [R] How to speed up R with version 2.9.2? Its under 5

[R] Slope between two points

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

[R] How to use multiple time series in R?

Dear All, I have temperature series for 10 different sea levels, where there are 3000 observations for each level. I am planning to use multiple time series modelling on these series, with/ without correlation structure for the residuals. Is there any package that can be used to do this operat

[R] Estimation in a changepoint regression with R

Dear All, I'm trying to do the estimation in a changepoint regression problem via R, but never found any suitable function which might help me to do this. Could someone give me a hand on this matter? Thank you. __ R-help@r-project.org mailing li

[R] How to find the interception point of two linear fitted model in R?

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: ##

[R] Decomposition method

Dear All, I'm doing an  imputation on the missing data and is looking for a decomposition method in R. Could somone advice me the way to do this? Thank you Fir __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Additive decomposition method

Dear All, I'm doing an  imputation on few missing data and is looking for an additive decomposition method in R. Could someone please advice me the way to do this? Thank you Fir __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

Re: [R] [R-SIG-Finance] Additive decomposition method

Hi, Thank you for the function. I was trying to use the decompose function, but the program cannot run as there were several missing values in my data frame. Could someone advice me the way to handle on this matter? Thanks F __ R-help@r-proje

[R] tapply function

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 fucntio

[R] A comparison between span in lowess method and h in sm.regression

Dear All, I have fitted a bivariate and an additive model on a monthly temperature from 1987-2007 via nonparametric regression and would like to compare between these models via F-test, but having difficulty in determining the possible values of span with lowess method. For instance, these two

[R] Delta Kronecker

Hi, Could some give some ideas on how to compute the spatiotemporal covariance matrix by a sum of Kronecker products in R. Is there any special function that can be used? Cheers. Firdaus [[alternative HTML version deleted]] __ R-he

[R] Spatiotemporal correlation function

Hi, I'm trying to compute  the spatiotemporal correlation matrix by using Delta Kronecker products of spatial and temporal correlation matrix  in R, but didn't find any delta Kronecker's operator in R. The operators in matrix such as multiplication, addition, eigen values/vector and etc is easi

[R] Correlation structure in AR(2) and ARMA(p,q)

Hi, Let theta be the coefficient in AR(1) model, the components of the correlation structure is equal to 1 for times s = t, and theta^|t-s| for t > s. Could someone do me a favour to give some ideas on the correlation structure in AR(2) model and correlation structure of ARMA(p,q), as well? T

[R] Function in R for computing correlation matrix and covariance matrix

Hi, At present, i have two distinct and real values for the coefficient, which is   required in AR(2) model. Based on my revision, for distinct and real values of the coefficients in AR(2) model, the correlation structure separated by lag h can be computed by p(h) = a*z1^(-h) + b*z2^(h), where p

[R] p-value for the parameter in ARIMA model with R

Dear All, I decided to use an AR(1) model for the residual series and trying to find the p-value for each parameter by using arima command in R, but i never find it from the output. The output gives me the parameter and mean's value, its standard error, estaimated variance, AIC and loglikehood,

[R] standard error and p-value for the estimated parameter in AR model

Dear All, I used an  AR(1) model to explain the process of the stationary residual and have used an 'ar' command in R. From the results, i tried to extract the standard error and p-value for the estimated parameter, but unfortunately, i never find any way to extract  it from the output. What

[R] Inverse matrix

Dear All, Let A is a matrix which is: A <- matrix(c(1,2,3,4),nrow=2) How could we find the inverse of A? I try to use ginv(A), but it didn't worked. Thanks Fir __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Interception point between two lines

Dear All, Let mod2 and mod3 are two regression equations representing two distinct lines and i'm keen to find the intreception point between these two lines and the following are part of the codes.     m1 <- as.matrix(rbind(coef(mod2), coef(mod3))) a <- cbind(c(1,1), -m1

[R] How to seperate date and time into different columns?

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.

[R] Color intervals in image function (image.plot)

Hi, The script below is my current coding in order to produce a contour plot of temperature across altitude and time. In my case, time,altitude and temperature are represented by x, y and z variables. ## Brazilan.Pallete <- colorRampPalette(c("blue",

[R] Color intervals in image.plot function

Hi, The script below is my current coding in order to produce a contour plot of temperature across altitude and time. In my case, time,altitude and temperature are represented by x, y and z variables. ## Brazilan.Pallete <- colorRampPalette(c("blue"

[R] SARIMA model

Dear All, Could someone please advice me the appropriate package for fitting the SARIMA model? Thanks Fir __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/

[R] NA with lmList

Dear All,   I was trying to use the lmList function to get the lmList graphic but  have a problem creating the graphic.  My data has missing values and an error occurred  as stated below. ### library(nlme) > a   schoolid spring score childi

[R] Kalman filter

Dear All, Could anyone give me a hand to suggest few packages in R to running Kalman prediction and filtration ? Thanks Fir __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://ww

[R] How to change font size in plot() function

Dear All, Could someone please advice me the way to change the size of the title and x and y-label in plot() function. I've tried to use 'font' and 'font.main' call in plot() function, but it didn't make any changes in terms of the size. Thanks Fir

[R] best subset regression in R

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]] __ R-help@r-project.org mailing list

[R] Box-Cox transformation in R

Hi,   Could any one please help how I can transform data based on Box-Cox Transformations in R.   Any helps will be much appreciated.   thanks, Kagba [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.c

[R] cluster analysis on extreme event

Dear all, I'm modelling extreme rainfall,particularly those that lie above a threshold & was searching for a suitable package in R which may enable a cluster analysis on those extreme events and would really appreciate for any suggestions. Thanks, Fir __

[R] smoothing parameter in locfit package

Hi, In the locfit package, could someone please let me know the automatic selection of smoothing parameter if Gauss kernel density function is used as weight function? thanks Fir [[alternative HTML version deleted]] __ R-help@r-proj

Re: [R] smoothing parameter in locfit package

s the value of the bandwidth is fixed (by default) in the package, unless we specify our own value? Cheers, Fir - Original Message From: "Liaw, Andy" To: FMH ; r-help@r-project.org Sent: Mon, May 10, 2010 4:00:51 PM Subject: RE: [R] smoothing parameter in locfit package

[R] Optimal choice of the threshold u in Peak Over Threshold (POT) Approach

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 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

[R] How to find points of intersection

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 __ R-help@r-project.org mailing list https://stat.ethz.ch/m

Re: [R] How to find points of intersection between harmonic function and a line

. Thank you, Fir - Original Message From: Uwe Ligges To: FMH Cc: r-help@r-project.org Sent: Tue, February 22, 2011 1:11:41 PM Subject: Re: [R] How to find points of intersection On 22.02.2011 12:27, FMH wrote: > Dear All, > > I'm looking an appropriate way in R to co

[R] How to find points of intersection between harmonic function and a line

Hi, Sorry for the very short explanation about the problem of intersection. I have a wave pattern monitored from the heart beat in a particular interval of times. Apart fom that, there is a line with positive slope (e.g: y = x+2) which lies across the wave and intersect on a number of points.

[R] How to find points of intersection between harmonic function and a line

Hi, Sorry for the very short explanation about the problem of intersection. I have a wave function monitored from the heart beat in a particular interval of times. Apart fom that, there is a line with positive slope (e.g: y = x+2) which lies across the wave and intersect on a number of point

[R] Extremum index from a sampling time series data

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  fin

[R] goodness of fit of a mixed model and function in R

Dear All,   Could anybody please advice me the way to check the goodness of fit of the linear mixed-effects model and the suitable function in R to do so, and thank you in advance for your assistance.   Fir [[alternative HTML version deleted]]

[R] R2 for a mixed-effects model with AR(1) error structure

Dear All, The following equation is a linear mixed-effects  model with linear trend and AR(1) error structure, y = B0 + B1x + bo + b1x + e; e~AR(1)  where  y is a response, x is the predictor, B0 and B1 are fixed effects and b0 and b1 are random effects. Coud someone please advice me a fun

[R] Dispersion Index in POT

Dear All, I would like to plot the dispersion index against lambda and dispersion index against threshold. Appreciate if someone could help me to extract dispersion index in POT package. Pls also let me know the way to calculate it manually for a skewed distribution. Thank you [[alterna

[R] Continuous Wavelet tranform in Biwavelet Package

Dear All, I have 15 minutes of temperature data for 3 years and would like calculate and plot the CWT using wt command in Biwavelet package. Appreciate if someone could tell me the way to define it in wt command. Thank you, [[alternative HTML version deleted]] __