For Loop pour tracer la courbe de survie dans Shiny

I am using the diabetic data set from the survival package in R. Within that dataset, is the column age. If I were to plot time-to-event curves comparing the patients' age, I would have way too many survival curves. I want it so that when I select the "age" button for my input, it splits the age column into 4 groups, so that when I plot the survival curve, I only have 4 curves. I keep getting an error though. Am I writing the loop correctly?

ui <- navbarPage(title = "Diabetic Retinopathy",
  ### Page 1 ###
  tabPanel(title = "Survival Curves",
  sidebarLayout(
  sidebarPanel(
  radioButtons(inputId = "comparison",
   label = "Select variable you wish to compare",
   choices = c("laser", "age", "eye", "risk"),
   selected = "laser")
  ),
  mainPanel(
  tabsetPanel(
  tabPanel(title = "Plot", 
   plotOutput(outputId = "curve")),
  tabPanel(title = "Summary", 
   verbatimTextOutput(outputId = "details"))
  )
  )
  )
  )
 )
server <- function(input, output, session) {
 ### Page 1 ###
 fit <- reactive(
 surv_fit(as.formula(paste("Surv(time, status) ~ ", 
  if(input$comparison == "age") {
  for(i in diabetic$age) {
   if(i < 15) {
   diabetic$age[i] <- "lessthan15"
   } else if(i >= 15 & i < 30) {
   diabetic$age[i]<- "fifteento30"
   } else if(i >= 30 & i < 45) {
   diabetic$age[i] <- "thirtyto45"
   } else {diabetic$age[i] <- "fortyfiveplus"}
  } 
  } else {return(input$comparison)})), data = diabetic)
 )
 
 output$curve <- renderPlot(
 ggsurvplot(fit())
 )
}
shinyApp(ui, server)

Image du message d'erreur ici


Solution du problème

dans surv_fit(as.formula(paste("Surv(time, status) ~ ",etc., deux choses se sont mélangées :

  • coller le prédicteur (laser, âge,...) à partir de la valeur deinput$comparison

  • manipuler la colonne d' âge de dataframe diabetic - qui ne renvoie cependant aucun terme de formule significatif

    • Vous pouvez ajouter une colonne age_class à vos données dès le début et l'utiliser comme option pour votre sélecteur de comparaison.

    PS : pour convertir le numérique (age) en facteur (age_class) sans if-else-skyscrapers :

    df$age_class <- cut(df$age, 
      breaks=c(-Inf, 15, 30, 45, Inf),
      labels=c("lessthan15","fifteento30",
      "thirtyto45","fortyfiveplus")
      )

    Commentaires

    Enregistrer un commentaire

    Posts les plus consultés de ce blog

    Comment signer la connexion OKEx API version 5 avec websockets ?

    Bien que cette question ait reçu une réponse pour les versions précédentes de l'API OKEx utilisant REST, elle ne l'a pas été pour la dernière version 5 de l'API utilisant des websockets. Les docs sont ici. Je reçois l'erreur suivante {"event":"error","msg":"Invalid sign","code":"60007"} , il doit donc y avoir un problème avec l'algorithme de chaîne de signature, mais je n'arrive pas à identifier où je fais une erreur. import hmac import json import time import hashlib import asyncio import websockets passphrase = "XXXX" secret_key = b"XXXX" api_key = "XXXX" timestamp = int(time.time()) print("timestamp: " + str(timestamp)) sign = str(timestamp) + 'GET' + '/users/self/verify' total_params = bytes(sign, encoding= 'utf-8') signature = hmac.new(bytes(secret_key, encoding= 'utf-8'), total_params, digestmod=hashlib.sha256).digest() sign...
    Lire la suite

    La fonction GCP Cloud pour écrire des données dans BigQuery s'exécute avec succès, mais les données n'apparaissent pas dans la table BigQuery

    Image
    J'exécute la fonction cloud suivante. Il s'exécute avec succès et indique que les données ont été chargées dans la table. Mais lorsque j'interroge BigQuery, aucune donnée n'a été ajoutée. Je ne reçois aucune erreur et aucune indication que cela ne fonctionne pas. from google.cloud import bigquery import pandas as pd def download_data(event, context): df = pd.read_csv('https://covid.ourworldindata.org/data/ecdc/full_data.csv') # Create an empty list Row_list =[] # Iterate over each row for index, rows in df.iterrows(): # Create list for the current row my_list =[rows.date, rows.location, rows.new_cases, rows.new_deaths, rows.total_cases, rows.total_deaths] #print(my_list) # append the list to the final list Row_list.append(my_list) ## Get Biq Query Set up client = bigquery.Client() table_id = "<project_name>.raw.daily_load" table = client.get_table(table_id) print(client) print(table_id) print(table) errors = client.insert_rows(...
    Lire la suite

    Erreur Symfony : "Une exception a été levée lors du rendu d'un modèle"

    j'essaie d'inclure <link href="{{ asset('css/mystyle.css') }}" rel="stylesheet"/> dans mon fichier twig et le rendre. Mais cela me donne une exception PHP non interceptée Twig_Error_Runtime : "Une exception a été levée lors du rendu d'un modèle ("Le fichier manifeste de l'actif "/../public/build/manifest.json" n'existe pas.") ça marche quand j'utilise link href="css/mystyle.css" rel="stylesheet"/> . Voici mon contrôleur: <?php namespace App\Controller; use Symfony\Component\Routing\Annotation\Route; use Symfony\Bundle\FrameworkBundle\Controller\AbstractController; use Symfony\Component\HttpFoundation\Response; class WelcomeController extends AbstractController { /** * @Route("/", name="welcome") */ public function index() { return $this->render('welcome/index.html.twig'); } } Voici mon modèle de brindille : <!doctype html> ...
    Lire la suite