diff --git a/.gitignore b/.gitignore index 8f9ec2c..f672d6b 100644 --- a/.gitignore +++ b/.gitignore @@ -6,4 +6,6 @@ fitdistrplus.Rproj inst/doc/ revdep fitdistrplus-github.Rproj -share/todolist.html \ No newline at end of file +share/todolist.html +AUTHORS +LICENSE diff --git a/share/pb-exists-envir-primarycensored.R b/share/pb-exists-envir-primarycensored.R new file mode 100644 index 0000000..cb90b92 --- /dev/null +++ b/share/pb-exists-envir-primarycensored.R @@ -0,0 +1,47 @@ + +install.packages( + "primarycensored", + repos = "https://epinowcast.r-universe.dev" +) + +library(primarycensored) +?fitdistdoublecens + +set.seed(123) +n <- 1000 +true_mean <- 5 +true_sd <- 2 +pwindow <- 2 +swindow <- 2 +D <- 10 +samples <- rprimarycensored( + n, rnorm, + mean = true_mean, sd = true_sd, + pwindow = pwindow, swindow = swindow, D = D +) +#etrange comme tirage +table(samples) + +delay_data <- data.frame( + left = samples, + right = samples + swindow +) + +fit_norm <- fitdistdoublecens( + delay_data, + distr = "norm", + start = list(mean = 0, sd = 1), + D = D, pwindow = pwindow +) + +head(delay_data) + +curve(pprimarycensored(x, pnorm, mean = 0, sd = 1, D=D, pwindow=pwindow), from=-5, to=10) +lines(ecdf(samples)) + +#c'est tres voisin de fitdistcens... +fitdistcens(delay_data, "norm") + +fitdistcens(delay_data, dnorm) + +