similar to: (no subject)

Displaying 20 results from an estimated 30000 matches similar to: "(no subject)"

2017 Sep 14
1
Help understanding why glm and lrm.fit runs with my data, but lrm does not
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
2013 May 12
1
Multinomial-Dirichlet using R
Hi: I have asked this question on Cross-Validated. So it might be a cross posting but havent received any responses to it. I am trying to see which distribution will best fit the data I am working on. The dataset is as following: Site Nausea headache Abdominal Distension 1 17 5 10 2 12
2013 Apr 15
2
regression with paired left-censored data
HI I am trying to analyse data which is left-censored (i.e. has values below the detection limit). I have been using the NADA package of R to derive summary statistics and do some regression. I am now trying to carry out regression on paired data where both my X and Y have left-censored data within them. I have tried various commands in R: rega = cenreg(Cen(conc, cens_ind) ~ Gp_ident))? with
2017 Sep 14
0
Help understanding why glm and lrm.fit runs with my data, but lrm does not
> 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
2017 Sep 14
3
Help understanding why glm and lrm.fit runs with my data, but lrm does not
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 <-
2013 Jan 10
2
Titles - main and subtitle won't plot with errbar
Hi, I'm struggling with errbar graphics. I'm trying to plot an x-y graph with correct labelling, however can't seem to get main and sub to show on my graph. They do work when I use title(main="," etc...., but this will make it look at lot messier,I'll have to blank out ylab=" " , and I need to try and get the titles to update automatically according to my
2009 Apr 03
2
Schoenfeld Residuals
Dear All, Sorry to bother you again. I have a model: coxfita=coxph(Surv(rem.Remtime/365,rem.Rcens)~all.sex,data=nearma) and I'm trying to do a plot of Schoenfeld residuals using the code: plot(cox.zph(coxfita)) abline(h=0,lty=3) The error message I get is: Error in plot.window(...) : need finite 'ylim' values In addition: Warning messages: 1: In sqrt(x$var[i, i] * seval) : NaNs
2010 May 10
1
R algorithm/package for creating spatial autocorrelation of uniformly distributed landscape values
Dear all: I would like to create a landscape of environmental values that follow a uniform frequency distribution and also have spatial autocorrelation in the landscape. I was wondering if there is an algorithm and/or package out there that creates autocorrelation of values that are distributed according to a non-normal frequency distribution. Any suggestions are greatly appreciated. Thank you,
2009 Jun 09
2
Isolating a single plot from plots produced simultaneously
Dear R-Help, I am using the 'mfp' package. It produces three plots (as I am using the Cox model) simultaneously which can be viewed together using the following code: fit <- mfp(Surv(rem.Remtime,rem.Rcens)~fp(age)+strata(rpa),family=cox,data=nearma,select=0.05,verbose=TRUE) par(mfrow=c(2,2)) plot(fit) They can be viewed separately but the return key must be pressed after each graph
2010 Dec 17
2
newbie question on str output
Hi! Here is some output from an str command: > str(CO2) Classes ‘nfnGroupedData’, ‘nfGroupedData’, ‘groupedData’ and 'data.frame': 84 obs. of 5 variables: $ Plant : Ord.factor w/ 12 levels "Qn1"<"Qn2"<"Qn3"<..: 1 1 1 1 1 1 1 2 2 2 ... $ Type : Factor w/ 2 levels "Quebec","Mississippi": 1 1 1 1 1 1 1 1 1 1 ... $
2018 May 13
0
(no subject)
> On May 12, 2018, at 9:42 AM, malika yassa via R-help <r-help at r-project.org> wrote: > > > hello > for exampl, i have this programme > # Generating data which are right truncated > library(DTDA) > library(splines) > library(survival) > n<-25 > X<-runif(n,0,1) > V<-runif(n,0.75,1) > for (i in 1:n){ > while (X[i]>V[i]){ >
2012 Oct 08
1
weighted cumulative distribution with ggplot2
Dear all, I am trying to draw a weighted cumulative distribution (as defined here http://rss.acs.unt.edu/Rdoc/library/spatstat/html/ewcdf.html) with ggplot2 however the syntax temp<-qplot(X,weight=weight,data=data,stat = "ecdf", geom = "step",colour=factor(year)) seems not to produce exactly the right figure (the values seems higher at some points)... I am wrong in the
2005 Feb 16
4
Passing colnames to graphics title
Hi, Just a quick query - if I'm creating a function to produce a number of histograms per page of output (one per column from a matrix), how can I pass the column name of the matrix into the title (or indeed to form part of the x-axis label)? TIA, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax:
2017 Dec 25
2
LSD-test
The model should be class aov or lm and my model class is aovlist. tried tidy from broom library but did not work. To make it class aov, I had to remove the error term; model <- aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen),data=Rotationdata_R) Ahmed Attia, Ph.D. Agronomist & Soil Scientist On Mon, Dec 25, 2017 at 7:38 PM, David Winsemius <dwinsemius at
2011 May 05
1
Confidence interval for difference in Harrell's c statistics (or equivalently Somers' D statistics)
Dear All, I am trying to calculate a 95% confidence interval for the difference in two c statistics (or equivalently D statistics). In Stata I gather that this can be done using the lincom command. Is there anything similar in R? As you can see below I have two datasets (that are actually two independent subsets of the same data) and the respective c statistics for the variables in both cases.
2009 Jul 31
1
Displaying function arguments using a Windows R console
I am relatively new to R, and would appreciate any suggestions you may have. I noticed on a Mac the functions' arguments are listed at the bottom of the R console. Is it possible to add such a feature to a Windows R console? I have Windows XP if that helps. I know function arguments can be found using args(...), but I was wanting to have something more automatic, like what I saw on the Mac
2019 Apr 11
1
failed: read(/var/run/dovecot/dns-client)
??????? Original Message ??????? On Thursday, April 11, 2019 9:05 AM, Aki Tuomi <aki.tuomi at open-xchange.com> wrote: > > On 11 April 2019 11:02 Laura Smith via dovecot dovecot at dovecot.org wrote: > > ??????? Original Message ??????? > > On Thursday, April 11, 2019 12:55 AM, John Fawcett via dovecot dovecot at dovecot.org wrote: > > > > > On 11/04/2019
2012 Sep 29
3
Removing lower whisker in boxplot to see the effects of the high values
Good Afternoon- I was wanting to alter the boxplot to remove the lower whisker, both the whisker line and staple just on the lower end. Is there a way to do this? As my code is currently: boxplot(log_loads~ind,data=nfmaum, horizontal=TRUE, notch=T, outline=FALSE, whisker=0, main="Maumee River Near Future Climate Scenarios", ylab="Log Load",xlab="Climate Scenarios")
2005 Mar 01
3
Reconstructing Datasets
Hi, Is it possible to recreate "smoothed" data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages and don't seem to find anything relevant. Thanks in advance, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596
2019 Apr 11
1
failed: read(/var/run/dovecot/dns-client)
??????? Original Message ??????? On Thursday, April 11, 2019 9:01 PM, John Fawcett via dovecot <dovecot at dovecot.org> wrote: > On 11/04/2019 10:02, Laura Smith via dovecot wrote: > > > ??????? Original Message ??????? > > On Thursday, April 11, 2019 12:55 AM, John Fawcett via dovecot dovecot at dovecot.org wrote: > > > > > On 11/04/2019 00:51, Laura Smith