What is the R code for estimating copula parameters of BB1 with dim=2? And what's the code for gof test of BB1?

Kindly assist with R code for BB1 copula. Text books and research articles provide codes for clayton, gumbel, frank, normal and t copulas. However, I can't find code for BB1.

For example, this is the code I found:

gf <- gofCopula(normalCopula(dim = 2), as.matrix(mydata), N = 50). ## I obtained results here


But I get an error when I run:

gfBB1 <- gofCopula(BB1Copula(dim = 2), as.matrix(mydata), N = 50)


Error in BB1Copula(dim = 2) : unused argument (dim = 2)

theta <- selectedCopula$$par delta <- selectedCopula$$par2