search for: unimod

Hi, Why is the following retuning a nodset of length 0: > library(XML) > test <- xmlTreeParse( > "http://www.unimod.org/xml/unimod_tables.xml",useInternalNodes=TRUE) > getNodeSet(test,"//modifications_row") Thanks for any hint. Joh

Dear All, A method for the estimation is univariate unimodal densities (with unknown mode) is described in "Statistical Inference under Order Restrictions" by Barlow et al.. Would anyone know whether there is an R-implementation (preferably with reference) for the estimation of univariate discrete unimodal densities (with unknown mode)? Thank...

Hola a tod en s, estoy intentando averiguar el número de modas en una distribución. Para ello utilizo diptest. Mi duda es que no acabo de entender cuando la información suministrada por los test suponen la existencia o no de unimodalidad/multimodalidad. Una parte de la salidad de diptest es la que pego a continuación (el resto esta en el fichero adjunto con las distribuciones kernels y las soluciones gráficas). *Distribución #1(suponía que hay más de una moda)* Hartigans' dip test for unimodality data: Datos8$V2...

Dear All, Does anybody know if there is a functionality in R to break histograms that show a clear bi-modal (or multi-modal) distribution into a series of unimodal histograms that added up result in the original histogram? I was thinking of using QQ-plots (for which tools are available in R), and then observing the number of times the observed quantiles cross the 1:1 line, but this only gives an indication of how many "peaks" the current histogra...

Dear all i could not estimate the optima value or range value in unimodal plot in glm please help me out thanking you regard madan _________________________________________________________________ Hotmail® has a new way to see what's up with your friends. orial_WhatsNew1_052009 [[alternative HTML version deleted]]

Dear R users, I am interested in uni- bi- multimodality tests, for analysing reaction times data. I was lead to Hartigan's dip test (Ann. Statistics, 13, 1985, pp. 70-84, Applied Statistics, 34, 1985, 320-325). Not being a programmer I am unable to translate the Fortran code given in ref. 2 into a R function. I'd be glad to learn that someone already did it, or has devised a better

...SpamAssassin 2.63 (2004-01-11) on hypatia.math.ethz.ch X-Spam-Level: **** X-Spam-Status: No, hits=4.2 required=5.0 tests=MSGID_FROM_MTA_HEADER,RCVD_IN_BL_SPAMCOP_NET autolearn=no version=2.63 Dear r-help list, my question is about ordination technics: To decide using a linear model response or a unimodal model response, in Canoco (software for ordination technics) you can calculate a so-called "length of gradient". (from Canoco Help file:) > The length of gradient is a measure of how unimodal the species > responses are along an ordination axis. It is the range of the > sample...

Hi all, I am using Hartigan & Hartigan's [1] "dip test" of unimodality via the diptest package in R. The function dip() returns the value of the test statistic but I am having problems calculating the p-value associated with that value. I'm hoping someone here is familiar with this process and can explain it. In the original article there is an example u...

Has a test for bimodality been implemented in R? Thanks, John NIWA is the trading name of the National Institute of Water & Atmospheric Research Ltd.

Dear List, I have a set of distributions recorded at an equal interval of time and I would like to plot them as series of horizontal histograms (with the x-axis representing time, and y-axis representing the bins) since the distribution shifts from unimodal to multimodal in several occasions. What I would like to see is something close to a violinplot, but I do not want a kernel density estimate... Any suggestions or advice will be great! Thanks in Advance, Enrico

I'm not sure where to begin with this, but I was wondering if someone could refer me to an R package that would test to see if a distribution fits a bimodal distribution better than a unimodal distribution. Thanks, Andrew [[alternative HTML version deleted]]

Hi, I am trying to use the optimize function to optimize a function. The results I am getting don't agree with what I compute on my own and when I look at the graph of f(x) = 100 + ((x-10)**2 + (x-10)) * cos(x-10), where -10 <= x <= 10 in gnuplot. I suspect I am making a mistake in the usage of the R optimize function, perhaps someone could point out where? >

