Re: [R] How to decrease size of points?

2020-09-30 Thread Medic
Yes, I just wanted to decrease the STARTING point size! Sorry for not being able to formulate Thank you very much, Rui and Avi! That's all for now. I need to comprehend information. Medic __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and

Re: [R] How to decrease size of points?

2020-09-30 Thread Medic
№1 Medic: The code works as I want, but the points (circles) on the plot are too big. How to decrease them? Where to insert (for instance) size = 0.8 for points (circles) on plot? p1 <- p + geom_point(aes(size = Stage), alpha = 1/3) + xlab ("X") + ylab("Y") + geom_smooth()

Re: [R] How to decrease size of points?

2020-09-30 Thread Medic
The code works as I want, but the points (circles) on the plot are too big. How to decrease them? Where to insert (for instance) size = 0.8 for points (circles) on plot? p1 <- p + geom_point(aes(size = Stage), alpha = 1/3) + xlab ("X") + ylab("Y") + geom_smooth() Stage is factor, x and y - continuo

[R] How to decrease size of points?

2020-09-30 Thread Medic
The code works as I want, but the points are too big. How to decrease them? (Where to insert: size = 0.8?) p1 <- p + geom_point(aes(size = Stage), alpha = 1/3) + xlab ("X") + ylab("Y") + geom_smooth() __ R-help@r-project.org mailing list -- To UNSUBSCRIB

Re: [R] Date format

2020-05-10 Thread Medic
Many Thanks!!! > cpolw...@chemo.org.uk: > Your X axis is plotting mydata not date? > Use aes(x=date __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://

Re: [R] Date format

2020-05-10 Thread Medic
I took a SAMPLE CODE (for Connected scatterplot) from the R gallery and applied to MY DATA, but got: "Don't know how to automatically pick scale for object ..." P.S. 1) R ver. 4.0 (Yes, Jeff); 2) Attached: mydata_dput (1 КБ) SAMPLE CODE library(ggplot2) library(dplyr) library(hrbrthemes) data <-

[R] Date format

