search for: supremum

Hi Given a real number x, I want to know how accurately R can represent numbers near x. In particular, I want to know the infinum of exactly representable numbers greater than x, and the supremum of exactly representable numbers less than x. And then the interesting thing is the difference between these two. I have a little function that does some of this: f <- function(x,FAC=1.1){ delta <- x while(x+delta > x){ delta <- delta/FAC } return(delta*FAC) } But this c...

...package truncgof and run the example for function ad.test. I got the following output: set.seed(123) treshold <- 10 xc <- rlnorm(100, 2, 2) # complete sample xt <- xc[xc >= treshold] # left truncated sample ad.test(xt, "plnorm", list(meanlog = 2, sdlog = 2), H = 10) Supremum Class Anderson-Darling Test data: xt AD = 3.124, p-value = 0.12 alternative hypothesis: two.sided treshold = 10, simulations: 100 So I cannot reject the hipothesis (at a standard confidence level) that the original sample comes from a lognormal distribution (as it is the case). But let us t...

Hi, I'm trying to understand the complete linkage method in hclust. Can anyone provide a breakdown of the formula (p9 of the pdf documentation) or tell me what the "sup" operator does/means? thanks in advance Tom [[alternative HTML version deleted]]

...o test if my dataset is a lognormal distribution or not (Histogram shows a nice normal distribution in log scale but I want to check) 2) two subsets from this dataset have same mean or not (like "t test" of normal distribution) What I tried are 1) ad.test{truncgof} of R, which is a "Supremum Class Anderson-Darling test" but I got error message. What I had problems are 1) Memory error happen with error message of "cannot allocate vector of size 1.0GB" and I cannot increase the computer spec. 2) I don't know what is "t test" of lognormal. What I am thinking...

Dear all, I use the Andersen plot to check for proportional hazards assumption for a factor (say x) in the Cox regression model and obtained a straight line that pass through the origin. However, the formal test done by the R-function cox.zph, which is based on the plot of Schonefeld residuals against time, indicates that proportional hazards assumption is violated. Further, a plot of the score

...ad.test. I > got the following output: > > set.seed(123) > treshold <- 10 > xc <- rlnorm(100, 2, 2) # complete sample > xt <- xc[xc >= treshold] # left truncated sample > ad.test(xt, "plnorm", list(meanlog = 2, sdlog = 2), H = 10) > > > Supremum Class Anderson-Darling Test > > data: xt > AD = 3.124, p-value = 0.12 > alternative hypothesis: two.sided > > treshold = 10, simulations: 100 > > > So I cannot reject the hipothesis (at a standard confidence level) that > the original sample comes from a lognormal dis...

Hi, I seek to identify those points in/outside a multidimensional convex hull (geometry::convhulln). Any suggestions? Background just in case I'm going down a really wrong road: Given an observed data set with one dependent/observed variable (Y) and multiple (3 to 10) independent/design variables (X1, X2, ...) I want to increase the number of points by interpolating. I'm using

...? Competing risks Model Test for nonparametric terms Test for non-significant effects sup| hat B(t)/SD(t) | p-value H_0: B(t)=0 (Intercept) 0 0 random 0 0 Test for time invariant effects supremum test p-value H_0: B(t)=b t (Intercept) 0 0 random 0 0 int (b(t)-g(t,gam))^2dt p-value H_0:constant effect (Intercept) 0 0 random 0...

CRAN (and crantastic) updates this week New packages ------------ * allan (1.0) Alan Lee http://crantastic.org/packages/allan Automates Large Linear Analysis Model Fitting * andrews (1.0) Jaroslav Myslivec http://crantastic.org/packages/andrews Andrews curves for visualization of multidimensional data * anesrake (0.3) Josh Pasek http://crantastic.org/packages/anesrake This