similar to: Changepoint analysis--is it possible to attribute changpoints to explanatory variables?

Hello Does anyone know how to get the p values for changepoints in the "changepoint" package? There is guidance in the CRAN vignette https://cran.r-project.org/web/packages/changepoint/changepoint.pdf but I can't make it work to get the p value out... Thanks Nick Wray [[alternative HTML version deleted]]

Hi, all, I'd like to find one or more changepoints in a dataset eg for data like 0.9 1 0.8 0.8 0.85 0.3 0.9 0.9 0.9 1.1 0.3 is unlikely a changepoint but for data like 0.9 1 0.8 0.8 0.85 0.3 0.2 0.2 0.1 0.2 0.3 0.3 is a change point I am wondering is there any function can do changepoint analysis? Ying Fan

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.

Hi All, I'm using the following code to produce some stacked bar graphs. *setwd("C:\\Users\\Tinus\\Documents\\NMMU\\R\\Seamounts")* *SChla <- read.csv("SM_Chla_data.csv")* * * *#Extract mean values from data file* * * *Coral <- SChla[185:223,] #Reduce SChla to Coral only* *coral <- with(Coral , aggregate(cbind(Pico, Nano, Micro), list(Depth),FUN=mean))*

Dear R Community, I am trying to optimize a water quality model that I am using. Based on conversations with others more familiar with what I am doing I plan to implement DEOptim to do this. The water quality model is interfaced through a GUI. I have the input file necessary to alter parameters and run the model as a text file. To do the optimization I have figured out the general procedure

Sorry for disturbing you with another newbie question! I have a data frame about coastal waters quality parameters: for some parameters (e.g. NH3) I have only 1 observation for each sampling station and each sampling date, while in other cases (chlorophyll) I have 1 obs for each meter-depth for each station and date. How can I select only the max chlorophyll value for each station/date? example

Hi, I would be very grateful for some advice. I have read the help pages for Date, strptime, etc. All examples seem to use some version of day month year as date format. However I have Weekly composite data so ONLY want to input the dates as Week.Year (eg 35.2011). strptime seems to show this is possible using %W for week (UK convention) and %Y for year. My data is in a df called chlorophyll and

Hello. I have made many pngs in the past, but started having trouble with them yesterday. I just downloaded R 2.13.0 this morning, and was able to make some pngs earlier this morning. Now I am trying to remake them, but get the same error message I got yesterday: > u='in' > w=10.5 > h=8 > png('winbugs changepoint graphs ind hibs gt 400

Hi, I am trying to run a 1D nutrient-phytoplankton-zooplankton model in R using the package 'deSolve'. The code is shown below: DEPTH = seq(2.5, 147.5, 5) NPZ = function(t, state, params){ with(as.list(params), { P <- state[1:NB] Z <- state[(NB + 1): (2*NB)] N <- state[(2*NB + 1): (3*NB)] F.I = function(z, hr){ I0 = function(hr){

I am trying to understand my population abundance data and am looking into analyses of change point to try and determine, at approximately what point do populations begin to change (either decline or increasing). Can anyone offer suggestions on ways to go about this? I have looked into bcp and strucchange packages but am not completely convinced that these are appropriate for my data. Here is

I am having difficulty understanding the syntax of the nlme() function for nonlinear mixed models. The data frame is called Marouane.chlorophyll. The model involves a dependent variable (Absorb) and an independent variable (Ch.surf), which are both numeric variables in the data frame. The data are hierarchically grouped as Espece, Plante nested within Espece, and Feuille nested within Plante

Hello The other say I posted a question about altering confidence limits in the seqMK() function in the pheno package. I didn't get any bites but in any case I've been trying to use this function. It seems to me to give very odd results - sometimes identifying 8 or 9 changepoints in a relatively short vector of maybe 40 elements, and at other times not giving any changepoints at all,

Dear users I am fitting a Generalized Additive Mixed Models (gamm) model to establish possible relationship between explanatory variables (water temperature, dissolved oxygen and chlorophyll) and zooplankton data collected in the inner and outer estuarine waters. I am using monthly time-series which are auto-correlated. In the case of the inner waters, I have applied satisfactoryly (by

Hello, I am a new user to R and have a(proabably very basic) question relating to a boxplot with eleven sites in a river system where i measured Chlorophyll a concentrations. The eleven sites are made up of three sites above a point source pollution site; three sites below this pollution source and five sites in a different unaffected river. The following code is for a simple plot of the data

Dear gracious R community, I would like to produce charts of phytoplankton biomass changes through time. Each species has a line, and the biomass varies in mirror form along the line for each species along the X time axis. Here is an example of what I'd like to do: http://www.er.uqam.ca/nobel/r30240/Succession.jpg Any suggestions? Thanks David Bird UQAM, Montreal

Hi everyone. I have a long label that I would like to split. I found that I could use "strwrap" for simple text. However, this is not working with this label: str = expression(paste("< 20 ?m phytoplankton ","(cells ? ",mL^-1,")")) plot(...., ylab = strwrap(str,20),...) I suspect this is because I'm using "expression" for form my label.

Hello does anyone know how to set the confidence level within the seqMK() function in the pheno package. It seems to be set automatically at 0.05 and there seems to be neither an input function to set a different level, as there are with some changepoint functions, nor an output to give one the p-value of the changepoint identification which you get with functions ie func$p_value etc ... Thanks

Hello all, I can't seem to figure how to use a greek character in expression() in plot() labels without adding a space. So for example below when plotting this out x<-1:10 plot(x,x^2, xlab=expression(Chlorophyll~italic(a)~mu~g~cm^-2)) the axis label read as μ g cm^-2 because I have space there with a tilda. But if I remove the tilda then my units are mug cm^-2. Can anyone recommend a

DeaR R-fellows, I got a lot of help from you concerning the generation of heat maps, so I was able to write and work with a script for this purpose at last. Some visitors of our Ecological Station asked about the generation of the graphs and so I decided to write down a small description. Is it adviseable to publish the link to the article in an R-wiki (which?) so I can give back a little to the

Hi All, I have been trying to get to plot my PCA legend outside of the PCA plot, but success still alludes me. Can you guys please advise how I can achieve this. I used locater() to obtain coordinates for below the Comp.1 axis. Using these coordinates the legend disappears. Below is the code for the PCA and legend. Thanks in advance for the help. Regards Tinus r.cols <-