similar to: running strucchange?

hello R-world, according to the strucchange package .pdf, "all procedures in this package are concerned with testing or assessing deviations from stability in the classical linear regression model." i'd like to test/assess deviations from stability in the Poisson model. is there a way to modify the strucchange package to suit my purposes, or should i use be using another

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:

Greetings - Am hoping to use the strucchange package to look for structural breaks in some messy regression data. A series of preliminary analyses indicate that BLUE for these data will involve some weighting the data (estimates of a particular population parameter) by a function of the variance of the estimate (say, inverse of the variance). While I've gone through the docs for

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

When I run the code below on Macintosh and Windows, the plot comes out fine. However, on Linux, the png generated is invalid from R console, and loading strucchange crashes rkward. Is this a known issue on Linux and, if so, is there a workaround? Many thanks! require(strucchange) data("RealInt") bp.ri <- breakpoints(RealInt~1, h=15) summary(bp.ri) fac.ri <- breakfactor(bp.ri,

Hello, Im using strucchange package in R software in order to apply Bai and Peron (1998, 2003) structural break tests to a set of n=1671 observations with a constant term (no AR terms). For that purpose I have read several papers, for instance Validating Multiple Structural Change Models An Extended Case Study, in which its aim is to replicate the results from Bai and Perron (2003) in R

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

Hi Achim, I'd like to try to run the strucchange package on R. However, it seems that the package can only run on systems running Debian and not Windows XP. Is this true? Or is there a windows version? I downloaded the strucchange file but they seem to have the dot deb extension and seem to be designed to run with Debian. If there is a windows version, could you please let me know how to

Hello, I am just studying the following example from vignette: strucchange-intro, contineousely ending up in an error. This is the given code: 1. library(strucchange) 2. data(USIncExp) 3. if (!"package:stats" %in% search()) library(ts) 4. USIncExp2 <- window(USIncExp, start = c(1985, 12)) A.Modelling: coint.res <- residuals(lm(expenditure ~ income, data = USIncExp2))

I want to apply Nyblom-Hansen test with the strucchange package, but I don't know how is the correct way and what is the difference between the following two approaches (leeding to different results): data("longley") # 1. Approach: sctest(Employed ~ Year + GNP.deflator + GNP + Armed.Forces, data = longley, type = "Nyblom-Hansen") #results in: # Score-based CUSUM

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

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

Hello, to understand better how efpFunctional works, I'm trying to construct my own functionals. But concerning already existing functionals I have some questions. With maxBB it is clear: functional = list(comp = function(x) max(abs(x)), time = max), with rangeBB: functional = list(time = function(x) max(x)-min(x), comp = max), with meanL2BB, if I understood correctly: functional =

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

Hi, I am using the Strucchange package in R to test for structural change in regression coeffcient. Given a model y = b0 + b1*X + b2*Z, the Fstats test whether there is a change in both b1 and b2 over a time period. Is there any way where I can restrict the test to hold b2 constant and test for break in only b1? That is, instead of a pure structural change, could I test for partial structural

Hi, I am a new user to R. I am using strucchange to generate breakpoints: -------------------------------------------------------------------------------------- > res <- gbreakpoints(GDP.new ~ 1,data=a,h=2,breaks=5) > print(res) Optimal 6-segment partition for `lm' fit: Call: gbreakpoints(formula = GDP.new ~ 1, data = a, h = 2, breaks = 5) Breakpoints at observation number:

It's two days I try to answer to a mail of Achim Zeileis, the maintainer of strucchange package, using the mail address: zeileis at ci.tuwien.ac.at. I got the following error: ----- The following addresses had permanent fatal errors ----- <zeileis at ci.tuwien.ac.at> (reason: 550 5.7.1 <zeileis at ci.tuwien.ac.at>... Access denied) Can anyone tell me how to reach him?

Hola! I am using package strucchange, and encounters the following: > bp <- breakpoints(diesel90 ~ regress -1, h=NULL) Error in La.chol2inv(x, size) : element (14, 14) is zero, so the inverse cannot be computed The obvious problems have been checked, that is, the model matrix is of full rank. What can be causing this? I can send some data if that can be of help. Kjetil Halvorsen

Hello R users: I'm trying now to apply the package strucchange to see whether there is a structural change in linear regression. I have noted the following problem that arises in my case with recursive-based CUSUM: generic function recresid() in efp() generates an error, since (probably) it cannot compute the inverse matrix of (X^(i-1)^T)*(X^(i-1)) at each step (i-1), because the matrix

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