similar to: confidence intervals around p-values

Dear list, How can I add the mean and standard deviation to each of the boxplots using the example provided in the boxplot function? boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset = supp == "VC", col = "yellow", main = "Guinea Pigs' Tooth Growth", xlab = "Vitamin C dose mg",

Dear R users I?m aware that the package ?modest? is useful to find the mode in an array. However, I?d like to know if someone has translated the ?mode? function built-in in MATLAB into R language. This function finds the most frequent value in an array (http://www.mathworks.com/help/techdoc/ref/mode.html). Best Fer

Dear R people I ask again 1. Is there a published reference presenting the normal score transformation? Is there a published paper (in any field) using that transformation in the analysis of data? And this is a new question 2. The ?outliers? library has a function called ?rm.outlier?. It offers the option of either removing one outlier (the largest or the smallest) or replacing the

Dear people I've learnt that by using the "boxplot.stats" command in the "grDevices" library I can get the 5-number summaries of a boxplot, plus other important information, like the confidence interval around the median. I'm interested in knowing the actual formula to used in that package to calculate that confidence interval. Can someone help me with this? Cheers,

Hi, I'm regressing various body dimensions upon body mass using the 'lmodel2' function, as I'm keen to obtain both OLS and RMA slope values. I also wish to create a plot of the regressions, with the 95% confidence interval of both the slope and intercept. I know how to plot 95% ci bands of the OLS slope using lm with the 'predict' function and 'matlines'. Does

Hello, First of all, when I tried to use function mic I got an error. mic(cbind(C, D)) Error in mic(cbind(C, D)) : could not find function "mic" So I've changed your function myCor and all went well, with a warning relative to BCa intervals. myCor <- function(data, index){ mine(data[index, ])$MIC } results=boot(data = cbind(C,D), statistic = myCor, R = 2000)

Dear R-Experts, Here below is my R code (reproducible example) to calculate the confidence intervals around the spearman coefficient. ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) cor(C,D,method= "spearman") library(boot) myCor=function(data,index){ cor(data[index, ])[1,2] } results=boot(data=cbind(C,D),statistic=myCor, R=2000)

You need: myCor <- function(data, index){ mine(data[index, ])$MIC[1, 2] } results=boot(data = cbind(C,D), statistic = myCor, R = 2000) boot.ci(results,type="all") Look at the differences between: mine(C, D) and mine(cbind(C, D)) The first returns a value, the second returns a symmetric matrix. Just like cor() David L. Carlson Department of Anthropology Texas A&M

Hi We are getting the following error message on Samaba 2.2.0 on AIX 4.3.3:- lib/debug check_log.size (437) File/var/samba/log.smbd failed - using console write_socket_data: write failure. Error = Broken Pipe Can you help please. Robin Smith Technical Director Reflex Data Systems Ltd CONFIDENTIALITY NOTICE This email and any attachments should be read only by

Hi Rui, Many thanks. The R code works BUT the results I get are quite weird I guess ! MIC = 0.2650 Normal 95% CI = (0.9614, 1.0398) The MIC is not inside the confidence intervals ! Is there something wrong in the R code ? Here is the reproducible example : ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) library(minerva) mine(C,D)$MIC

Dear R users: I am using lmer fo fit binomial data with a probit link function: > fer_lmer_PQL<-lmer(fer ~ gae + ctipo + (1|perm) -1, + family = binomial(link="probit"), + method = 'PQL', + data = FERTILIDAD, + msVerbose= True) The output look like this: > fer_lmer_PQL Generalized linear mixed model fit

Hi Anton, Thank you very much for your helpful hint. Thunderbird clearly wants to move Public/Test to the Trash of the user who subscribed the Public folder. Question is how to solve this from a Dovecot point of a view, so a user can also delete folders? 2020200448[10053c700]: proposed url = Public/Test folder for connection INBOX has To Wait = FALSE can run = FALSE 2020200448[10053c700]:

Hi, unfortunately I found a bug in Dovecot's ACL handling for shared mailboxes. It turns out Dovecot doesn't enforce lower casing the privileged username to whom the mailbox should be shared to. This results in a invalid configuration. Users get confused, since they passed on a valid email address in their ACL setup. /usr/local/www/default/mail/test at

Hi, I want to create upper and lower 95% confidence intervals for a p-p plot of an empirical distribution with a theoretical gamma distribution. This is my code: x<-rgamma(100,shape=2, rate=1) # empirical data fitdistr(x,"gamma") # fit a gamma distribution dist<-pgamma(x,shape=1.9884256 ,rate=0.8765314 ) # fitted distribution, using the loglikelihood estimated parameters

Dear All, Apologies if you have a seen a question like this from me before. I am hoping that if I re-word my question more carefully someone may be able to offer more help than the last time I asked something similar. I am using R 2.9.2 and Windows XP. I am trying to determine if there is a statistically significant difference between two c-statistics (or equivalently D statistics). In Stata

Hola Carlos: Te adjunto un ejemplo de aplicación: las funciones (he borrado los path de las funciones y las ordenes "source()" que las carga ) y un ejemplo para ejecutarlas para las opciones que tengo implementadas con la tabla de datos birthwt del paquete"MASS": - Descriptiva de todas las variables de una tabla. - Análisis univariado de todas las variables de una tabla

Hello all. Anybody here had success with install and running RC Reflex XTR simulator in linux using Wine?? I manage to do the install, but failed to install the Reflex usb interface, therefore was unable to control the models. Thanks in advance for any tips, Jorge --------------------------------- Sponsored Link Don't quit your job - take classes online and earn your degree in 1

Estimado Enrique Ramos Yo podría decir ¿y data.table?. Hay muchas alternativas Javier Rubén Marcuzzi De: Enrique RAMOS via R-help-es Enviado: martes, 31 de mayo de 2016 14:03 Para: R-help-es en r-project.org Asunto: Re: [R-es] sumar una variable con cast yo de nuevo, ahora se me presento otro problema en la base de datos del ejemplo solo tenia unas cuantas columnas mi base de datos tiene mas

Hi, I'm the developer of KSambaPlugin. Is there a way to find out if samba was compiled with SSL / LDAP support or not ? I want to disable the options in the case that samba doesn't support them. Thanks, Jan Sch?fer

Hi, Unfortunately it looks like the parameter: "imap_client_workarounds = tb-extra-mailbox-sep" does not work. Neither inside nor outside "protocol imap {}" configuration section. ACLs are set correctly. Manual "DELETE Public/Test" via telnet is working as expected. Is this a known bug? See: - https://help.ubuntu.com/community/Dovecot#Thunderbird -