2020-05-09 Thread Medic
I took a SAMPLE CODE (for Connected scatterplot) from the R gallery and applied to MY DATA, but got: "Error in as.Date.numeric(mydata$date) : 'origin' must be supplied". P.S. I can not understand ?as.Date() SAMPLE CODE library(ggplot2) library(dplyr) library(hrbrthemes) data <- read.table("https:

Re: [R] Survival analysis

2020-04-26 Thread Medic
cpolwart wrote: > the event can only happen once > 1. Yes, I'm not interested in repeat events. 2. I just had an interest in learning the code for survival analysis in case there are: event, left censored (not left truncating) and right censored, nothing more! But I no longer want to burden anyone

Re: [R] Survival analysis

2020-04-26 Thread Medic
Very grateful for the all comments! My data contains: • left censored • right censored • events (interval censored does not contain!) (P.S. I understood, that the code with "type = 'left'" is not suitable, because is ONLY for left-censored.) I wanted to get the appropriate code for my so mix dat

Re: [R] Survival analysis

2020-04-17 Thread Medic
On 2020-04-17 20:06, Medic wrote: > I can't understand how to do a survival analysis (?Surv ()) when some > event occurred before the start of observation (left censored). If I > understand correctly, there are two methods. I chose a method with: 1) > time from the start of treat

[R] Survival analysis

2020-04-17 Thread Medic
I can't understand how to do a survival analysis (?Surv ()) when some event occurred before the start of observation (left censored). If I understand correctly, there are two methods. I chose a method with: 1) time from the start of treatment to the event and 2) the indicator of the event. I did (i

[R] Split

2019-12-23 Thread Medic
Dear Burt, you gave a very elegant solution. Many thanks! Jeff, I understand your solution, thank you very much for your time! Colleague Milos, patronage is exactly what I need. I hope for your further guidance! (Rcmdr is not enough for some purposes.) Dear Ivan, I used your solution! It's the most

[R] Split

2019-12-23 Thread Medic
I have mydata$var #this is ONE group of patients And I would like to get median and ICR of mydata$var. How can I get this? With summary (mydata$var)! Ok! And now I would like to get THE SAME, but for TWO group: male and female (which are contained in the group mydata$var) How can I get this? F

[R] Split

2019-12-23 Thread Medic
I have mydata$var and I have mydata$group #two group I would like to split mydata$var by mydata$group #to get var1 and var2 And then get summary (var1, var2) #this is my finite aim How to encode it all? __ R-help@r-project.org mailing list -- To UNS

[R] choose randomly

2019-12-18 Thread Medic
rt Thank you, Jim, you always come to the rescue! Jim Lemon mydata$Temperature[sample(1:N,N) should do the trick. You will just get a pseudo-randomly shuffled set of the same values. Jim Medic wrote: > > Variable temperature: > > mydata$temperature > > has N value

[R] choose randomly

2019-12-18 Thread Medic
Variable temperature: mydata$temperature has N values. With what code to сhoice (without return) n values from them RANDOMLY? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] Creating map in R

2019-11-22 Thread Medic
Thank you very much! > From: Michael Hannon > ... following might be useful in this case: > https://stat.ethz.ch/mailman/listinfo/r-sig-geo __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

[R] Creating map in R

2019-11-21 Thread Medic
Creating map in R. I have all the files (with different extensions) for the country I need. The problem is that I don’t know how to use them in R. I began according to information on the Internet, and immediately ran into difficulties (see below): `library("tmaptools") geo <- read_shape("Rom.shp",

[R] Read shp file

2019-11-21 Thread Medic
Help me. pls, to read .shp file. `library("tmaptools") geo <- read_shape("Rom.shp", as.sf = TRUE) This function is deprecated and has been migrated to github.com/mtennekes/oldtmaptools` I have to turn to another function, but I get an unclear message `library(raster) geo <- shapefile ("Rus.shp")

[R] pairs.panels ()

2019-11-20 Thread Medic
Hi, Sarah, Thank you not just for the answer, but for the teaching (!) answer! Yes, cex.labels is what I need! __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

[R] pairs.panels ()

2019-11-19 Thread Medic
Point me, please, the parameter for changing the font size for the variable name in pairs.panels () __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://w

Re: [R] Change Y-axis labels

2019-09-05 Thread Medic
For David Winsemius. As always, You help out! Immensely grateful! __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.ht

[R] Change Y-axis labels

2019-09-04 Thread Medic
The Y scale is divided (by default) as: 0.0 ... 0.2 ... 0.4 ... 0.6 ... 0.8 ...1.0 But I would like so: 0 ... 20 ... 40 ... 60 ... 80... 100 (with rotating axis labels) When I use par function (marked as comment here) it turns out correctly for ONLY ONE picture?! Help me, please. (This is the code

Re: [R] plot (cox)

2019-04-15 Thread Medic
in your plot() statement (see ?par for > > > details), and then use axis() to make anything you'd like (see ?axis > > > for details). > > -- > > >> Medic wrote: > > >> In this code: > > >> plot (cox, col

Re: [R] plot (cox)

2019-04-15 Thread Medic
Thanks, but too hard for me Sarah Goslee : > You can presumably use xaxt="n" in your plot() statement (see ?par for > details), and then use axis() to make anything you'd like (see ?axis > for details). ------ >> Medic wrote: >>

[R] plot (cox)

2019-04-15 Thread Medic
In this code: plot (cox, col=1:2, xscale=1, xlab="OS", ylab="Probability") the X scale is divided (by default) as: 0 ... 50 ... 100 ... 150 ... 200 And I would like so: 0 ... 12 ... 24 ... 36 ... 48. I looked ?plot(cox), but did not understand what argument is responsible for this. Pls, he

[R] Thanks! Re: Kaplan-Meier plot

2019-01-17 Thread Medic
Bert Gunter: "Have you consulted ?plot.survfit ? " Marc Schwartz "The 'mark.time' argument for plot.survfit() is FALSE by default." Great thanks, Bert, for explanation in which documentation to view information! Thank you very much, Marc! Yes, specification required: plot (km, mark.time=TRUE) for

[R] Kaplan-Meier plot

2019-01-17 Thread Medic
According to the guidelines (if I'm not mistaken), the code below is sufficient (without any specification) to give Kaplan-Meier curves with censored data markings on Kaplan-Meier curves. But in my case censored data don't appears on the curves?! library(survival) mydata<-read.csv (file="C:/mydata

Re: [R] SURVDIFF()

2019-01-12 Thread Medic
Dear David, you are (as always) come to the rescue! With your wonderful pedagogical talent of explanation! Many thanks for the support! My warmest and sincere wishes for the new year! __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see h

[R] SURVDIFF()

2019-01-12 Thread Medic
How to note (in code) a few (!) adjusting covariates for cox regression. I had an example for one covariate, and tried (according to my own understanding) two variantes of code (pls, see below), and got ... a different p-value. What is the right code? Many thanks!!! (1) survdiff (Surv(survt,statu

[R] THANKS! Re: Installation...

2018-12-24 Thread Medic
1. Dear John, thank you for the important point! 2. Dear Duncan, thank you for the clear clarification! __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide h

[R] Installation Re: Problems with library...

2018-12-24 Thread Medic
automatically) 2) installing all packages only in RStudio (R will take them automatically) 3) installing all packages in R and in RStudio (i.e. in each of them) Thanks, Medic __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https

[R] Problems with library(survival)

2018-12-23 Thread Medic
> install.packages("survival") package ‘survival’ successfully unpacked and MD5 sums checked > library(survival) Error: package or namespace load failed for ‘survival’ in get(Info[i, 1], envir = env): lazy-load database 'C:/Program Files/R/R-3.5.1/library/lattice/R/lattice.rdb' is corrupt In addit

Re: [R] HISTOGRAM

2018-11-10 Thread Medic
Rui Barradas, thank you for your prompt response, your code will be useful to me in the future! Rick Bilonick ("your data appear to be categorical"), thank you very much for your comment (I would have to more correctly express my task). JIM Lemon, THANKS!!! THIS IS EXACTLY what I needed! ___

[R] HISTOGRAM

2018-11-09 Thread Medic
What would be the correct code (simplest version) (without gplot()) for histogram (with 7 bars), which would include 7 names of bars under the X-axis. The data are: name number ds6277 lk 24375 ax46049 dd70656 az216544 df 220620 gh641827 (I'm attaching mydata.r, making

[R] cox model

2018-11-03 Thread Medic
I need a R-code for a situation that is well described in the sas help. I would be very grateful for the help! "Time-dependent variables can be used to model the effects of subjects transferring from one treatment group to another. One example of the need for such strategies is the Stanford heart t

[R] survival analysis question

2018-10-31 Thread Medic via R-help
When I run the code from the book Kleinbaum "Survival analysis" (7. Running an extended cox model. Page 646) I got an ERROR. What do you think is the reason? I attach a file with few lines of dataset "addicts-dput.r". > library(cmprsk) > addicts = read.csv ("addicts.csv")     id clinic status su