similar to: R strucchange question -- robust regression

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:

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

Inspired by colouring the dots of box-whisker plots I am trying to also fill the boxes (rectangles) with different colours. This seems not to work as I expected. Looking at the help page of panel.bwplot it says: 'fill - color to fill the boxplot'. Obviously it is only intended to fill all boxes with only one colour? Nevertheless the following example shows, that 'fill' from

Wed, Apr 18, 2018 at 08:43:15PM CEST, sridhar.samudrala at intel.com wrote: >On 4/18/2018 2:25 AM, Jiri Pirko wrote: >> Wed, Apr 11, 2018 at 09:13:52PM CEST, sridhar.samudrala at intel.com wrote: >> > On 4/11/2018 8:51 AM, Jiri Pirko wrote: >> > > Tue, Apr 10, 2018 at 08:59:48PM CEST, sridhar.samudrala at intel.com wrote: >> > > > This provides a

Wed, Apr 18, 2018 at 08:43:15PM CEST, sridhar.samudrala at intel.com wrote: >On 4/18/2018 2:25 AM, Jiri Pirko wrote: >> Wed, Apr 11, 2018 at 09:13:52PM CEST, sridhar.samudrala at intel.com wrote: >> > On 4/11/2018 8:51 AM, Jiri Pirko wrote: >> > > Tue, Apr 10, 2018 at 08:59:48PM CEST, sridhar.samudrala at intel.com wrote: >> > > > This provides a

Hello All, I am working on graph object using IGRAPH package wanted to do Bonacich Power. This is my graph object. The file 'Graph.RData' (4.2 MB) is available for download at http://dropbox.unl.edu/uploads/20090424/cfe4fcb854bb17f2/Graph.RData Graph size Vertices: 20984 Edges: 326033 Directed: FALSE No graph attributes. Vertex attributes: name. No edge attributes. When I use

Wed, Apr 11, 2018 at 09:13:52PM CEST, sridhar.samudrala at intel.com wrote: >On 4/11/2018 8:51 AM, Jiri Pirko wrote: >> Tue, Apr 10, 2018 at 08:59:48PM CEST, sridhar.samudrala at intel.com wrote: >> > This provides a generic interface for paravirtual drivers to listen >> > for netdev register/unregister/link change events from pci ethernet >> > devices with the

Wed, Apr 11, 2018 at 09:13:52PM CEST, sridhar.samudrala at intel.com wrote: >On 4/11/2018 8:51 AM, Jiri Pirko wrote: >> Tue, Apr 10, 2018 at 08:59:48PM CEST, sridhar.samudrala at intel.com wrote: >> > This provides a generic interface for paravirtual drivers to listen >> > for netdev register/unregister/link change events from pci ethernet >> > devices with the

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

I am using 'read.csv' in V2.0.1 to read in a CSV file with the colClasses option and am getting an error from 'scan' when it encounters a non-numeric value for a 'numeric' column, i.e. > ds <- read.csv(in_file, nrows=irow, row.names=NULL, colClasses=zclass, comment.char="") Error in scan(file = file, what = what, sep = sep, quote =

I need to make plots similar to those produced by the s-plus rayplot function but can't seem to find it in R. These 'vector maps' plot a ray or vector at each specified location. Is there something similar in R ? --Rich Richard Kittler AMD TDG 408-749-4099

Does the 'groups' option on qqmath just color the points differently in the main distribution or does it actually overlay separate quantile plots for each subset? I would like to be able to do the latter. --Rich Richard Kittler AMD TDG 408-749-4099

I am trying to coerce a data frame column from character to datetime using strptime but keep getting an error because the length of the coerced object is always 9. What am I doing wrong here: ................................................................. > ds <- cbind(1:2, c("02/27/92 23:03:20", "02/27/92 22:29:56")); ds [,1] [,2] [1,]

Is there a way of controlling the font size and alignment of a footnote in an xyplot, or alternatively of using 'mtext' to place a footnote at the bottom of a graph? --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

I am stuck on trying to get the Lattice xyplot to output a separate PNG file each time through my 'for' loop. The files get produced but are empty. Here is the code. I'm running 1.9 on Windows. BTW is there a more efficient way of creating the separate output files than looping over the levels and subsetting? ......................................................... lev <-

Does anyone know where I can get R code for plotting the Brown , Durbin and Evans cumsum procedure ( 1975 ) ? I wrote my own code but I am a little worried that my confiodence bands may not be correct ( I find the formula in the original paper confusing and S+Finmetrics has a formula but that formula implies that there should be 4 lines as far as I can tell ) so I would like to see someone

Hi, does anybody know of a nice test to detect trend turning points in time series? Possibly with reference? Thanks, joerg

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24),

Dear list, I'm migrating a project from Matlab to R, and I'm facing a relatively complicated problem for nls. My objective function is below: >> objFun <- function(yEx,xEx,tEx,gamma,theta,kappa){ yTh <- pdfDY(xEx,tEx,gamma,theta,kappa) sum(log(yEx/yTh)^2) } The equation is yTh=P(xEx,tEx) + noise. I collect my data in: >> data <-