اختبار ماتلاب كولموجوروف-سميرنوف

Question

أنا أستخدم MATLAB لتحليل بعض بيانات علم الأعصاب، وقمت بتوزيع الفواصل الزمنية بين النقاط وتناسبها بشكل أسي.بعد ذلك، أردت التحقق من هذا التوافق باستخدام اختبار Kolmogorov-Smirnov مع MATLAB.

يتم تخزين بيانات المسامير العصبية فقط في ناقل من المسامير.ال spikes المتجه هو متجه 111 × 1، حيث يكون كل إدخال هو متجه آخر.كل إدخال في ثي spikes يمثل المتجه تجربة.يختلف عدد المسامير في كل تجربة.على سبيل المثال، spikes{1} هو [1x116 double], ، مما يعني أن هناك 116 سنبلة.التالي يحتوي على 115 سنبلة، ثم 108، وهكذا.

الآن، أفهم أن اختبار kstest في MATLAB يتطلب بعض المعلمات.قمت بإدخال البيانات في أول واحد، لذلك أخذت كل الفواصل الزمنية بين النقاط وأنشأت متجه صف alldiffs الذي يخزن جميع الفواصل الزمنية بين النقاط.أرغب في ضبط CDF الخاص بي على ذلك من أجل ملاءمة الوظيفة الأسية:

test_cdf = [transpose(alldiffs), transpose(1-exp(-alldiffs*firingrate))];

لاحظ أن الأسي النظري (الذي أتناسب معه البيانات) هو r*exp(-rt) أين r هو معدل إطلاق النار.أحصل على معدل إطلاق النار حوالي 0.2.الآن، عندما أجمع كل هذا معًا، أقوم بتشغيل kstest:

[h,p] = kstest(alldiffs, 'CDF', test_cdf)

ومع ذلك، فإن النتيجة هي قيمة p حسب ترتيب 1.4455e-126.لقد حاولت إعادة test_cdf باستخدام إحدى الطرق الأخرى الموجودة في وثائق موقع Mathworks:

test_cdf = [transpose(alldiffs), cdf('exp', transpose(alldiffs), 1/firingrate)];

وهذا يعطي نفس النتيجة بالضبط!هل الملاءمة فظيعة؟لا أعرف لماذا أحصل على مثل هذه القيم p المنخفضة.الرجاء المساعدة!

أود أن أنشر صورة مناسبة، لكن ليس لدي سمعة كافية.

ملاحظة.إذا كان هناك مكان أفضل لنشر هذا، فأخبرني وسأعيد النشر.

Answer 1

فيما يلي مثال لبيانات مزيفة وطريقة أخرى لإنشاء CDF:

>> data = exprnd(.2, 100);
>> test_cdf = makedist('exp', 'mu', .2);
>> [h, p] = kstest(data, 'CDF', test_cdf)

h =

     0


p =

    0.3418

ومع ذلك، لماذا تقوم بإجراء اختبار KS؟

جميع النماذج خاطئة، وبعضها مفيد.

لا توجد خلية عصبية تماما من خلال عملية بواسون ومع وجود بيانات كافية، سيكون لديك دائمًا مؤشر ISI غير أسي بشكل ملحوظ، كما تم قياسه بواسطة اختبار KS.هذا لا يعني أنه لا يمكنك القيام بافتراض مبسط لـ ISI الأسي، اعتمادًا على الظواهر التي تحاول نموذجها.