Estou usando o pacote quickpsy em R. Gostaria de ter um parâmetro que dimensione o desvio padrão que o quickpsy calcula para a função de distribuição normal cumulativacum_normal_fun. Existe uma maneira de fazer isso? Não consigo passar os parâmetros iniciais porque o pacote quickpsy aparentemente está calculando os parâmetros (mean e std dev) usando a análise probit, portanto os parâmetros mudam a cada conjunto de dados. Esta é a minha primeira vez usando R, então talvez eu esteja perdendo algo óbvi

fit <- quickpsy(avg, Gain, nYes, n, grouping = (ID), prob = levels, xmin = 0, xmax = 2, fun = cum_normal_fun)

Eu olhei para o código fonte para quickpsy e vi que a função cum_normal_ chama apenaspnorm(), então tentei definir minha própria função como:

scaled_dev_fun <- function(x,p) { pnorm(x, p[1],p[2]) }

mas recebo a mensagem de erro:

Error in fitpsy(d, x, k, n, random, within, between, grouping, xmin, xmax,  : 
  parini (initial parameters) must be specified. 

Pelo que entendi, não posso especificar parâmetros iniciais paraparini porque isso impede o quickpsy de calculá-los por meio de análise probit, o que eu quero que ele faça. Qualquer conselho é muito apreciad