Create a Weibull Probability Paper Plot in R

Multi tool use
Multi tool use


Create a Weibull Probability Paper Plot in R



I am trying to create a base plot for a weibull probability plot. I have been exploring the survival package in R but have not found the exact tool I need yet. So far I have been able to recreate the plot I need by hand (which is terrible - tons of hard coding, it is not flexible at all and looks terrible).


survival



My guess is that there is a package to create this base graph, with the grid lines, I just haven't discovered it yet.



This image is what I need to graph, it is called "Weibull probability plotting paper"



weibull paper



Is there a way in ggplot2 or within the survival package (or anything other that base R graphics) to create this plot? From there, I can overlay the points.


ggplot2


survival



Please keep in mind the graphics I create in the end will need to be compatible with Shiny. Thank you.


Shiny



For reference this is the ugly plot I have managed to create using ggplot2 by hand.


ggplot2



current_plot



Also, please note I cannot provide a reproducible code example since this is a very complex problem, I am looking for a basic reproducible base plot.



UPDATE



I am looking for a what to produce this plot, here is some examples of how I created mine by hand, please note it is not fully reproducible.



Here is how I set up the y-axis tick marks:


yticks <- c(log(-log(1-0.0001)), log(-log(1-0.001)), log(-log(1-0.005)), log(-log(1-0.010)), log(-log(1-0.05)), log(-log(1-0.10)),
log(-log(1-0.50)), log(-log(1-0.90)), log(-log(1-0.99)))



Here are the y-labels:


ylabs <- c('0.01','0.1','0.5','1','5','10','50','90','99')



Here is how to create the minor grid lines:


yminorticks <- c(
log(-log(1-0.00001)),log(-log(1-0.00002)),log(-log(1-0.00003)),log(-log(1-0.00004)),log(-log(1-0.00005)),log(-log(1-0.00006)),log(-log(1-0.00007)),log(-log(1-0.00008)),log(-log(1-0.00009)),
log(-log(1-0.0001 )),log(-log(1-0.0002)),log(-log(1-0.0003)),log(-log(1-0.0004)),log(-log(1-0.0005)),log(-log(1-0.0006)),log(-log(1-0.0007)),log(-log(1-0.0008)),log(-log(1-0.0009)),
log(-log(1-0.001)),log(-log(1-0.002)),log(-log(1-0.003)),log(-log(1-0.004)),log(-log(1-0.005)),log(-log(1-0.006)),log(-log(1-0.007)),log(-log(1-0.008)),log(-log(1-0.009)),
log(-log(1-0.01)),log(-log(1-0.02)),log(-log(1-0.03)),log(-log(1-0.04)),log(-log(1-0.05)),log(-log(1-0.06)),log(-log(1-0.07)),log(-log(1-0.08)),log(-log(1-0.09)),
log(-log(1-0.10)),log(-log(1-0.20)),log(-log(1-0.30)),log(-log(1-0.40)),log(-log(1-0.50)),log(-log(1-0.60)),log(-log(1-0.70)),log(-log(1-0.80)),log(-log(1-0.90)))





You just want a plot with log scales? When asking for help, you should include a simple reproducible example with sample input and desired output that can be used to test and verify possible solutions.
– MrFlick
Jul 2 at 15:38





I cannot provide a reproducible example since I do not have a way of creating this plot. That is my question, how to create this base plot in a simple and reproducible manner. Thanks.
– user9842212
Jul 2 at 15:40





Just an aside on your minor grid line code: log, and subtraction are all vectorized. vectorization means you can do, e.g., log(-log(1 - seq(0.01, 0.09, by = 0.01))) instead of c(log(-log(1-0.01)),log(-log(1-0.02)),log(-log(1-0.03)),log(-log(1-0.04)), ... )
– Gregor
Jul 2 at 16:05



log


log(-log(1 - seq(0.01, 0.09, by = 0.01)))


c(log(-log(1-0.01)),log(-log(1-0.02)),log(-log(1-0.03)),log(-log(1-0.04)), ... )





You can generate yminorticks with yminorticks <- c(log(-log(1-outer(1:9, 10^(-5:-1))))) rather than writing all that out. Though you probably want to be using log10 rather than log
– MrFlick
Jul 2 at 16:06



yminorticks


yminorticks <- c(log(-log(1-outer(1:9, 10^(-5:-1)))))


log10


log





I think what you want to do for ggplot2 is define your own transformation with the scales package and use that with the trans= argument for scale_x_continuous. Here's an example where I do something similar to make a sign-restoring cox-box transformation github.com/ellisp/frs-r-package/blob/master/pkg/R/…
– Peter Ellis
Jul 3 at 2:55


trans=


scale_x_continuous









By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

WOy4 zyP 5nxRRH9xfWNZlmd
-
bYBdCeRAdkyHSh pV Vv3tMm2nFrA 3gIica6GRcYpMAKaFCjuJLqhAj4ZTSNoIASR7ibZmQJtMJ ToRUGpg,5E8

Popular posts from this blog

PHP contact form sending but not receiving emails

PHP contact form sending but not receiving emails I have created a contact form using PHP to send the data. The form seems to send fine as there are no errors and success message appears, however I am not receving the emails into the designated inbox. I am using PHPMailer as I originally tried just using the php 'mail' command which I now understand is a bit hit and miss. I cannot see why I would not be receving the emails so I would be very grateful for any help that could be given. I'm quite new to PHP so please be patient with me :) <?php use PHPMailerPHPMailerPHPMailer; use PHPMailerPHPMailerException; $msg = ""; if (isset($_POST['submit'])) { require 'phpmailer/src/Exception.php'; require 'phpmailer/src/PHPMailer.php'; require 'phpmailer/src/SMTP.php'; function sendemail ($to, $from, $body) { $mail = new PHPMailer(true); $mail->setFrom($from); $mail->addAddress($to); $ma...
Read more

Do graphics cards have individual ID by which single devices can be distinguished?

Do graphics cards have individual ID by which single devices can be distinguished? I have several graphics cards of the same manufacturer and same model and I want to distinguish them not only by the pcie socket they reside in but by some individual value I can read from the hardware. Is there such value and if so what would it be called and how could I access it? Right now I'm using modern nvidia cards (10xx series) and nvidia inspector.but I'd prefer if there was any value that was also applicable to amd. I have programming and scripting knowledge. I am running Windows 10 1 Answer 1 This is a powershell example to pull information from video cards. The DeviceID is supposed to be unique but only if the manufacturer supports Windows Display Driver Model (WDDM). (I copied this directly from docs.microsoft.com) $strComputer = "." $colItems = get-wmiobject -class "Win32_VideoC...
Read more

Create weekly swift ios local notifications

Create weekly swift ios local notifications I can create daily notifications but not weekly ones using swift 3 ios 10 User Notifications framework. Daily works fine but when I add in the .weekday parameter it doesn't send me a notification. The code is: for i in 1 ... 7 { let center = UNUserNotificationCenter.current(); let content = UNMutableNotificationContent(); content.title = "Title"; content.body = "content"; content.sound = UNNotificationSound.default(); var dateComponents = DateComponents(); dateComponents.weekday = i; // hours is 00 to 11 in string format if (daynight == "am") { dateComponents.hour = Int(hours); } else { dateComponents.hour = Int(hours)! + 12; } // mins is 00 to 59 in string format dateCom...
Read more