Good afternoon,
Can someone help me with this problem? I am trying to create a linear regression with this data and in some variables I have the following problem:
Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
Can you please help me? Thank you very much for your help.
Version:1.0 StartHTML:0000000107 EndHTML:0000038820 StartFragment:0000000127 EndFragment:0000038802
library(tidyverse)
#> Warning: package 'tidyverse' was built under R version 3.6.3
#> Warning: package 'ggplot2' was built under R version 3.6.3
#> Warning: package 'tibble' was built under R version 3.6.3
#> Warning: package 'tidyr' was built under R version 3.6.3
#> Warning: package 'dplyr' was built under R version 3.6.3
#> Warning: package 'forcats' was built under R version 3.6.3
library(dplyr)
library(tibble)
library(psych)
#> Warning: package 'psych' was built under R version 3.6.3
#>
#> Attaching package: 'psych'
#> The following objects are masked from 'package:ggplot2':
#>
#> %+%, alpha
library(flextable)
#> Warning: package 'flextable' was built under R version 3.6.3
#>
#> Attaching package: 'flextable'
#> The following object is masked from 'package:purrr':
#>
#> compose
library(Gmisc)
#> Warning: package 'Gmisc' was built under R version 3.6.3
#> Loading required package: Rcpp
#> Loading required package: htmlTable
#> Warning: package 'htmlTable' was built under R version 3.6.3
library(Hmisc)
#> Warning: package 'Hmisc' was built under R version 3.6.3
#> Loading required package: lattice
#> Loading required package: survival
#> Loading required package: Formula
#>
#> Attaching package: 'Hmisc'
#> The following object is masked from 'package:psych':
#>
#> describe
#> The following objects are masked from 'package:dplyr':
#>
#> src, summarize
#> The following objects are masked from 'package:base':
#>
#> format.pval, units
library(kableExtra)
#> Warning: package 'kableExtra' was built under R version 3.6.3
#>
#> Attaching package: 'kableExtra'
#> The following objects are masked from 'package:flextable':
#>
#> as_image, footnote
#> The following object is masked from 'package:dplyr':
#>
#> group_rows
library(knitr)
#> Warning: package 'knitr' was built under R version 3.6.3
library(prettydoc)
#> Warning: package 'prettydoc' was built under R version 3.6.3
library(gtsummary)
#> Warning: package 'gtsummary' was built under R version 3.6.3
#>
#> Attaching package: 'gtsummary'
#> The following object is masked from 'package:flextable':
#>
#> as_flextable
library(eeptools)
#> Warning: package 'eeptools' was built under R version 3.6.3
library(rstudioapi)
#> Warning: package 'rstudioapi' was built under R version 3.6.3
library(htmlTable)
library(lubridate)
#>
#> Attaching package: 'lubridate'
#> The following object is masked from 'package:base':
#>
#> date
library(nnet)
library(questionr)
#> Warning: package 'questionr' was built under R version 3.6.3
#>
#> Attaching package: 'questionr'
#> The following objects are masked from 'package:Hmisc':
#>
#> describe, wtd.mean, wtd.table, wtd.var
#> The following object is masked from 'package:psych':
#>
#> describe
library(sjPlot)
#> Warning: package 'sjPlot' was built under R version 3.6.3
library(readr)
de <- read_csv("C:/Users/juanp/Desktop/de.csv",
col_types = cols(status.dead = col_factor(levels = c("1",
"0")), status.icu = col_factor(levels = c("1",
"0"))))
#> Warning: Missing column names filled in: 'X1' [1]
structure(de)
#> # A tibble: 981 x 9
#> X1 PATIENT.NUMBER status status.icu status.dead lockdown_date Age_category
#> <dbl> <dbl> <chr> <fct> <fct> <chr> <chr>
#> 1 1 0 Disch~ 0 0 Before lockd~ (50 to 59 y~
#> 2 2 1 Disch~ 0 0 After lockdo~ (50 to 59 y~
#> 3 3 2 Disch~ 0 0 After lockdo~ [18 to 40 y~
#> 4 4 3 Disch~ 0 0 After lockdo~ Above 80 ye~
#> 5 5 3 Disch~ 0 0 After lockdo~ Above 80 ye~
#> 6 6 4 Disch~ 0 0 After lockdo~ [18 to 40 y~
#> 7 7 5 Disch~ 0 0 After lockdo~ (50 to 59 y~
#> 8 8 6 Disch~ 0 0 After lockdo~ [18 to 40 y~
#> 9 9 7 Disch~ 0 0 After lockdo~ (50 to 59 y~
#> 10 10 8 Disch~ 0 0 After lockdo~ (60 to 69 y~
#> # ... with 971 more rows, and 2 more variables: OBESITY <dbl>, morbidity <chr>
t3.2 <-
glm(de$status.dead ~ de$Age_category , data= de, family = binomial) %>%
tbl_regression(exponentiate = TRUE) %>%
bold_p(t = 0.10) %>%
bold_labels() %>%
italicize_levels()
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
#> Warning in regularize.values(x, y, ties, missing(ties)): collapsing to unique
#> 'x' values
t3.2
| Characteristic |
OR 1 |
95% CI 1 |
p-value |
| de$Age_category |
|
|
|
| (40 to 49 years] |
— |
— |
|
| (50 to 59 years] |
9.18 |
1.54, 175 |
0.042 |
| (60 to 69 years] |
4.41 |
0.99, 30.5 |
0.073 |
| (70 to 79 years] |
0.55 |
0.19, 1.45 |
0.2 |
| [18 to 40 years) |
6.09 |
1.38, 42.1 |
0.029 |
| Above 80 years |
0.10 |
0.04, 0.22 |
<0.001 |
| Under 18 years |
958247 |
0.00, 3558377584379648918482064260406646080066682662420688868642226602864626064426646442264480046204228684402864062280648644 |
> 0.9 > > |
|
1 OR = Odds Ratio, CI = Confidence Interval |
|
|
|