similar to: Suggests and examples

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

Hello, I'm trying to rebuild a package (using R version 3.3.2 (2016-10-31)) that has not seen any changes since it was last built without problems in previous R versions. It fails during vignette building with: ---<--------------------cut here---------------start------------------->--- $ R CMD build diveMove * checking for file ?diveMove/DESCRIPTION? ... OK * preparing ?diveMove?: *

Hi All: Has anyone analyzed Illumina 550K data using DNAcopy package? I have around 2300 samples. According to Venkatraman and Olshen 2007 paper, it needs about 25 min to run a single sample for Affymetrix 100K. I am afraid it will take too long to analyzed my data. Anybody has an idea? Thanks! Best, -- View this message in context:

Hello R Developers I am using DNACopy package http://bioconductor.org/packages/1.9/bioc/html/DNAcopy.html I am not able to figure out how to (if at all possible) to modify output format, namely, I am getting the following: "ID" "chrom" "loc.start" "loc.end" "num.mark" "seg.mean" Is it a way to get also median and standard deviation?

Hi, I am trying to read my TDR data into R using the readTDR function for the diveMove package. > seal <- readTDR("file location and name here", dateCol=1, depthCol=3, > speed=FALSE, subsamp=1, concurrentCols=4:5) But I keep getting the following error: > Error: all(!is.na(time)) is not TRUE All my columns to have values in them (there are no empty records) The manual

Hi guys, I've got a problem where if I set permissions on a folder(Admin) "setfacl -R -d -m u::rwx,g::wrx,o::- Admin/" and "setfacl -m o:- Admin" I get the following. mail:/data/samba/shared # getfacl Admin/ # file: Admin # owner: BCP+administrator # group: samba user::rwx user:samba:rwx group::rwx group:BCP+admin:r-x mask::rwx other::--- default:user::rwx

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.

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, Doing 'R CMD check diveMove' is now throwing this message: Data codoc mismatches from documentation object 'sealLocs': Variables in data frame 'sealLocs' Code: id.time.class.lon.lat Docs: class id lat lon time with: R> sessionInfo() R version 2.10.1 (2009-12-14) x86_64-pc-linux-gnu locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C

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, I'm complety new to "R" and have a question. I've read through a couple of manuals but I'm having a problem with getting something run properly. I'd like to attempt to use the "strucchange" package with some sample data however I'm having trouble understanding the proper syntax of the commands from which to do so. I basically want to take

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

Hi All, I am using DNAcopy package in R for copy number analysis of 500K chip. The final output which I get from DNA copy is too big to be printed in a file. So I am getting an error as "reached getOption("max.print") -- omitted 475569 rows " Can somebody please provide me the pointers with how to increase the limit for max.print . Thanks, Pooja DISCLAIMER\

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

Hi folks, I'm trying to estimate bias-corrected percentile (BCP) confidence intervals on a vector from a simple for loop used for resampling. I am attempting to follow steps in Manly, B. 1998. Randomization, bootstrap and monte carlo methods in biology. 2nd edition., p. 48. PDF of the approach/steps should be available here: https://wyocoopunit.box.com/s/9vm4vgmbx5h7um809bvg6u7wr392v6i9 If

OvermindDL1 <overminddl1 at gmail.com> writes: > On Tue, Aug 17, 2010 at 5:25 PM, Chris Lattner <clattner at apple.com> wrote: >> On Aug 17, 2010, at 3:42 PM, OvermindDL1 wrote: >> >>> >>> On Tue, Aug 17, 2010 at 4:09 PM, David A. Greene <greened at obbligato.org> wrote: >>>> OvermindDL1 <overminddl1 at gmail.com> writes:

>> >> >> Hi All, >> >> After getting licences for Skype for asterisk a while ago I finally got >> around to setting up a server with two channels and setting up a bcp on >> the skype end. >> >> The idea behind this is the following: >> >> Users can dial into the PBX, get authenticated and only after> >>

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

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))