Re: [R] Error message when using 'optim' for numerical maximum likelihood

Good morning, How can I resolved error message New_X with convergence 1.Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] restricted cubic spline in FGR function

Dear Thomas, I want to use evaluate effect of Age using restricted cubic form in the FGR function as Fgr.crr <- FGR(Hist(time, event) ~ rcs(Age_years), data=dat) It provides error. " Error in parse(text = termtext, keep.source = FALSE): 1: response ~ rcs(Age_years Do I need to

Re: [R] exporting data to stata

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-boun...@r-project.org] On Behalf Of rosario scandurra

[R] selectFGR vs weighted coxph for internal validation and calibration curve- competing risks model

? Thank you very much for all your inputs and suggestions Regards Amalraj raja -Original Message- From: Ronald Geskus [mailto:statist...@inter.nl.net] Sent: 21 March 2018 04:01 To: r-help@r-project.org Cc: Raja, Dr. Edwin Amalraj <amalraj.r...@abdn.ac.uk> Subject: Re: [R] sel

[R] selectFGR - variable selection in fine gray model for competing risks

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

Re: [R] Competing risks - calibration curve

L, estimates=TRUE, bw=FALSE, aics=0, what="observed-predicted", tol=1e-12, maxdim=5) plot(cal, subtitles=FALSE) plot(cal) cal.km <- calibrate(mod, u=1, cmethod='KM', m=10, B=10, pr=FALSE) plot(cal.km) cal.hare = calibrate(mod, u=1, cmethod='hare', m=20, B=20) plot(cal.h

[R] Competing risks - calibration curve

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

Re: [R] Differences between SPSS and R on probit analysis

is akin to the SPSS method, and thus produces nearly-identical results. You may still see very, very minor differences in chi square goodness of fit, and 95% CI of the doses/concentrations, etc. but this is due to differences in rounding under the hood of the software. Hope this helps! Edwin R

[R] Embed fonts in an R graph

appreciated. Edwin http://r.789695.n4.nabble.com/file/n4709707/Capture.png -- View this message in context: http://r.789695.n4.nabble.com/Embed-fonts-in-an-R-graph-tp4709707.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r

[R] Rcpp cpp11 and R CMD build

... OK * cleaning src * installing the package to process help pages --- * installing *source* package ‘fluEvidenceSynthesis’ ... ** libs g++ -I/usr/share/R/include -DNDEBUG -I/home/edwin/R/x86_64-pc-linux-gnu-library/3.2/Rcpp/include -I/home/edwin/R/x86_64-pc-linux

Re: [R] Rcpp cpp11 and R CMD build

Thank you! I was missing the SystemRequirements. I guess it could be useful to add this to the example given here: http://gallery.rcpp.org/articles/simple-lambda-func-c++11/ Cheers, Edwin On Wed, 24 Jun 2015 at 17:50 Charles Determan cdeterma...@gmail.com wrote: Hi Edwin, If you look

[R] sprintf doesn't care of escape characters

Hi. I am having trouble with something that should be simple. I am unable to get sprintf using escape sequences: sprintf(a\nb) [1] a\nb sprintf(a\bc\d) [1] a\bc\d But it seems to need them any way: sprintf(a\bcd) Error: unexpected symbol in sprintf(a\bcd Any suggestion on how to solve

Re: [R] sprintf doesn't care of escape characters

characters, but if you use cat, or output to a file, you will see that the result is correct: cat( sprintf(a\nb)) a b cat(sprintf(a\bc\d)) abcd On Mon, Jul 30, 2012 at 1:29 PM, Edwin Helbert Aponte Angarita helber...@gmail.com wrote: Hi. I am having trouble

[R] read.xls question

Hello, I've been using: tmp.df = read.xls(filename, stringsAsFactors = FALSE) to read in my files. Even though I get the There were 50 or more warnings thing, for the most part most of the data is read in correctly. However, there are a few select rows where there are values but they are

[R] undefined S4 class in parallel computing at snowfall

Dear All, I have a question of how to export S4 class specification to clusters/workers in parallel computing. The package I used is snowfall. The problem is reproducible as follows. Any hint is greatly appreciated. Edwin Sun === begin === library(snowfall) sfInit(parallel = TRUE, cpus

Re: [R] I cannot get species scores to plot with site scores in MDS when I use a distance matrix as input. Problems with NA's?

copping my TinnR script. It contains not only the R script but also the errors returned in the R console. I apologize that the script is still not really reproducible but I just do not know how to create a sample data frame for it. Thanks again for your help. Edwin. ##Load required packages

Re: [R] I cannot get species scores to plot with site scores in MDS when I use a distance matrix as input. Problems with NA's?

, it seems to be able to handle NAs and non-dummy coded character variables metaMDS(daisy(df, metric=gower)) Edwin Lebrija Trejos wrote Hi, First I should note I am relatively new to R so I would appreciate answers that take this into account. I am trying to perform an MDS ordination

[R] I cannot get species scores to plot with site scores in MDS when I use a distance matrix as input. Problems with NA's?

= type, display = display, :Species scores not available“ Can anyone guide me to get the trait (“species”) scores to plot together with my species (“site”) scores? Thanks in advance, Edwin __ R-help@r-project.org

[R] Help with R plots

Hi all, Mind my ignorance and complete newbiness, but I am very new to R and am currently learning it through my stats course at university. My question is regarding adding titles to plots. I'll post the code below, and then explain my issue: hist(gender$HeartRate) title(Histogram: HeartRate

[R] strange problem with strptime and date variable

Hello all, I am manipulating a large database with 70,000 records. strptime generates a date variable but R treats some of the values as NA. I attach a simple example below. I have spent hours on this problem. Any hint would be greatly appreciated. Many thanks, Edwin Sun # ===start

Re: [R] strange problem with strptime and date variable

David - Thanks very much. That is a great hint. I believe that is the solution. Edwin Sun -- View this message in context: http://r.789695.n4.nabble.com/strange-problem-with-strptime-and-date-variable-tp3677178p3677205.html Sent from the R help mailing list archive at Nabble.com

[R] Piecewise distribution function estimation with Generalized Pareto for tail

, Edwin Sun -- View this message in context: http://r.789695.n4.nabble.com/Piecewise-distribution-function-estimation-with-Generalized-Pareto-for-tail-tp3649961p3649961.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] How to Plot Two Curves Into One Page

(vectorOfxValues, FooBarAsMatrix, lwd=2, lty=3, col=c(red, blue)) legend(0.6,0.6,c('Default','Probabilistic'), col=c('red','blue'),lwd=3) Regards, Edwin -- Dr. Edwin Groot, postdoctoral associate AG Laux Institut fuer Biologie III Schaenzlestr. 1 79104 Freiburg, Deutschland +49 761-2032945

Re: [R] how to control ticks

to control ticks. Dear Yogesh, I spray my clothing with No-Bite, and that controls ticks quite well. :-) Edwin plot(file$time, file$ch4*1000, ylim=c(1500,1700), xaxt='n', xlab= NA, ylab=NA,col=blue,yaxs=i,lwd=2, pch=10, type=b)# axis(1, at=1:12, labels=c('J','F','M','A','M','J','J','A','S','O

Re: [R] Help Please!!!!!!!!!

file, I would delete the white space before None, save the file, and use the following function to read your data into a data.frame: read.delim(smokedata.txt) This assumes you used a tab character between each field. HTH, Edwin -- Dr. Edwin Groot, postdoctoral associate AG Laux Institut fuer

[R] how to update my own function

Hello all, I wrote a function with many arguments. Then I need to call it many times with changes on some arguments only. Is there any way to write a function or have a method to update it, like the relationship between lm() and update()? Many thanks, Edwin Sun This is the sample

Re: [R] how to update my own function

Thank you all the reply. The use of sys.call() as suggested by Ducan works pretty well. This is good as the outputs from my function are large. The use of apply family functions as suggested by Erik is good for small amounts of outputs. Edwin Sun -- View this message in context: http://r

Re: [R] dnorm and qnorm

Thank you all for the great help. I think the optimize function and approach solves my problem well. Edwin Sun -- View this message in context: http://r.789695.n4.nabble.com/dnorm-and-qnorm-tp3040427p3041962.html Sent from the R help mailing list archive at Nabble.com

[R] dnorm and qnorm

(pdf.xb) [1] 0.628605 Many thanks, Edwin -- View this message in context: http://r.789695.n4.nabble.com/dnorm-and-qnorm-tp3040427p3040427.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https

Re: [R] store matrix in an arrary

call, that is the source of your error message. It means you tried to store something in s[] that goes beyond its dimensions. Edwin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Tal Galili wants to stay in touch on LinkedIn

Tal: Invite everyone on R-Help, why not? Please check your contact e-mail address! Edwin -- On Mon, 20 Sep 2010 02:14:17 -0700 (PDT) Tal Galili via LinkedIn mem...@linkedin.com wrote: LinkedIn Tal Galili requested to add you as a connection on LinkedIn

Re: [R] Multiple CPU HowTo in Linux?

of the packages I have used so far. The HPC package guides lacked something in the practical aspects of their usage. lapply - mclapply at the beginning of my script? Well, I never would have thought of such a thing. Thanks! I might be back on the list when I run out of physical RAM ;-) Edwin

[R] Multiple CPU HowTo in Linux?

--- Thanks for your help, Edwin -- Dr. Edwin Groot, postdoctoral associate AG Laux Institut fuer Biologie III Schaenzlestr. 1 79104 Freiburg, Deutschland +49 761-2032945 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Multiple CPU HowTo in Linux?

, multiprocessing and cluster processing. Edwin On Tue, 14 Sep 2010 10:15:42 -0400 Johnson, Cedrick W. cedr...@cedrickjohnson.com wrote: ?multicore perhaps On 09/14/2010 10:01 AM, Edwin Groot wrote: Hello all, I upgraded my R workstation, and to my dismay, only one core appears to be used

[R] memory problem for scatterplot using ggplot

for additional information: I'm running my R script on my 32-bit machine: Ubuntu 9.10, hardware: AMD Athlon Dual Core Processor 5200B, memory: 1.7GB. Many thanks in advance. Kind Regards, -- Ir. Edwin H. Sutanudjaja Dept. of Physical Geography, Faculty of Geosciences, Utrecht University

Re: [R] ggplot2: How to change font of labels in geom_text

I have the same problem and I wonder if there is any answer from the community. Thanks. -- View this message in context: http://r.789695.n4.nabble.com/ggplot2-How-to-change-font-of-labels-in-geom-text-tp991579p2286671.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] GARCH estimation with exogenous variables in the mean equation

function that can handle exogenous variables in estimating GARCH? Thank you a lot. Edwin -- View this message in context: http://n4.nabble.com/another-GARCH-problem-tp859307p1819640.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r

[R] navigation panel with help

notice there is an index link at the bottom of the page. By the way, I did not make any change on my browser. Is this a change for this version? Thank you for your help. Edwin Sun -- View this message in context: http://n4.nabble.com/navigation-panel-with-help-tp1395663p1395663.html Sent from the R

[R] ggplot legend for multiple time series

- data.frame(year=1875:1972, level=LakeHuron) ggplot(huron, aes(year)) + geom_line(aes(y=level+5, colour=above)) + geom_line(aes(y=level-5, colour=below)) + scale_colour_manual(Direction, c(above=black, below=black)) Thanks, Edwin Changyou Sun, Ph.D. Associate Professor Natural Resource Policy

Re: [R] ggplot legend for multiple time series

Hello - Thank you so much for the help. It works perfectly. I guess that as many have pointed out, ggplot is a great package but there is a lack of documentation and examples. Edwin Sun baptiste auguie-5 wrote: Hi, I don't understand why you used scale_manual_colour if you want only

Re: [R] IP-Address

- o[sort.list(ip[kk,o], method=radix, na.last=TRUE)] } df[o, ] } ) print(res) test elapsed 1 peda 4.12 2 peda2 4.08 3 hb 0.28 4 hb2 0.25 On Sun, May 31, 2009 at 12:42 AM, Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: edwin

