search for: aitchison

...proach (based on EPA’s guidance found in http://www.epa.gov/oswer/riskassessment/pdf/ucl.pdf). This looks like it should be straight forward, but I am need to calculate an MVUE for the population mean and an MVUE for the population variance, which requires a value (g_n) from a table A7, found in Aitchison and Brown (1969): The lognormal distribution. I have looked across the RSiteSearch and can not seem to find a function that will give me g_n or the MVUE for mean and variance of lognormal distribution. Is there an R function that will give me g_n or will calculate an MVUE for the population mean a...

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 2490 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010301/996dcd36/attachment.bin

The compositional data xi=(x_i1,x_i2,...,x_in), for each fixed i , xij>0, and sum(xij)=1; I want to compare the mean( u_i) of several groups i.e. H0: u_1=u_2=...=u_N or H0: u_11=u_21=...=u_N1 Are there any ANOVA tpye tools to do this work in R? Thanks, WEN S Q [[alternative HTML version deleted]]

...--------------------- for anyone interesed in SOM I found Samuel Kaski's thesis "Data exploration using self-organizing maps" fairly accessible http://www.cis.hut.fi/~sami/thesis.ps and there is a wealth of material at http://www.cis.hut.fi/research/som-research/som.shtml John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

The (slow) SOM code can be found on my web page (towards the bottom) http://peabody.vanderbilt.edu/depts/psych_and_hd/faculty/wallern/ On Thu, 1 Mar 2001, John Aitchison wrote: > > Is there an implementation of the SOM ('Self Organizing Map') procedure in > R ? > > I am aware of the implementations of Sammon mapping, multidimensional > scaling and, somewhat peripherally, principal curves and projection > pursuit .. but not SOM per se...

...; inc function(a){ eval.parent(a<-a+1,1) } > z<-1 > inc(z) [1] 2 > z [1] 1 > no good Is what I want to do possible? (I realize that it is possibly inadvisable as there would be no obvious way of knowing from the function definition that it modified its arguments) thanks John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

...through the slideshow at http://www.dmg.org/ but wonder if anyone else has examined PMML and can give an opinion as to whether this (PMML) is worthwhile opening up as a discussion thread. If not, I will look at it in a day or so when I can get some energy and time and report back. thanks John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

...hat rnorm(n) was what was intended, not that it would make any difference in practice in this case. Is this a standard practice? .. that in passing in a vector where a scalar is expected it is as if you passed in the length of that vector as the argument (unless the vector has length 1) John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

...ot;i386-pc-mingw32" $arch [1] "x86" $os [1] "Win32" $system [1] "x86, Win32" $status [1] "" $major [1] "1" $minor [1] "2.1" $year [1] "2001" $month [1] "01" $day [1] "15" $language [1] "R" John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

...mingw32" $arch [1] "x86" $os [1] "Win32" $system [1] "x86, Win32" $status [1] "" $major [1] "1" $minor [1] "2.1" $year [1] "2001" $month [1] "01" $day [1] "15" $language [1] "R" John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

...m_slope Bandwidth(s): 29.22422 2.500559e-24 3.111467 class_unsup_pc_iso Bandwidth(s): 0.2304616 Bandwidth Type: Fixed Log Likelihood: 1531.598 Continuous Kernel Type: Second-Order Gaussian No. Continuous Vars.: 3 Unordered Categorical Kernel Type: Aitchison and Aitken No. Unordered Categorical Vars.: 1 > predict(npu_dens,newdata=origin[training_cols_select])) [1] 7.762187e+18 9.385532e+18 6.514318e+18 7.583486e+18 6.283017e+18 [6] 6.167344e+18 9.820551e+18 7.952821e+18 7.882741e+18 1.744266e+19 [11] 6.653258e+18 8.704722e+18 8.631365e+18 1.87605...

...e 'internal' (or 'working') in function names? What about p.value ? is there some reason for this, as opposed to pvalue ? In the oo languages with which I am most familar, dot is used to denote hierarchy eg text.font.style .. is there a similar intention at work here? John Aitchison -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

Dear R-users, I have an outcome variable and I'm unsure about how to treat it. Any advice? I have spending data for each county in the state of California (N=58). Each county has been allocated money to spend on any one of the following four categories: A, B, C, and D. Each county may spend the money in any way they see fit. This also means that the county need not spend all the money that

A few of you have helpfully sent me references to WMF formats. There is also some code to _read_ it at http://www.csn.ul.ie/~caolan/docs/libwmf.html However, the issue is to _write_ WMF on R under Unix. I am told (by someone who has tried) that the problem is the if you write a WMF file according to the published specs (which were once issued by Microsoft) you find that e.g. Word does not read

Hello all, I'm dealing with geochemical analyses of some rocks. If I use the full composition (31 elements or variables), I can get reasonable separation of my 6 sources. Then when I go onto do LDA with the 6 groups, I get excellent separation. I feel like I should be reducing the variables to thos that are providing the most discrimination between the groups as this is important

Folks, There appears to be a bug(?) in step() when used to screen logistic models. The problem appears to be specific to 1.2.1 (or maybe also 1.2.0?), as the proper behavior was observed in earlier versions. When I did the same on the surrogate log linear model, everything seemed okey. The data involved was the detergent data found in earlier editions of MASS, as given below, > detg1

Dear List, I am looking for a possibility to fit a mixture model under R using maximum likelihood estimation. Venables and Ripley describe a solution working under S+ (in MASS, 3. ed., p. 263) which requires the D system function and deriv3. This solution does not seem to be portable to R or at least I do not realise how. Is there anyone who a) knows how one could make the MASS-method run under

...e: Thu, 13 Jun 2002 08:21:50 +0200 --The supplementary (table of sum one) is fixed for 1.5.1. --Detection code for the first problem has been added to 1.5.1 which will stop --the crash, but the underlying cause is still open. * PR# 1716 * Subject: rdiscrete in e1071 fails when n==1 From: "John Aitchison" <jaitchis@hwy.com.au> Date: Fri, 28 Jun 2002 10:28:17 +1000 * PR# 1729 * Subject: problem with qq( ) From: Jarno.Tuimala@Helsinki.Fi Date: Tue, 2 Jul 2002 10:37:10 +0200 (MET DST) * PR# 1741 * Subject: groupedData constructor from a function From: dieter.menne@menne-biomed.de Date: Thu,...

...e: Thu, 13 Jun 2002 08:21:50 +0200 --The supplementary (table of sum one) is fixed for 1.5.1. --Detection code for the first problem has been added to 1.5.1 which will stop --the crash, but the underlying cause is still open. * PR# 1716 * Subject: rdiscrete in e1071 fails when n==1 From: "John Aitchison" <jaitchis@hwy.com.au> Date: Fri, 28 Jun 2002 10:28:17 +1000 * PR# 1729 * Subject: problem with qq( ) From: Jarno.Tuimala@Helsinki.Fi Date: Tue, 2 Jul 2002 10:37:10 +0200 (MET DST) * PR# 1741 * Subject: groupedData constructor from a function From: dieter.menne@menne-biomed.de Date: Thu,...

I would like to use R to perform a logistic regression with about 30,000 independent variables. That's right, thirty thousand. Most will be irrelevant: the intent is to use the regression to identify the few that actually matter. Among other things, this calls for giving glm() a colossal "y ~ ..." formula with thirty thousand summed terms on its right hand side. I build up the