Warning message: In sqrt(VS) : NaNs produzidos - análise modificada de Mann Kendall

Boa Tarde!

Estou trabalhando com a análise modificada de Mann Kendall para verificar tendências em séries temporais com a função mmkh do pacote "modifiedmk" Essa análise é uma adaptação para séries temporais autocorrelacionadas. Funciona bem para dados com skewness e com outliers, então não é esse o problema. Meus dados não tem zeros, números negativos ou NAs.

Após rodar a análise com várias variáveis, encontrei esse Warning em uma delas

Warning message:
In sqrt(VS) : NaNs produzidos

e o output da análise não me forneceu o parâmetro Z corrigido nem o valor de p que deveria fornecer.

O que pode ser esse erro? Procurei bastante no Google e não encontrei solução.

Acrescentando um detalhe que achei relevante:

O output da análise retorna um número referente ao tamanho efetivo da amostra e nas variáveis que apresentam o warning, esse número é negativo. Achei estranho.....

cv

Oi, Patrícia.

Pelo warning dado, alguma etapa da função faz o cálculo de raiz quadrada (o sqrt() que parece) e dá errado. Realmente é estranho que isto aconteça mesmo sem termos NAs, zeros ou números negativos nos dados.

Você pode dar um summary nos dados (ou apenas nas varáveis relevantes para a função que você está usando) e colar aqui o resultado? Minha primeira desconfiança é de que o R pode ter lido alguma variável incorretamente, então vamos começar daí e avançar mais caso o problema seja outro.

O resultado do summary é esse:

> summary(dados$chl)
  Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
 0.800   3.015   4.881   5.450   6.865  30.090

E dando str ele reconhece como numérico.

Oi, de novo, Patricia.

Procurando um pouco na internet, vi que este tipo de situação (relacionada à função sqrt retornar NaNs) costuma estar ligado ao método por trás gerar números que são muito próximos de zero e/ou geração de números negativos por erros de arrendondamento. Não encontrei uma situação com a função mmkh, mas achei situações similares com estas, por exemplo:
https://stackoverflow.com/questions/30267039/how-to-avoid-nans-due-to-rounding-errors
https://stackoverflow.com/questions/22463835/why-is-rmvnorm-function-returning-in-sqrtevvalues-nans-produced-what-i

O problema aqui, que pode ser complicado de resolver, é que os números gerados podem fazer parte de uma etapa bem "interna" dos cálculos da função... Mas eu ainda daria um chute, se ele não for comprometer muito os seus dados: seria possível você utilizar apenas os valores maiores do que 1? Vi que o menor valor nos seus dados é 0.8, e se você não for perder muitas informações descartando os valores abaixo de 1, pode valer à pena tentar.