Re: [R] IP-Address

waclaw.marcin.kusnierc...@idi.ntnu.no wrote: edwin Sendjaja wrote: Hi VQ, Thank you. It works like charm. But I think Peter's code is faster. What is the difference? i think peter's code is more r-elegant, though less generic. here's a quick test, with not so surprising

Re: [R] IP-Address

wow =P Thank you guys. I will use the hb2 , because it is the fastest one. Eddie wow! :) vQ Henrik Bengtsson wrote: library(gsubfn) library(gtools) library(rbenchmark) n - 1 df - data.frame( a = rnorm(n), b = rnorm(n), c = rnorm(n), ip = replicate(n,

Re: [R] IP-Address

Maybe you can make a package for something like this. (IP sort,etc). Like in Perl, there is a module called Sort::Key::IPv4 I think, not only me has the same problem for sorting e.g IP? Eddie Not really, just the old saying that any piece of code can be made twice as fast (which often holds

Re: [R] IP-Address

Hi Allan, Thank you for your reply. I have more than 4 Rows. I don't believe than it's a good solution, if I need to write down all of the values. I have about 200 Rows. Do you have any other solution?. Is there any package that I can use? Edwin Date: Fri, 29 May 2009 09:59:40 +0100

Re: [R] IP-Address

Hi Peter, This works as I expected. Thank you very much. Eddie Date: Fri, 29 May 2009 11:15:37 +0200 From: p.dalga...@biostat.ku.dk To: all...@cybaea.com CC: r-help@r-project.org Subject: Re: [R] IP-Address Allan Engelhardt wrote: IP addresses are very (very!) difficult to

Re: [R] IP-Address

Hi VQ, Thank you. It works like charm. But I think Peter's code is faster. What is the difference? Eddie Date: Fri, 29 May 2009 11:44:15 +0200 From: waclaw.marcin.kusnierc...@idi.ntnu.no To: p.dalga...@biostat.ku.dk CC: r-help@r-project.org Subject: Re: [R] IP-Address Peter

[R] IP-Address

Hi, Is there any way to sort a tabel with a colum with IP-address? table: id rank color status ip 138 29746 yellow no 162.131.58.26 138 29746 red yes 162.131.58.16 138 29746 blue yes 162.131.58.10 138 29746 red no 162.131.58.17 138 29746 yellow no 162.131.58.14 138 29746 red no

Re: [R] help with ggplot2 -- ggpoint function missing?

Thanks, that helped! I didn't realize there was another version of the book. On Wed, Apr 1, 2009 at 2:51 PM, Tobias Verbeke tobias.verb...@openanalytics.be wrote: Hi, I'm trying to follow the ggplot introduction here: http://had.co.nz/ggplot/ggplot-introduction.pdf I've installed ggplot2

[R] Regression

Hello, I have problem I have a, and b in regression then I can't plot x,y with (and) a, b lines can you help me ? thx [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] My Problem

Hello, My name is Edwin, I come from INDONESIA I have problem I creating function then I have many calculation like this xx-function(){ a-sd() b-beta1.hat c-beta2.hat data.entry(a,b,c) } then i have function too, almost same yy-function(){ d-sd() e-beta1.hat f-beta2.hat data.entry(d,e

[R] Project Robust Linier Regresssion

Hello, I'm EDWIN, I create (make) GUI, with call many function but I don't know why when I call function I can't. if without function, Yes I can.. can you help me ? can you make this, become true with full code? Can you help me to create data.entry with interface LM - BETA1.HAT - BETA2.HAT

[R] Package

Hello, My name is edwin, I come from Indonesia, Can you help me, I want package which is contain lqs and lmsreg. I need fast... Can you help me to create bootstrap regression ?? Thx Edwin [[alternative HTML version deleted]] __ R-help@r

[R] Large Dataset

? Thank you in advance, Edwin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Large Dataset

Hi Simon, Thank for your reply. I have read ?Memory but I dont understand how to use. I am not sure if that can solve my problem. Can you tell me more detail? Thanks, Edwin type ?memory into R and that will explain what to do... S - Original Message - From: Edwin Sendjaja

Re: [R] Large Dataset

I think he meant: ?Memory edwin When I do it on a Mac installation I get: Help for the topic memory was not found. Is that a Linux-specific function? Or perhaps you meant to type: ?Memory Which does produce useful information. -- David Winsemius sessionInfo() R version 2.8.0

Re: [R] Large Dataset

Hi Ben, Using colClasses doensnt improve the performace much. With the data, I will calculate the mean, min, max, and standard deviance. I have also failed to import the data in a Mysql Database. I dont have much knowledge in Mysql. Edwin Edwin Sendjaja edwin7 at web.de writes: Hi

Re: [R] Large Dataset

Bellow, you can see the R data. But this stucks even in first line (read.table..). I dont know how to calculate this and write the result into a new table. Edwin data - read.table(test.data) data - subset(data, (data$Zusatz!=60) (data$Zusatz!=0)) list

[R] How to calculate words in column?

Hi, I have a table for an 1 week exam result for many classes in school, like this: Day Class_ID TestResult 1 Monday1 Paper Passed 2 Tuesday1 Oral Passed 3 Friday 1 Paper Passed 4 Monday3

Re: [R] Creating Tabel for Mean, Deviance

Hi, Thank you Stefan, I think, this package will do what I want. Edwin Edwin Sendjaja schrieb: Hi, Can someone help me how to create a table for Min, Max, Mean, Deviance. Try the summaryBy function of the doBy package. The package has an excellent vignette (=documentation). hth

[R] Creating Tabel for Mean, Deviance

Hi, Can someone help me how to create a table for Min, Max, Mean, Deviance. I have input data like this: Person Age Child 1 517 2 2 517 0 3 517 1 4 5202 5 520 3 6 710 1 7 721

Re: [R] trying to superimpose a line plot onto a histogram

- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: July 5, 2008 8:43 PM To: Edwin Lei Subject: Re: [R] trying to superimpose a line plot onto a histogram hi: you can do the same thing as below but just replace density(x) with whatever you want to plot ? On Sat, Jul 5, 2008

Re: [R] trying to superimpose a line plot onto a histogram

Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: July 5, 2008 10:22 PM To: Edwin Lei Subject: RE: [R] trying to superimpose a line plot onto a histogram i'm no R expert but if you send EXACTLY what you want to do to the list, i bet someone can help. if not tomorrow, then i

[R] trying to superimpose a line plot onto a histogram

correspond to the break points in the histogram. Next I used the command points(x,y,type=l) But for some reason, the line plot is shifted to the right and doesn't line up with the histogram. Thanks for the help! Edwin Lei data.pdf Description: Adobe PDF document

Re: [R] trying to superimpose a line plot onto a histogram

Thanks for the quick response. But the points I'm trying to superimpose isn't exactly the density of the data. Is there any other way to do it? From: milton ruser [mailto:[EMAIL PROTECTED] Sent: July 5, 2008 6:56 AM To: Edwin Lei Cc: r-help@r-project.org Subject: Re: [R] trying to superimpose

[R] 3-D-Plot

horisontally. I might wrong with this. Thank a lot in advance Edwin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented

Re: [R] 3-D-Plot

Hi Uwe, Thanks for your answer. What is the different between rgl and scatterplot3d? I dont need a graphik like vulcano. I just need 3D-dot-plot. Am Donnerstag, 17. April 2008 09:03:17 schrieb Uwe Ligges: Edwin Sendjaja wrote: Hello, I've got some problems. I hope someone can help me

Re: [R] 3-D-Plot

, heise.de,etc If no, can you tell me little how to modify the code to get this. Thank you. another question: How can I change the width line of the axis (x,y,z). The default line is for me too thin. Thank you in advance, Edwin Sendjaja Am Donnerstag, 17. April 2008 14:43:39 schrieb Uwe Ligges

Re: [R] Mapping file-Legend- from 2 Files/Tables

, table2, by.x=User_ID, by.y=User_ID) result-data$Name (for example) Thank you. Kind Regards, Edwin Sendjaja Am Samstag, 12. April 2008 19:40:03 schrieb Edwin Sendjaja: Does anyone know how to solve this. I am so desperate. I'd be terribly grateful for any help. Am Samstag, 12

Re: [R] Mapping file-Legend- from 2 Files/Tables

Does anyone know how to solve this. I am so desperate. I'd be terribly grateful for any help. Am Samstag, 12. April 2008 05:19:16 schrieb Edwin Sendjaja: Hello, I have got 2 Tables from different files: Table 1 (lets say file: Salesman.data) | ID | User_ID

[R] Legend position outside

Hello, Is it possible to get the legend box outside the graphic? Kind regards, Edwin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

[R] Legend position outside

Hello, Is it possible to get the legend box outside the graphic? Kind regards, Edwin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Legend position outside

Hello Gred, I try to read R-Reference, but I cant understand it. What does par (xpd=NA)? and par('usr')) mean? Thank you. Kind regards, Edwin Am Freitag, 11. April 2008 19:19:24 schrieben Sie: Yes, here is one way: plot(1:10, pch=1:2) par(xpd=NA) tmp.u - par('usr') legend

[R] Mapping file-Legend- from 2 Files/Tables

it replaces only 1 values. Kind Regards, Edwin Sendjaja __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

[R] Hex Decimal Convert

0x47780439 0x 0x 0 --- R-Table: 6565 1199047737 0 0 0 Kind Regards, Edwin __ R-help@r-project.org mailing list https

[R] Hex Decimal Convert

0x47780439 0x 0x 0 R-Table: 6565 1199047737 0 0 0 Kind Regards, Edwin __ R-help@r-project.org mailing list https

Re: [R] Hex Decimal Convert

Dear John, Thank you. Is there any possibility to get original stored number printed. Because i have another coloum like: Protocol TCP This is gonna cause probleme (as you notice before). I dont really understand what you mean with a list. i am new with R. Thanks, Edwin Am Sonntag, 6

[R] Is there simple code for this simple financial time series task?

My code below makes a data frame with columns for date, time, and price. Time on each date runs from 1 to 4. I'd like to add a new column ts$closingprice, which would have the closing price for that date. To find the closing price, I'd like to take the price in the row having the greatest time