similar to: simultaneous confidence intervals for multinomial proportions: sample size

Dear all, Does someone know an R function implementing the method of Sison and Glaz (1995) (see full ref below) for computing Simultaneous confidence intervals for multinomial proportions? As alternative method, I think to boostrap the mean of each proportion and get in that way confidence interval of the mean. I observed 21 times a response that could be one out of 8 categories

While using the Hodges Lehmann Mean in DescTools (DescTools::HodgesLehmann), I found that it generated incorrect answers (see <https://github.com/AndriSignorell/DescTools/issues/97> https://github.com/AndriSignorell/DescTools/issues/97). The error is driven by the existence of tied values forcing wilcox.test in Base R to switch to an approximate algorithm that returns incorrect results - see

Hi R users! Is there a function that extracts the simultaneous confidence intervals for multinomial proportions as described by Sison and Glaz 1995? or anyone else for that matter? I have seen that SAS has macro for this ( http://www.jstatsoft.org/v05/i06) and i was wondering if R had something similar? Thank you and have a great day ahead! ---------------------------------

Greetings, I am having a problem with DescTools::Quantile (a function computing quantiles from weighted samples): # these sum to one probWeights = c( 0.0043, 0.0062, 0.0087, 0.0119, 0.0157, 0.0204, 0.0257, 0.0315, 0.0378, 0.0441, 0.0501, 0.0556, 0.06, 0.0632, 0.0648, 0.0648, 0.0632, 0.06, 0.0556, 0.0501, 0.0441, 0.0378, 0.0315, 0.0257, 0.0204, 0.0157, 0.0119, 0.0087,

Hello, I have three numerical variables and I would like to test if their correlation is significantly different. I have seen that there is a package that "Test the difference between two (paired or unpaired) correlations". [https://www.personality-project.org/r/html/paired.r.html] However, there is the need to convert the correlations to "z scores using the Fisher r-z

It looks like a homework assignment. It also looks like you didn't read the documentation carefully enough. The 'len.out' argument in seq is solely for specifying the length of a sequence. The 'quantile' function omputes the empirical quantile of raw data in the vector 'x' at cumulative probabilit(y)(ies) given in the weights' argument, with interpolation I'm

Hi Everyone, Following the previous discussion on optimizing *wilcox functions, Andreas Loeffler brought to my attention a few other bugs in `wilcox` family functions. It seems like these issues have been discussed online in the past few months, but I haven?t seen discussion on R-devel...unless I missed an email, it seems like discussion never made it to the mailing list. I haven?t seen any bug

Hello, I did some nice grafics using plot 3d and scatter3d. Is there a, hopefully not too complicated way, to transform these dynamic, three-dimensional and interactive grafics in a kind of executable file? I want to show and send them via e-mail to projekt partners who don`t use GNU R and who are not used to do programming. So the result should be quite comfortable to execute for them in

¿Alguien sabe si el test de reacciones extremas de Moses está escrito en algún paquete de R? Gracias de antemano.

Thank you, but this now sounds more difficult: what would be the point in having these ready-made functions if I have to do it manually? Anyway, How would I implement the last part? On Thu, Mar 23, 2023 at 1:23?AM Ebert,Timothy Aaron <tebert at ufl.edu> wrote: > > If you are open to other options: > The null hypothesis is that there is no difference. > If I have two equations

Hello, I took my first steps in dynamic reports with Gnu R and used sweave(). I therefore run Sweave() with an example of Friedrich Leisch, starting like this: \ documentclass [ a4paper ]{ article } \ title { Sweave Example 1} \ author { Friedrich Leisch } \ begin { document } \ maketitle and so on. It worked very well but when I tried to compile the latex file I get this error message: !

Hi!   I want to print the output all together with a single column name   s21<-c(1:1000); var21<-lapply(s21,function(x){    ns<-rnorm(78,8,9);    n<-length(ns);    Mn<-mean(ns)    Sn2<-var(ns)    return(cbind(x,Mn,Sn2)); }); var21 but my code is giving me somewhat like the following [[1]] x   Mn          Sn2   [1,] 1 7.86 10.56540 [[2]] x   Mn          Sn2  

Hi all, I have the following code in my application layout: <div id="utility"> <% if logged_in? %> <%= link_to "Sign Out", logout_path %> <% else %> <%= link_to_function "Sign in", "$(''signin'').toggle();$(''signin_link'').toggle()", :id => "signin_link" %>

I am doing a cluster analysis with hclust. I want to get hclust to output the Hotelling's T squared statistic for each cluster so I can evaluate is data points should be in a cluster or not. My research to answer this question has been unsuccessful. Does anyone know how to get hclust to output the Hotelling's T squared statistic for each cluster? Mike [[alternative HTML version

I have a two part question one about statistical theory and the other about implementations in R. Thank you for all help in advance. (1) Am I correct in understanding that Heteroscedasticity is a problem for Generalized Additive Models as it is for standard linear models? I am asking particularly about the GAMs as implemented in the mgcv package. Based upon my online search it seems that some

I present you with a function that solves a problem that has bugged me for many years. I think the problem may be general enough to at least consider adding this function, or a revamped version of it, to the 'stats' package, with the other multinomial functions reside. I'm using R to export data to text files, which are input data for an external model written in C++. Parts of the

Definitely yes I just tried a one line latex document and it bombed with \begin{ document } but is fine with \begin{document} . I think the spaces within the {} must all be removed. Latex does not like spaces John Kane Kingston ON Canada > -----Original Message----- > From: jdnewmil at dcn.davis.ca.us > Sent: Thu, 07 Apr 2016 05:42:57 -0700 > To: paladini at trustindata.de, r-help

Lets say I have an * at my business, with 7960 SIP phones. All the sip phones are registered using their extension number (like 305), but I would also like to put my SIP URI on my business card and in a name format, not an extension number (like lee.goodman), so that the SIP URI would read lee.goodman@asterisk.company.com. How would I set this up in extensions.conf? I got

the method is described in the article:goodman leo A.,a modified multiple regression approch to analysis of dischotomous variables",american sociological review 33(hebruary):28-46 thank you in advance:)

The Goodman/Kruskal gamma is a nice descriptive rank-order correlation statistic, often used in psychology. It is nice because it is easy to understand. It takes all pairs of values of each variable and asks whether they are congruent (S+ is the number in the same order for both variables) or discordant (S-, opposite ranking). The statistic is (S+ - S-)/(S+ + S-). It is like tau except for the