محاسبه آزمون من کندال و شیب سن در ماکرو اکسل (MAKESENS_2_0_beta)
پیش تر، به بررسی ماکروی محاسبه من کندال و شیب سن پرداخته شده بود (این لینک). در حال حاضر نسخه دوم این ماکرو منتشر شده است که به پیوست تقدیم می گردد:
این ماکرو برای تشخیص و تخمین روند در سری های زمانی از مقادیر سالیانه شیمیایی جو توسط Nordiska Ministerrådet طراحی شد. با استفاده از این ماکرو ارزیابی کیفیت هوا در کشورهای بالتیک در یک پروژه تحقیقاتی صورت پذیرفت. این روش نیاز به حداقل ۱۰ مقدار در سری زمانی دارد. جهت مطالعه بیشتر می توانید صفحه about در این فایل را مطالعه بفرمایید.
با سلام فایل متاسفانه دانلود نمیشه .
سلام و درود بررسی شد فایل قابل دانلود است روی لینک زیر کلیک کنید:
http://climatology.ir/wp-content/uploads/2019/04/MAKESENS_2_0_beta.xls
با سلام وقتتون بخیر
آزمون پتیت (Petit Test) رو با گدام نرم افزار میشه انجاکش داد؟؟ (برای تحلیل تغییرات سریع سریهای زمانی)
خیلی ممنون میشم راهنماییم کنید.
سلام و درود. حقیقت بنده تاکنون این آزمون را محاسبه نکرده ام. اما پیشنهاد می کنم لینک زیر را مشاهده کنید.
https://www.mathworks.com/matlabcentral/fileexchange/60973-pettitt-change-point-test-for-univariate-time-series-data
لازم به ذکر است ام فایل زیر به نظر برای محاسبه همان آزمون در متلب به کار می رود. این موضوع را بررسی نمایید.
<
_______________________________
% This code is used to find the change point in a univariate continuous ti
me series
% using Pettitt Test.
%
%
% The test here assumed is two-tailed test. The hypothesis are as follow:
% H (Null Hypothesis): There is no change point in the series
% H(Alternative Hypothesis): There is a change point in the series
%
% Input: univariate data series
% Output:
% The output of the answer in row wise respectively,
% loc: location of the change point in the series, index value in
% the data set
% K: Pettitt Test Statistic for two tail test
% pvalue: p-value of the test
%
%Reference: Pohlert, Thorsten. "Non-Parametric Trend Tests and Change-Point Detection." (2016).
%
function a=pettitt(data)
[m n]=size(data);
for t=2:1:m
for j=1:1:m
v(t-1,j)=sign(data(t-1,1)-data(j,1));
V(t-1)=sum(v(t-1,:));
end
end
U=cumsum(V);
loc=find(abs(U)==max(abs(U)));
K=max(abs(U));
pvalue=2*exp((-6*K^2)/(m^3+m^2));
a=[loc; K ;pvalue];
return >
سلام. شما از نرم افزار anclime هم می توانید برای این آزمون استفاده کنید
آقای دکتر ممنون از کمکتون، لطف کردین
فقط یک سوال دیگه داشتم در نرم افزار anclim که معرفی کردید، فرمت داده های ورودی رو به چه شکل در فایل txt. وارد کنیم.
من هر چقدر تلاش کردم نتونستم داده هامو به نرم افزار معرفی کنم.
سلام. به نظر اگر داده های ماه و سال را به فرمت تکس وارد کنید قابل خواندن باشد
سلام.خسته نباشید.ببخشید برای محاسبه شاخصهای اقلیمی میخام ازمون من کندال استفاده کنم.مثلا برای حدقل دمای روزانه.من داده ها روزانه رو دارم.باید میانگین بگیرم و به داده ها ماهانه تبدیلشون کنم؟
سلام. بله. ابتدا داده ها را به ماهانه تبدیل کنید. سپس روند بگیرید. داده ها را به ماه به ماه در هر سال تبدیل کنید تا بتوانید برای ۱۲ ماه روند بگیرید
باسلام.من بخوام نقشه روند رو ترسیم کنم میتونم از اکسل برای من کندال استفاده کنم یا خیر؟
سلام
داده هایتان را در اکسل جایگزین کنید و با تنظیم کردن سال و ابتدای و انتهای سال می توانید در همین فایل اکسل روند را بدست آورید.
درغیر این صورت آیا میتوانید نرم افزار دیگری غیر از R و Matlab معرفی کنید؟
می توانید به سایت نرم افزار های اقلیمی اردبیل مراجعه کنید. جناب همتی نرم افزار هایی جهت برآورد من کندال قرار داده اند
موفق باشید
سلام وقت بخیر
وقتی داده های ماهانه را جایگزین میکنم و اطلاعات را وارد میکنم بعد که روی دکمه calculated کلیک میکنم برای من چیزی محاسبه نمیکنه دلیلش چی میتونه باشه؟
نویسنده: علیرضا کربلائی
دکتری اقلیم شناسی سینوپتیک دانشگاه خوارزمی تهران