You can replicate if you’d like. Go to https://vaers.hhs.gov/data/datasets.html?
Then download “All Years Data” and unzip.

Here is the R script:

v21 = read.csv(“2021VAERSVAX.csv”)
v20 = read.csv(“2020VAERSVAX.csv”)
a21 = read.csv(“2021VAERSDATA.csv”)
a20 = read.csv(“2020VAERSDATA.csv”)
v = rbind(v21, v20)
a = rbind(a21, a20)
b = merge(a,v, on=’VAERS_ID’)
length(unique(b[b$DIED == ‘Y’ & b$VAX_TYPE == ‘COVID19’ & b$VAX_MANU == ‘MODERNA’ & b$VAX_DOSE_SERIES == 1 & b$NUMDAYS < 30,]$VAERS_ID))
length(unique(b[b$DIED == ‘Y’ & b$VAX_TYPE == ‘COVID19’ & b$VAX_MANU == ‘MODERNA’ & b$VAX_DOSE_SERIES == 2 & b$NUMDAYS < 30,]$VAERS_ID))
length(unique(b[b$DIED == ‘Y’ & b$VAX_TYPE == ‘COVID19’ & b$VAX_MANU == ‘PFIZER\\BIONTECH’ & b$VAX_DOSE_SERIES == 1 & b$NUMDAYS < 30,]$VAERS_ID))
length(unique(b[b$DIED == ‘Y’ & b$VAX_TYPE == ‘COVID19’ & b$VAX_MANU == ‘PFIZER\\BIONTECH’ & b$VAX_DOSE_SERIES == 2 & b$NUMDAYS < 30,]$VAERS_ID))

#subtract those who died from dose 2 from the dose 1 totals as there is a #vaers entry for dose 1 for all dose 2 deaths

x = c(30,60,100,200)
y = c( 394/113, 710/113, 545/75.74, 767/75.74)

vnames = c(“Pfizer Dose 1”, “Pfizer Dose 2”, “Moderna Dose 1”, “Moderna Dose 2”)

plot(x,y, xlab=”Cumulative mRNA dose (μg)”, ylab=”Fatality Rate per 1m doses < 30 days”, main=”Cumulative mRNA Dose vs. Fatalities”)
text(x,y,labels=vnames,cex=1,pos=3)

lx = log(x)

plot(lx,y, xlab=”Cumulative mRNA dose (μg)”, ylab=”Fatality Rate per 1m doses < 30 days”, main=”log Cumulative mRNA Dose vs. Fatalities”)

fit <- lm(y~lx)
abline(fit, col=2)