Dear all, I new to r and I would like your help. I want to explore the patterns (unimodal, monotonically increased/decreased) of species richness~altitude using GAM in R. Although I run the gam function in mgcv package I do not know how to manually define knots and degrees of freedom. Any help would be greatly appreciated. Spyros -- View this message in context: http://r.789695.n...

Hi all, I have the data from the posterior distribution for some parameter. I want to find the 95% credible interval. I think "t.test(data)" is only for the confidence interval. I did not fine function for the Bayesian credible interval. Could some one suggest me? Thanks [[alternative HTML version deleted]]

...xa taken 1:11 at a time is 2048, all combos of 12 taken 1:12 is 4096, 13 taken 1:13 = 8192...etc etc such that when I reach about 25 taken 1:25 the number of combos is 33554432 and I get a memory error. I have found that the number of combos of x taxa taken from a pool of n creates a very kurtotic unimodal distribution,... x<-vector("integer",20) for (i in 1:20) {x[i]<-choose(20,i)} plot(x) ...but have found that limiting the number of samples for any community size to 1000 is good enough for the further analyses I wish to do. My problem lies in sampling all possible combos witho...

...ocess_cpu_clocks.o ./boost/chrono/src/thread_clock.o ./pwiz/data/msdata/Version.o ./pwiz/data/common/MemoryIndex.o ./pwiz/data/common/CVTranslator.o ./pwiz/data/common/cv.o ./pwiz/data/common/ParamTypes.o ./pwiz/data/common/BinaryIndexStream.o ./pwiz/data/common/diff_std.o ./pwiz/data/common/Unimod.o ./pwiz/data/msdata/SpectrumList_MGF.o ./pwiz/data/msdata/DefaultReaderList.o ./pwiz/data/msdata/ChromatogramList_mzML.o ./pwiz/data/msdata/examples.o ./pwiz/data/msdata/Serializer_mzML.o ./pwiz/data/msdata/Serializer_MSn.o ./pwiz/data/msdata/Reader.o ./pwiz/data/msdata/Serializer_MGF.o ....

...e, Sep 15, 2009 at 9:01 AM, David Winsemius <dwinsemius at comcast.net >>>> wrote: >>>> I get very different output from the two versions of Mac OSX R as >>>> well. The 32 bit version puts out a histogram that has an expected, >>>> almost symmetric unimodal distribution. The 64 bit version >>>> created a >>>> bimodal distribution with one large mode near 0 and another smaller >>>> mode near 10E+37. Postcript output attached. >>>> >>> >>> >>> >>> -- >>> Ani...

...for any input Rob Knell School of Biological Sciences Queen Mary, University of London 'Phone +44 (0)20 7882 7720 http://www.qmw.ac.uk/~ugbt794 http://www.mopane.org Here's the code: switchline<- function(X,Y,Intmin,Intmax,Intgap,Slopemin,Slopemax,Slopegap) { # X = unimodal variable (e.g. filename$bodysize) # Y = bimodal variable (e.g. filename$hornsize) # Intmin = Minimum switchline intercept # Slopemin = Minimum switchline slope # Intmax = Maximum switchline intercept # Slopemax = Maximum switchline slope # Intgap...

Dear R-users Do you know of an easy way in R of performing "isotonic (unimodal)" regression ala BBBB, Barlow, Bartolomew, Brenner and Brunk. best regards Helgi -- Helgi Tomasson FAX: 354-552-6806 University of Iceland PHONE:354-525-4571 Faculty of Economics and Business Administration email:h...

Hi, anyone knows about any functions in R can get multidimensional integration not over a multidimensional rectangle (not adapt). For example, I tried the following function f(x,n)=x^n/n! phi.fun<-function(x,n) { if (n==1) { x }else{ integrate(phi.fun, lower=0, upper=x, n=n-1)$value } } I could get f(4,2)=4^2/2!=8, but failed in f(4,3)=4^3/3! Thanks Best, Lynette