similar to: p-value in changepoint package

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,

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.

Hello, I would like to evaluate the relationship between flows and phytoplankton abundance (or Chlorophyll a concentrations) using a changepoint analysis.? Specifically, I have two study questions: Study Question 1: Are there certain flow thresholds that result in spikes in phytoplankton abundance? Study Question 2: Are the duration of certain flows important for phytoplankton abundance (e.g.,

Hello I have lots of folders of individual Scottish river catchments on my uni One Drive. Each folder is labelled with the river name eg "Tay" and they are all in a folder named "Scotland" I want to list the folders on One Drive so that I can cross check that I have them all against a list of folders on my laptop. Can I somehow use list.files() - I've tried various things

? Tue, 9 Jan 2024 16:42:32 +0000 Nick Wray <nickmwray at gmail.com> ?????: > she has a problem with R studio on her laptop Does the problem happen with plain R, without Rstudio? What's the student's sessionInfo()? > I have a screenshot which could email if anyone needs to see what it > looks like. I think that PNG screenshots are allowed on the mailing list, so it could

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 I'd like to create a simple topological map of a river system, ie all the waterways, where they join and points on them such as gauging stations so that it looks something like the London Tube map - where distances are not important but the relationships of the nodes are. I could probably Heath-Robinson something together but I wondered whether there was any kind of package already out

Hello I've successfully installed the package "Detpack" and remembered to call it via library() But when I try to use the functions i get this: chi2testuniform(vals, 0.05) Error in chi2testuniform(vals, 0.05) : could not find function "chi2testuniform" And yet this function (amongst other which also are not recognised) are all listed on the CRAN page:

Hello I am using k medoids in R to generate sets of clusters for datasets through time. I can plot the individual clusters OK but what I cannot find is a way of pulling out the co-ordinates of the individual points in the cluster diagrams - none of the kmed$... info sets seems to be this. Beneath is an example of a k medoid prog using the built in US arrests dataset - this is not the data I am

Hello As a postgrad I have been helping an undergraduate student with R coding but she has a problem with R studio on her laptop which I can't fix - basically when she runs a plot it appears without a y axis label with the black line plot frame hard against the plot window and the bottom of the plot, where you would expect to see the horizontal axis and the x axis label etc is completely

R help email Hello I?m trying to work my way through the rnrfa package RJ-2016-036.pdf (r-project.org) <https://journal.r-project.org/archive/2016/RJ-2016-036/RJ-2016-036.pdf> I?ve got it loaded OK and can access the data by the catalogue() command - gives me a tibble with names, grid refs etc allStations <- catalogue() allStations This works OK But when I try to run the example

I would just use fi <- file.info(dir(path, recursive=TRUE, include.dirs=TRUE)) path could be the OneDrive directory or Scotland (and is not needed if you're already in the directory you want). Then rownames(subset(fi, isdir)) will contain all the directories. Maybe you want to use grep or other machinery to thin it out. Regards, Jorgen Harmse.?

I did use "detpack" ie not with a capital detpack::chi2testuniform(vals,0.05) gives this: Error: 'chi2testuniform' is not an exported object from 'namespace:detpack' ?? Thanks Nick On Wed, 1 Feb 2023 at 16:29, Eric Berger <ericjberger at gmail.com> wrote: > Detpack or detpack? > > What happens when you try detpack::chi2testuniform(...) ? > > >

Hi, I am attaching a plot where you can see there are a few "jumps" (plots 1, 4, 5 and 6), due to incidents with the measuring sensors (basically someone touching the sensor). I need to revert those changes to have a plot without unreal measurements, so make those fragments go back to its original pattern before the jump. I have used the function cpt.mean {changepoints} so I can

Dear friends. I would like to link an object file made in C. Dyn.load rather wants to link a whole library and I only wanted to link this small object file. I followed instructions given for S - by Clive Loader in his software for changepoint handling and it is his object file I made from http://cm.bell-labs.com/cm/ms/departments/sia/clive/smcp.shar Any dirty tricks or easy oversights ? Yours

I am looking for R code to be able to fit a linear-linear piecewise model with person-specific changepoint. I have searched the web, but have not been able to locate any code. Below is my attempt at some code: chgpt = function(a1,a2,a3,gam,wave){ yht=numeric(10) y1=(wave <= gam)*(a1+(a2*wave)) y2=(wave > gam)*((a1+(a2-a3)*gam)+a3*wave) yhat=y1+y2 return(yht) } nl.dat <- nlme(y ~

Mmmh, seems as if I really should change my options as I seem to keep sending off empty bug-reports ;-/ Sorry guys. Here is the content that should have been in the last e-mail: `R CMD check nlme' fails on my machine. The final output in nlme-Ex.Rout is: > library(nlme) > data(Soybean) > fm1 <- nlme(weight ~ SSlogis(Time, Asym, xmid, scal), data = Soybean, + fixed =

Dear R-Help list, I have a nonlinear least squares problem, which involves a changepoint; at the beginning, the outcome y is constant, and after a delay, t0, y follows a biexponential decay. I log-transform the data, to stabilize the error variance. At time t < t0, my model is log(y_i)=log(exp(a0)+exp(b0)) at time t >= t0, the model is log(y_i)=log(exp(a0-a1*(t_i - t0))+exp(b0=b1*(t_i -

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