similar to: Automatically fix big jumps in one variable due to anomalies

Hello Does anyone know how to get the p values for changepoints in the "changepoint" package? There is guidance in the CRAN vignette https://cran.r-project.org/web/packages/changepoint/changepoint.pdf but I can't make it work to get the p value out... Thanks Nick Wray [[alternative HTML version deleted]]

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! How can I cope with overlapping or covered labels (covered by labels from other data points) in plots? Martina Renninger [[alternative HTML version deleted]]

Dear R-help Members, I am wondering if anyone think of the optimal way of computing for several numeric variable the difference between 2 levels of a factor. To be clear let's generate a simple data frame with 2 numeric variables collected for different subjects (ID) and 2 levels of a TIME factor (time of evaluation)

Dear list member, My question is related to input file format to an Anova from car package. Here is an example of what I did: My file format is like this (and I dislike the idea that I will need to recode it): Hormone day Block Treatment Plant Diameter High N.Leaves SH 23 1 1 1 3.19 25.3 2 SH 23 1 1 2 3.42 5.5 1 SH 23 1 2 1 2.19 5.2 2 SH 23 1 2 2 2.17 7.6 2 CH 23 1 1 1 3.64 6.5 2 CH 23 1 1 2

Hi, I have a function that generates a set of data but I am having problems determining the parameters using the nls fitting procedure. #### "MH"<-function(field,diameter,mu=10e-7,sig=0.1,Ms=100,chi=0){ #variables mu, sig, chi, Ms #input: field and diameter #all in CGS rho <- 5 kb <- 1.38e-16 t <- 300 length.d<-length(diameter) length.H<-length(field)

I am trying to define groupings from levels of factor variables and this the warning message that R give "& not meaningful for factors". The nature of my task is this. I have a variable stage which has the levels (1B, 2A, 2B) - these are the AJCC TNM stages of cancer, and another variable diameter with factor levels ("=< 4", "4 - 6.5, > 6.5; limit values are

Hi, Maybe this is more a programming questions than a specific R-project question, but maybe there is someone who can point me in the right direction. I have a picture of cirkels which I took with a digital camera. Now I want to use the diameter of the cirkels on the picture for analysis in R. I can use pixmap to import the picture, but how do I find the outside cirkels and calculate the

Dear all, I am making a barplot as following: barplot(c(1,2,3,5,2,3,1),names.arg=c("100","200","300","400","500","600","700"),xlab="diameter",ylab="flow",main="some title",space=0.1) I am also trying to add a probability density curve, however using lines(density(c(1,2,3,5,2,3,1))) does not give

Hi all, I got one problem with comparing strings like if any string is like "*RIGHT, EPICARDIUM: FOCUS, GRAY-WHITE, SINGLE, APPROX 0.6 CM IN DIAMETER*." and i have to compare "*GRAY-WHITE*" with the above string or otherwise i have to compare " *TUMOR BENIGN*" this string with "*MEDULLRY TUMOR BENIGN,TYP PHEOCHROMOCYTOMA*" i

Dear list I'm trying to solve the following problem since 2 days with no success. Could someone help a newbie, please ? I have a dependant variable which is diameter of shell, and two factors which are level on the beach (3 levels 'B', 'E' and 'H') and a hydrodynamics indice (3 levels 1, 2 and 3) To study diameter variations according to these factors, I may use

Hi all, I would like to test the homoegeneity of variances between several linear model for some analysis of covariance. It seems that the Bartlett's test is a good test to use but I am having problem using with linear model and I cannot find any examples on the internet. There are some examples for comparisons of variances but not linear models. If I take the hellung data set, which is

I have a data set with the following structure (with many more obs.): var1 expt day diameter 1 1 2 0.5 1 1 3 0.9 1 1 4 1.3 1 1 5 1.7 1 2 2 0.3 1 2 3 0.5 1 2 4 0.9 1 2 5 1.6 2 1 2 0.7 2 1 3 1.2 2 1 4 1.6 2 1 5 2.3 2 2 2 0.4 2 2 3 0.8 2 2 4 1.6 2 2 5 3.2 I can get separate regression analysis for each level of var1 and expt with the command: by(data3, data3$var1:data3$expt, function(x)

Dear all, I have a function to optimize for a set of parameters and want to set a constraint on only one parameter. Here is my function. What I want to do is estimate the parameters of a bivariate normal distribution where the correlation has to be between -1 and 1. Would you please advise how to revise it? ex=function(s,prob,theta1,theta,xa,xb,xc,xd,t,delta) { expo1=

Data into R from Excel csv file xd<-read.csv("court.dates1.txt",as.is=T, header = F) > str(xd) 'data.frame': 5 obs. of 1 variable: $ V1: chr "6/6" "5/27" "5/16" "5/2" ... >xd V1 1 6/6 2 5/27 3 5/16 4 5/2 5 4/29 cdates <- as.Date(xd, format = " %m/ %d") Error in as.Date.default(xd, format = " %m/

Dear all, The 'make check' step fails for the pacakge mva on IBM AIX. The tail of the Rout log file looks like: > for(factors in 2:4) print(update(Harman23.FA, factors = factors)) Call: factanal(factors = factors, covmat = Harman23.cor) Uniquenesses: height arm.span forearm lower.leg weight 0.170 0.107 0.166

I'm experiencing some confusion with the ellipsis argument (...). In my .Rprofile, I have the following functions: stderr <- function(x, ...) { sqrt( var(x, ...) / length(x) ) } se <- stderr I can use tapply to calculate some means: > tapply( Diameter, factor(Region), mean, na.rm=TRUE ) I II III IV V 0.02896429

Dear R Development Team, I compile R-1.9.1 on AIX 5.2 under 2.9-aix51-020209, and xlf 7.1. In order to let R compile successfully under gcc 2.9, I have to change one C statement of file(RHOME//src/modules/X11/devX11.c) line 1768 from "retrun FALSE" to "retrun NULL", following is C code snapshot: newX11Desc * Rf_allocNewX11DeviceDesc(double ps) { newX11Desc *xd;

Hello, I want to fit a Gompertz model for tree diameter growth that depends on a 4 levels edaphic factor (?Drain?) and I don?t manage to introduce the factor variable in the formula. Dinc is the annual diameter increment and D is the Diameter. >treestab > Dinc D Drain [1,] 0.03 26.10 2 [2,] 0.04 13.05 1 [3,] 0.00 24.83 1 [4,] 0.00 15.92 4