similar to: Fitting Maximums of data series with cubic spline

Hi Users, I have daily series of data from 1962 - 2000, with the data for February 29th in leap years excluded, leaving 365 daily values for each year. I wish to convert the daily series to monthly series. How can I do this using the zoo package or any other package? Thanks ZABLONE OWITI GRADUATE STUDENT Nanjing University of Information, Science and Technology College of International

Dear Users, I have a monthly data for a number of years(1960-2007) for a number of stations and i wish to extract sesonal time-series for the months of March-May and October-November for very station. I have read this data with read.table in R with stations as columns and time (months) as rows. My attempt to aggregate with the zoo package using the function as.yearqtr failed since this sums

Dear users, I have several columns of data (each column containing monthly data for a particular year from january - december) . I would wish to combine the columns to get a Single column of continuous data as shown in (b) below. I have read this data as table in R (a) Data example Year 1903 1904 1905 1906 Jan 125.0 30.0 113.0 5.0 Feb 128.0 100.0 70.0 388.0

Dear users, I have monthly station data (44) stations data for 45 years which I have read in R using read.table. the data is in the format: Year Month Station1 Station2 ...................................... Station 44 (i.e the column names in the 1st row), I also have the latitude and longitude of the stations in a separate file in R (in the format : Station LAT LON). I wish to

Hi R users, What is the command for adding (or how do we add) country boundaries on r-spatial plots. Thanks ZABLONE OWITI GRADUATE STUDENT Nanjing University of Information, Science and Technology College of International Education Add: 219 Ning Liu Rd, Nanjing, Jiangsu, 21004, P.R. China Website: www.nuist.edu.cn ==================================================== [[alternative

Dear Users, I wish to know at what confidence level is the confidence interval provided in the Spectrum function (plot.spec) plots. The only information provided in the help regarding this is : a confidence interval will be plotted by plot.spec: this is asymmetric, and the width of the centre mark indicates the equivalent bandwidth. Attached is an example plot with the BLUE confidence

Dear all, I am plotting daily data from 1950 to 2000 and I wish to have the time in X-axis to indicate years and NOT days (i.e 1950 1951,......,2000 instead of 1,2, 3, ......). How do I set the x-axis to show years instead of days. ZABLONE OWITI GRADUATE STUDENT Nanjing University of Information, Science and Technology College of International Education Add: 219 Ning Liu Rd, Nanjing,

Hi, all I am trying to get the basis matrix and penalty matrix for natural cubic splines. In the "splines" package of R,"ns" can generate the B-spline basis matrix for a natural cubic spline. How can I get the basis matrix and penalty matrix for natural cubic spline. Thanks a lot! Lee [[alternative HTML version deleted]]

Hello everone,            Anyone who knows how to force a cubic smoothing spline to pass through a particular point?            I found on website  someone said that we can use "cobs package" to force the spline pass through certain points or impose shape           constraints (increasing, decreasing). However,  this package is using  B-spline and can only do linear and quadratic

Hallo, I'm facing a problem and I would really appreciate your support. I have to translate some Matalb code in R that I don't know very well but I would like to. I have to interpolate 5 point with a cubic spline function and then I expect my function returns the Y value as output a specific X value inside the evaluation range. Let's suppose that: 1- *X = [-10, -5, 0, 5, 10]* 2

Hello,   does anyone have an example on how to use restricted cubic splines function rcs within survfit.cph, if cph (Cox Proportional Hazard Regression) was done with restricted cubic splines (which I made to work)? Thank you. > [[alternative HTML version deleted]]

(corrected version of previous posting) I fit a GAM to turtle growth data following methods in Limpus & Chaloupka 1997 (http://www.int-res.com/articles/meps/149/m149p023.pdf). I want to obtain figures similar to Fig 3 c & f in Limpus & Chaloupka (1997), which according to the figure legend are "expected size-at-age functions" obtained by numerically integrating

I have used LME to fit a mixed effects model on my data. The data has 274 subjects with 1 to 6 observations per subject. Time is not linearly associated with the outcome, so I used ns to fit a natural cubic spline with 3 auto knots. Subject and the natural cubic time of spline are both treated as random effects. This model has run without any problem, but now I would like to plot trajectories for

I'm trying to use R to determine the quality of a cubic Bezier curve approximation of an elliptical arc. I know the four control points and I want to compute (x,y) coordinates of many points on the curve. I can't find anything in either the base distribution or CRAN that does this; all the spline-related packages seem to be about *fitting* piecewise Bezier curves to a data set.

I want to compare the fit of a quadratic model to continuous data, with that of a cubic spline fit. Is there a way of computing AIC from for e.g. a GAM with a smoothing spine, and comparing this to AIC from a quadratic model? Cheers ****************************************** Tom Reed PhD Student Institute of Evolutionary Biology 102 Ashworth Laboratories Kings Buildings University of

Greetings, I would like to use a cubic spline or smoother to model the fixed effects within nlme. So far the only smoother I have been able to get to run successfully in nlme is smooth(). I tried smooth.spline: fixed=list(lKa~1,lCL~smooth.spline(BSA, df=3)) the error I got was the following. Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid

I have a question about the basis functions of cubic regression spline in mgcv. Are there some ways I can get the exact forms of the basis functions and the penalty matrix that are used in mgcv? Thanks in advance! Yan [[alternative HTML version deleted]]

Hello,        Dr. Simon Wood told me how to force a cubic spline passing through a point. The code is as following. Anyone  who knows how I can change the code to force the first derivative to be certain value. For example, the first derivative of the constrained cubic spline equals 2 at point (0, 0.6).        I really appreciate your help!        Thanks!                 Best             Victor   

Dear Thomas, I want to use evaluate effect of Age using restricted cubic form in the FGR function as Fgr.crr <- FGR(Hist(time, event) ~ rcs(Age_years), data=dat) It provides error. " Error in parse(text = termtext, keep.source = FALSE): .... 1: response ~ rcs(Age_years Do I need to change any of the R code? Regards Amalraj Raja The University of Aberdeen is a charity

Hi, I'm using the cubic splines from termstrc package. I invoked the splines_estim function with a group of 43 bonds. It computes 6 knot points and returns values for alpha1 to alpha7. My question is how to use these alpha1 to alpha7 in the equation of yield? For example, if I'm trying to find the yield at, say, 12.25 years, which falls between 3rd & 4th knot points what should the