Чтобы пояснить мой комментарий выше, вы можете отредактировать функцию в R следующим образом:
qqcompValues <- edit(qqcomp)
В редакторе замените строку 111 на следующую:
return(data.frame(x=fittedquant, y=sdata))
Обратите внимание, что в приведенном выше редакторе предполагается, что вы используете параметр функции по умолчанию plotstyle = "graphics".
Затем вы можете получить значения QQ следующим образом:
qqValues <- qqcompValues(lognormalfit)
qqValues
x y
1 1026674 1050000
2 1158492 1100000
3 1247944 1230000
4 1327616 1300000
5 1407939 1450000
6 1497825 1459785
7 1613479 1654000
8 1820639 1888000