similar to: restricted cubic spline in FGR function

Dear Geskus, I want to develop a prediction model. I followed your paper and analysed thro' weighted coxph approach. I can develop nomogram based on the final model also. But I do not know how to do internal validation of the model and subsequently obtain calibration plot. Is it possible to use Wolbers et al Epid 2009 approach 9 (R code for internal validation and calibration) . It is

Dear All, I would like to use R function 'selectFGR' of fine gray model in competing risks model. I used the 'Melanoma' data in 'riskRegression' package. Some of the variables are factor. I get solution for full model but not in variable selection model. Any advice how to use factor variable in 'selectFGR' function. The following R code is produced below

Hi , library(foreign) write.dta(data1, "data1.dta") should work. The file will be saved in the working directory. Use getwd() to know the working directory. Best wishes Amalraj Raja -----Original Message----- From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of rosario scandurra Sent: 22 March 2018 07:47 To: r-help at r-project.org Subject: [R] exporting data to stata

Hi, I am new to R and I want to export data into Stata. Could somebody help with that? Thanks a lot. This is the code I am using: > setwd("D:/datasets/Seg-bcn/ESBD") > data1 <- readRDS("r17045_ESDB_Habitatges_BDD_V_1_0.rds") > library(foreign) > write.dta(data="data1", file = "D:/datasets/data1.dta") Error in write.dta(data =

On Thu, Mar 22, 2018 at 4:52 AM, Raja, Dr. Edwin Amalraj <amalraj.raja at abdn.ac.uk> wrote: > Hi , > > library(foreign) > write.dta(data1, "data1.dta") > > should work. I don't think so: > library(foreign) > example(svydesign) > write.dta(dstrat, "~/Downloads/foo.dta") Error in write.dta(dstrat, "~/Downloads/foo.dta") : The

Dear Raja, A Fine and Gray model can be fitted using the standard coxph function with weights that correct for right censoring and left truncation. Hence I guess any function that allows to perform stepwise regression with coxph should work. See e.g. my article in Biometrics https://doi.org/10.1111/j.1541-0420.2010.01420.x, or the vignette "Multi-state models and competing risks" in the

Hello, We previously had Win2008 R2 DCs which we have begun to replace with Win2012 servers (forest is still at 2008 level however). I was able to query DFS via rpcclient but now I find that I am unable to. I know that SMB3 was introduced in Win2012 but I believe that it should auto-negotiate earlier versions. I did however test from v4.1.9 as well, unfortunately that still fails. From Samba3

Hi, Sorry not to provide R-code in my previous mail. R code is below #install.packages("rms") require(rms) #install.packages("mstate") library(mstate) require(splines) library(ggplot2) library(survival) library(splines) #install.packages("survsim") require(survsim) set.seed(10) df<-crisk.sim(n=500, foltime=10, dist.ev=rep("lnorm",2),

Dear R users, I am new to R and wanted to apply competing risk methods in my research work. I used the R code given by Zhang et al in his paper 'Nomogram for survival analysis in the presence of competing risks published in Ann Trans Med 2017:5(20):403. I am struggling with getting calibration curve thro' internal validation. I am happy to receive suggestion in the coding as well

Dear R-List, I have a dataframe X containing factor f and numeric variable x1, x2, ... I want to create a new dataframe (or possibly a matrix) that gives statistics (e.g. sum) for the variables x1, x2, ... in each group defined by factor f. What is the simplest way to do this? I tried: attach(X) Z <- data.frame(f=levels(f),x1=as.vector(lapply(split(x1,f),sum))) and stumbled on: Error

I cobbled this together based on a few previous posts and some tinkering, nothing fancy (my first bit of Ruby infact) but it seams to work for my purposes. My questions is actually why can I call <%= h(child.age(child.dob)) %> in a view but not <%= h(child.age_year_month(child.dob)) %>? Both work fine in the console. Thanks for any insight, Hugh child.rb class Child < Person

Hi! The postscript device in 0.90 does not produce symbols (pch=) as it should. E.g.: postscript(file="test.ps",width=8,height=6.4,horizontal=FALSE) plot(1:10,1:10,pch=1:5) legend(1,9,letters[1:5],pch=1:5) dev.off() The result looks quite different compared to the X-windows version. The symbol 1 (pch=1) is much to small, while the other symbols seem to big. Is it a bug? How can I

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

Hi! Since installation of R 0.90 there are problems with the size of text and numbers also in the x11-device (I reported about similar problems with the postscript device recently). For numbers and letters a different size is used: plot(1:10,1:10,type="n") text(2,8,"F5F5F5F5F5") text(2,7,"F5F5F5F5F5",cex=2) The same problem is also found in the postscript

Hi! Since installation of R 0.90 there are problems with the size of text and numbers also in the x11-device (I reported about similar problems with the postscript device recently). For numbers and letters a different size is used: plot(1:10,1:10,type="n") text(2,8,"F5F5F5F5F5") text(2,7,"F5F5F5F5F5",cex=2) The same problem is also found in the postscript

Dear list, I want to perform a biplot, using customized titels for the x and y axis. Setting xlab="" and ylab="" resulted in an error, e.g.: > data(USArrests) > biplot(princomp(USArrests),xlab="",ylab="") Error in biplot.default(t(t(scores[, choices])/lam), t(t(x$loadings[, : length of dimnames[1] not equal to array extent > How do I

Dear R-Users! My aim is to produce boxplots without the outliers included in the plot. I started to write a function that looks something like: myboxplot <-function(x,fa) { bpdata <- boxplot(x~fa,plot=FALSE) bpnames <- names(bpdata) for (JJ in bpnames) { command <- paste("bpdata$",JJ,"$out <- numeric(0)",sep=""); eval(command) }

From: "Wolfgang Koller" <koller2@fgr.wu-wien.ac.at> Organization: University of Economics, Vienna Date: Tue, 23 Nov 1999 16:29:22 CET (+0100) MIME-Version: 1.0 Content-type: text/plain; charset=US-ASCII Content-transfer-encoding: 7BIT CC: r-bugs@biostat.ku.dk, r-help@stat.math.ethz.ch References: <153DF9230E8@fgr.wu-wien.ac.at> X-mailer: Pegasus Mail

Hi, I've searched in the wiki and in the mailinglist archives but haven't found anything about password character sets within the dovecot authentication deamon. My problem is that we have users with passwords containing scandinavian characters (???, umlauts) and the debug log shows that different clients send the password in different charsets. The passwords are stored in a Mysql

Dear R-developpers, currently the function window shows the following behaviour, if start or end is out of range: > x <- ts(1:14,start=c(1950,2),frequency=4) > y <- window(x,start=c(1950,1),end=c(1953,4)) Warning messages: 1: start value not changed in: window.default(x, ...) 2: end value not changed in: window.default(x, ...) Instead of this I would like it to fill with NAs the