similar to: logistic regression - glm() - example in Dalgaard's book ISwR

Dear R'ers, When I used S-plus i wrote a small program for a Mantel-Haenszel test for trend (I think it worked). Unfortunately I can't get it working in R. It appears as if my use of 'el' is the problem but I can't sort it out. Error in apply(array, c(, 2, 3), function(el) el * 1:s) : argument is missing, with no default Further down in the program I use 'el'

Dear R People: There is a function from the fBasics library to get the probability and quantiles for the hyperbolic probability function. Is there one that will estimate parms of the hyperbolic probability function from a data set, please? Thanks in advance! Sincerely, Erin Hodgess mailto: hodgess at gator.uhd.edu R Version 2.0.1 windows

library(glmpath) data(heart.data) # heart.data is a list, $y a vector, $x a matrix data <- data.frame(x=I(heart.data$x), y = heart.data$y) > data[1:2,] x.1 x.2 x.3 x.4 x.5 x.6 x.7 x.8 x.9 y 1 160 12 5.73 23.11 1 49 25.3 97.2 52 1 2 144 0.01 4.41 28.61 0 55 28.87 2.06 63 1 > dimnames(heart.data$x)[[2]] [1] "sbp"

Dear all, I am using the publically available GustoW dataset. The exact version I am using is available here: https://drive.google.com/open?id=0B4oZ2TQA0PAoUm85UzBFNjZ0Ulk I would like to produce a nomogram for 5 covariates - AGE, HYP, KILLIP, HRT and ANT. I have successfully fitted a logistic regression model using the "glm" function as shown below. library(rms) gusto <-

Philipp Benner reported a Debian bug report against r-cran-rpart aka rpart. In short, the issue has to do with how rpart evaluates a formula and supporting arguments, in particular 'weights'. A simple contrived example is ----------------------------------------------------------------------------- library(rpart) ## using data from help(rpart), set up simple example myformula <-

I am trying to figure out how to calculate p-values for the difference in prevalence of a risk factor between men and women. For example, I find that 277 out of 710 male patients and 125 out of 305 female patients have obesity, what is the p-value for their difference? If there is a package that can calculate this in bulk, I would appreciate to learn about it! Thank you [[alternative HTML

The mount log file of the volume would help in debugging the actual cause. On Tue, Jun 27, 2017 at 6:33 PM, Joel Diaz <mrjoeldiaz at gmail.com> wrote: > Good morning Gluster users, > > I'm very new to the Gluster file system. My apologies if this is not the > correct way to seek assistance. However, I would appreciate some insight > into understanding the issue I have.

> On Sep 14, 2017, at 12:30 AM, Bonnett, Laura <L.J.Bonnett at liverpool.ac.uk> wrote: > > Dear all, > > I am using the publically available GustoW dataset. The exact version I am using is available here: https://drive.google.com/open?id=0B4oZ2TQA0PAoUm85UzBFNjZ0Ulk > > I would like to produce a nomogram for 5 covariates - AGE, HYP, KILLIP, HRT and ANT. I have

Good morning Gluster users, I'm very new to the Gluster file system. My apologies if this is not the correct way to seek assistance. However, I would appreciate some insight into understanding the issue I have. I have three nodes running two volumes, engine and data. The third node is the arbiter on both volumes. Both volumes were operation fine but one of the volumes, data, no longer

Fixed 'maxiter' in the help file. Thanks. Please give the original source of that dataset. That dataset is a tiny sample of GUSTO-I and not large enough to fit this model very reliably. A nomogram using the full dataset (not publicly available to my knowledge) is already available in http://biostat.mc.vanderbilt.edu/tmp/bbr.pdf Use lrm, not lrm.fit for this. Adding maxit=20 will

This is necessary to handle nested traps to the hypervisor more than one deep. I''ve not seen an actually failure relating to this but I''m not quite sure how we''ve managed to get away with not doing it (I suppose multiply nested traps are uncommon). Signed-off-by: Ian Campbell <ian.campbell@citrix.com> --- xen/arch/arm/entry.S | 4 ++++ 1 files changed, 4

Hello, I have allocated 4 GB memory for a guest instance on the Xen hypervisor. But within the instance there is only 2 GB memory available. See the details, the guest is named pr-idm: # Allocated memory on the hypervisor [root@pr-hyp xen]# xm list Name ID Mem(MiB) VCPUs State Time(s) Domain-0 0 3846 8 r-----

Hello I would like to pass a model formula as an argument to the with.mids function from the mice package. The with.mids functon fits models to multiply imputed data sets. Here's a simple example library(mice) #Create multiple imputations on the nhanes data contained in the mice package. imp <- mice(nahnes) #Fitting a linear model with each imputed data set the regular way works

Hello all, This is the error message that I get. > hyp.res <- nls(log(y)~log(pdf.hyperb(theta,X)), data=dataModel, + start=list(theta=thetaE0), + trace=TRUE) 45.54325 : 0.1000000 1.3862944 -4.5577142 0.0005503 3.728302 : 0.0583857346 0.4757772859 -4.9156128701 0.0005563154 1.584317 : 0.0194149477 0.3444648833 -4.9365149150 0.0004105426 1.569333 :

Hello, I have problems running Xen (from Debian 7, i386 userspace) on Asus P3-AE5 system (M2V-TVM board) with Athlon 64 3000+ CPU (AM2, 1.8GHz). In short, Xen is not accessible from amd64 Linux kernels but works fine with 686 kernels. Long version: xen-system-i386 (4.1.4-3+deb7u1) works fine with linux-image-3.2.0-4-686-pae xen-system-amd64 (4.1.4-3+deb7u1) works fine with

Hi, I was wondering if I might be able to ask some advice about doing residual plots for the lmer function in the lme4 package. Our group's aim is to find if the expression staining of a particular gene in a sample (or "core") is related to the pathology of the core. To do this, we used the lmer function to perform a logistic mixed model below. I apologise in advance

Hi Bill, > ... Still, people would break the builds and let > things go for hours or days at a time. don't forget the time-zone effect. I regularly get build failures in the morning, presumably because someone in the US committed just before going to bed. I guess they are happily snoring away when the build-bots (and humans) start complaining! So when hours go by without a fix, it

Hi Ian, I try compile the bootwrapper for cubieboard2 like this: joshzhao@joshzhao-ThinkCentre-M58p:~/project/Xen/A20/boot-wrapper$ make cubieboard2_defconfig CROSS_COMPILE=arm-linux-gnueabihf- # # configuration written to .config # joshzhao@joshzhao-ThinkCentre-M58p:~/project/Xen/A20/boot-wrapper$ make CROSS_COMPILE=arm-linux-gnueabihf- make -C scripts/kconfig -f Makefile.bootwrapper

I'm plotting obesity rates (y-axis) vs Public Health Unit (x-axis) for the province of Ontario and would like to have the Public Health Unit names appear vertically rather than the default, horizontally. I'm actually using the 'barplot2' function in the {gregmisc} library ... I haven't been able to find a solution in either the barplot2 options or the general plotting

Could someone kindly help me with the following question: when I analyze microarray data I need to fit multiple linear regression models between genes and clinical patameters followed by estimation of the p-values. What's the solution to implement Westfall and Young's algorithm + resampling into the scheme: lm -> stepAIC -> anova. Actually permcor works fine for me in the case of