Eu medi dois parâmetros (carbono orgânico dissolvido DOC = y, e descarga = x). Quando essas duas variáveis são plotadas uma contra a outra, obtemos um loop de histerese (veja o exemplo de código e a figura).
Agora, para uma análise mais aprofundada, eu gostaria de determinar a área desse loop histérico. Eu descobri que isso pode ser feito usando o método de dardo em Monte Carlo. Este método diz que a área de uma área desconhecida é proporcional à área de um retângulo conhecido vezes os acertos no campo interno (o loop).
Meu problema agora é: como resolver o problema interno / externo usando R. Como posso desenhar um retangular com uma área conhecida e como posso destacar os acertos aleatórios dentro e fora do loop histerético?
Observe que estou aberto a qualquer outro método ...
Pesquisei em vários sites estatísticos, mas não consegui encontrar uma resposta. Qualquer ajuda direta ou ligação a outros sites / postagens é muito apreciada.
Data <- read.table("http://dl.dropbox.com/u/2108381/DOC_Q_hystersis.txt", sep = ";",
header = T)
head(Data)
plot(Data$Q, Data$DOC, type = "o", xlab = "Discharge (m3 s-1)", ylab = "DOC (mg C l-1)",
main = "Hystersis loop of the C/Q relationship")
fonte
Uma possibilidade seria esta: parece-me que o ciclo de histerese deve ser convexo, certo? Assim, pode-se gerar pontos e testar para cada ponto se faz parte do casco convexo da união do seu conjunto de dados e o ponto aleatório - se sim, está fora do conjunto de dados original. Para acelerar as coisas, podemos trabalhar com um subconjunto do conjunto de dados original, ou seja, os pontos que compõem seu próprio casco convexo.
Obviamente, a maneira R de fazer as coisas não usaria meu
for
loop, mas isso é fácil de entender e não é muito lento. Recebo uma área estimada de 0,703.EDIT: é claro, isso depende da suposta convexidade do relacionamento. Por exemplo, parece haver uma parte não-convexa no canto inferior direito. Em princípio, poderíamos estimar Monte-Carlo da área dessa área da mesma maneira e subtraí-la da estimativa original da área.
fonte