search for: meds

Hi, If I set the p.adjust="none", does it meant that the output p values from the pairwise.t.test will be the same as those from individual t.tests (set var.equal=T, alternative="t")? I actually got different p values from the two tests. See below. Is it supposed to be this way? Thanks Johnny > x [1] 61.6 52.7 61.3 65.2 62.8 63.7 64.8 58.7 44.9 57.0 64.3 55.1 50.0 41.0

Hello, I am using Dr. Harrell's design package to make a nomogram. I was able to make a beautiful one without stratifying, however, I will need to stratify to meet PH assumptions. This is where I go wrong, but I'm not sure where. Non-Stratified Nomogram:

Hi, I have newbie question. Suppose I have the following data: temp <- data.frame(type1 = c("male", "female", "male", "female", "female"), type2 = c("low", "med", "high", "low", "med"), a = c(1,2,4, NA, 3), b = .... [TRUNCATED] temp type1 type2 a b c 1 male low 1 5 0 2 female

Hello I am using sem to look at the direct effect of one variable on another but i am uncertain if i am progressing correctly. An example: covar1<-? matrix(c(0.4,-0.2,3,-0.2 , 0.3,-2 , 3 ,-2 , 60), nrow=3,byrow=T) rownames(covar1)<-colnames(covar1)<-c("endo","exo","med") path1<-matrix(c(? ? "exo -> endo",? "g1", NA,

The following two patches allow the parallel use of memsharing, xenpaging and xen-access by using an independent ring buffer for each feature. Please review. v2: - update mem_event_check_ring arguments, check domain rather than domain_id - check ring_full first because its value was just evaluated - check if ring buffer is initialized before calling mem_access_domctl/mem_paging_domctl

Estimados, tengo un grafico al que agregue una linea con abline(), este es el codigo ... library(Hmisc) # agregar una linea horizontal con el valor del angulo de contacto del control trellis.device(color = FALSE) xYplot(Cbind(ca.med[,3],ca.med[,5],ca.med[,6])~ca.med[,2], groups=Sol, data=ca.med, xlab=list(label="Concentration (mM)", cex=1.3), ylab=list(label="Contact Angle

Hello. I'm a newbie here. In my script (I name it readData.R), I wrote the followings: readData <-function(){ med = read.csv("medicalData.csv");} Then I tested the script by 'Source R Code' then on the command I typed 'readData()' then I typed 'med' to check if the variable contains the medical data but it returned 'Object med is not found'. What

Dear colleagues, What did I not understand ? ->my intention I want to create a new variable: In plain language: If someone is taking anithypertensive treatment (med.hyper==1) table(med.hyper) med.hyper 0 1 472 97 I want to subtract 5 mmHg (rr.dia.2m-5) from the measured diastolic blood pressure (rr.dia.2m) if not treated - the value of the measured diastolic blood pressure should

Hello, I'd like to remove the individual pairs of square brackets along with their content - plus the space directly behind it - from address strings such as this: [Swidsinski, Alexander; Loening-Baucke, Vera; Lochs, Herbert] Charite Humboldt Univ, Innere Klin, D-10098 Berlin, Germany; [Hale, Laura P.] Duke Univ, Med Ctr, Dept Pathol, Durham, NC 27710 USA I'd like get the

Dear all, This relatively routine analysis has left me frustrated and in a rut. I have a dataset (data1), which I subset in order to remove rows where HabitatDensity="Med". This dataset looks correct when I call it up, however, when I create a table out of the new subset (data2), my table continues to show the "Med" information as 0. This is a problem because I need a

Dear all, I'm trying to analyze data with the following structure: ind cat tx age 40.2 por fol peq vh 35 41.9 por fol med vh 35 68.9 por fol preov vh 35 71.5 por fol peq ser 37 67.5 por fol med ser 37 76.9 por fol preov ser 37 78.7 por fol peq otr 37 78.3 por fol med otr 37 82.1 por fol

Good afternoon, I would like to ask for help in controlling y-axis limits and labels in lattice doplots. Unfortunately, the problem is somewhat convoluted, please bear with the long explanation. I would like to create a 3x3 lattice of dotplots, say subject ~ count. The plot is conditioned on variables treatment and risk: subject ~ count | treatment + risk. In the experiment, not all subjects

I have a data frame that I plan on importing from a file (there are other columns of numeric data): > xy x y 1 1 1 2 1 1 3 1 1 4 1 1 5 2 2 6 2 2 7 2 2 8 2 2 > where x is a column of factors, and y is a column of factors, that have different levels, e.g., > x [1] low low low low hi hi hi hi Levels: low hi > y [1] med med med med medhi medhi medhi medhi Levels: med

Hello everyone, First i explain my aim. I would like to calculate the frecuency of temperature (datos.mx1, columns 3-6) for each month with the thresholds of table lf.med and lc.med. numero1 <- hist(subset(datos.mx1[,3], datos.mx1$Mes==1),plot=FALSE,breaks=c(min(datos.mx1[,3]),lf.med[1,2],lc.med[1,2],max(datos.mx1[,3]))) # First month and only the first weather station (column 3) I now

Confusing. Github says that is open. Ok. My mistake. The question remains why the Windows clients allow login for an expired certificate despite a correctly loaded CRL. What is the purpose of specifying the CRL in smb.conf? It seems to me that the smartcard login is not really reliable. Then my users still have to log in with password. For now, as long as 4.19 is not yet released. Hans

Hi! Does anyone know how to do the test for goodness of fit of a logistic model (in rms package) after running fit.mult.impute? I am using the rms and Hmisc packages to do a multiple imputation followed by a logistic regression model using lrm. Everything works fine until I try to run the test for goodness of fit: residuals(type=c("gof")) One needs to specify y=T and x=T in the fit. But

Dear R-listers, I have a little problem with a boxplot and I hope you can help me figuring it out. I'll try to make up some data to illustrate the issue. Sorry, if my procedures look naive, but these are my first steps in R. Any comments and/or suggestions are very welcome. let's create a vector var1: var1 <- rnorm(100) and 5 five logical vectors. In this case the vectors don't

Hi there, I recently posted about problems with winbind resolving names from ADS. We have an ADS Domain Controller serving the domain MED-DUS and an NT 4 PDC serving the domain MEDOIL. We want to use a samba server as memberserver in the MED-DUS domain. I joined the samba server to the MED-DUS domain. The domain join worked ok. All users (including MEDOIL) are visible. Next i started winbind.

Hello, Any advice or pointers for implementing Sobel's test for mediation in 2-level model setting? For fitting the hierarchical models, I am using "lme4" but could also revert to "nlme" since it is a relatively simple varying intercept model and they yield identical estimates. I apologize for this is an R question with an embedded statistical question. I noticed that a

Muchas gracias Freddy y Carlos ... estuve intentando con .() y con list(), para calcular la media y el error estandar al mismo tiempo en dos columnas, pero me arrojaba un error que no supe interpretar. Ahora ya funciona como sugiere Carlos. Muchas gracias de nuevo. Saludos, Eric. On 14/05/15 19:28, Carlos Ortega wrote: > Hola, > > La forma de hacerlo con data.table es esta: >