similar to: Strucchange gbreakpoints

Hi, I am new to R. I am using strucchange to get the breakpoints in time series dataset. So the problem I am facing is: I want to link the result generated by the breakpoints to further analysis (for eg. generating volatility for each group). The result is in following form: --------------------------------------- > res <- gbreakpoints(GDP.new ~ 1,data=a,h=2,breaks=6) > res

Dear R users, I'm trying to use the "strucchange" package to determine structural breaks in an ARFIMA model. Unfortunately I'm not so familiar with this topic (and worse, I'm a beginner in R), so I don't know exactly how to specify my model so that the "Fstats","sctest" and "breakpoint" functions to recognize it and to calculate the

Hi, I am new to R. It would be kind if I could get some help on this. I am using R to estimate Fstats but I am getting following error. a3 is annual GDP data from 1951 to 2010. > fs<- Fstats(ecm.model, from=1954, to = 1975,data=a3) Error in Fstats(ecm.model, from = 1954, to = 1975, data = a3) : inadmissable change points: 'from' is larger than 'to' In addition: Warning

Respected Sir I tried the strucchange My data is attached. However I tried the attached commands (last save.txt) to perform Bai Perron 2003... I t worked well but in the end it is giving warning that overlapping confidence interval... I am not sure how to proceed... Please Help Me Thanking You Ayanendu Sanyal -- Please have a look at our new mission and contribute into it (cut and paste the

Is it possible to use some form of robust regression with the breakpoints routine so that it is less sensitive to outliers? --Rich Richard Kittler Advanced Micro Devices, Inc. Sunnyvale, CA

Dear all, I'm trying to understand how the strucchange package is working and I have been looking at the examples given for the Fstats() function. The first example (Nile), shows one peak in the F-stats and one breakpoint is estimated, that can be plotted using the following code ## Nile data with one breakpoint: the annual flows drop in 1898 ## because the first Ashwan dam was built

Folks, I have a strange situation, which I may have isolated as a bug report. Or, it could just be that there's something about R that I don't know. :-) I have attached the data file and the program file but don't know whether these attachments will make it into the list. Here is my bugreport.R program -- ---------------------------------------------------------------------------

I am attempting to manually specify x axis value labels (NOT x axis label as in xlab) on a density plot. Generally speaking, I would like to erase the default x axis point labels and replace them with arbitrary x axis point labels relating to specific CDF points. I am sure this can be done, but it isn't obvious to me how. Has anyone done this before? My R code follows if it helps (new

HI i'm trying to detect breaks points in various flow time series, they all contains seasonality and trend my question is : i have to remove this seasonality and trend before apply the function breakpoints du package strucchange?? another question, the function breakpoints is similar to de Pettit tests ? or how does it realy works? THANKS!!!! DENISSE -- View this message in context:

If I try the breakpoints() function (strucchange package) with a minimum segment size = the number of regressors, there appears the following error message: "minimum segment size must be greater than the number of regressors" According to the documentation: "breakpoints implements the algorithm described in Bai & Perron (2003) for simultaneous estimation of multiple

Hello, we want to identify breakpoints (different phases) in environmental data, algae cell counts of three years with intervals between 7 and 30 days (N=40). We found that breakpoints(cells ~1) works great and identifies 5 very good breaks, however we are uncertain about these, because the data are unequally spaced. Is there a way to include the information about the measurement intervals,

Dear all, I'm trying to remove some text after the period (a decimal point) in the data frame 'hi', below. This is one step in formatting a table. So I would like e.g. "2.0" to become "2" and "5.3" to be "5.3", where the variable digordered contains the number of digits after the decimal that I would like to display, in the same order in which

I don't understand this behavior. Why does the every data point get trashed by data.matrix when there is one non-numeric element in the array? Thanks. > temp GDP CPIYOY 19540 2098.1 garbage 19632 2085.4 0.9 19724 2052.5 0.8 19814 2042.4 1.1 > data.matrix(temp) GDP CPIYOY 19540 4 4 19632 3 2 19724 2 1 19814 1

Dear R-Team, Am I right supposing that the "breakpoints()" function in the strucchange package is an implementation of the pure structural change model proposed by Bai and Perron (1997, 2003)? My question relates to a partial structural change model that Bai and Perron formulate in their 2003 paper, e.g. formulated as y = x' beta + z' delta_j + epsilon, where beta and delta

Hi! I'm trying to do a lm() test on three objects. My problem is that R protests and says that the variable lengths differ for one of the objects (Sweden.GDP.gap). But I have double checked that the number of observations are the same. All three objects should contain 9 observations but R only accepts 9 observations in two of the objects. The third must have 10! Very confusing because there

Hi to all, I am trying to run breakpoints() on a fairly large sample (>10.000 observations). The process is very slow, any idea on how to speed this up? I have tried the hpc="foreach" parameter, but this didn't work at all when I tried to run it on a smaller sample. breakpoints(x ~ x.l1 + x.l2 + X.l3 + x.l4 + x.l5 + x.l6 + x.l7 + x.l8 + y.l1 + y.l2 + y.l3 + y.l4 + y.l5 + y.l6

Can anyone please tell me how to get the r-squared output from a piecewise (segmented) regression using the strucchange package? Here is the R code I have tried thus far. library(lmtest) library(strucchange) data <- ts(c(rnorm(30), runif(30)), frequency = 12, start = c(2005, 01)) bpts <- breakpoints(data ~ 1) print(bpts) summary(bpts) coeftest(bpts) [[alternative HTML version

Dear R Users, I am attempting to use the strucchange package but get an error which is difficult to decipher. I get the following error with breakpoints: > bp<-breakpoints(regr[,1]~regr[,2]+regr[,3]) Error in my.RSS.table[as.character(i), 3:4] <- c(pot.index[opt], break.RSS[opt]) : nothing to replace with I was wondering if anyone has seen this issue before and off the top of

Good morning to all, I am encountering a blocking issue when using the function 'breackpoints' from package 'strucchange'. *Context:* I use a data frame, 248 observations of 5 variables, no NA. I compute a linear model, as y~x1+...+x4 x4 is a dummy variable (0 or 1). I want to check this model for structural changes. *Process & issues:* *First, I used function Fstats.* It

Dear all, I have the following question. I have to fit the hierarchical model for the hypothesis concern the individual-level effects by controlling for the individual -level attributes and national-level contextual effects on individuals by using R. O have to obtain the estimates of the impact of the second-level (national: GDP per capita) effects on individuals ( in this instance the impact