মিশ্র মডেলগুলির জন্য প্যারামেট্রিক, সেমিপ্রেমেট্রিক এবং ননপ্যারমেট্রিক বুটস্ট্র্যাপিং


9

নিম্নলিখিত গ্রাফ্ট এই নিবন্ধ থেকে নেওয়া হয়েছে । আমি বুটস্ট্র্যাপে নবাগত এবং R bootপ্যাকেজের সাথে রৈখিক মিশ্র মডেলের জন্য প্যারামিমেট্রিক, সেমিপ্রায়মেট্রিক এবং ননপ্যারমেট্রিক বুটস্ট্র্যাপিং বুটস্ট্র্যাপিং বাস্তবায়নের চেষ্টা করছি ।

এখানে চিত্র বর্ণনা লিখুন এখানে চিত্র বর্ণনা লিখুন এখানে চিত্র বর্ণনা লিখুন এখানে চিত্র বর্ণনা লিখুন এখানে চিত্র বর্ণনা লিখুন

আর কোড

আমার Rকোডটি এখানে :

library(SASmixed)
library(lme4)
library(boot)

fm1Cult <- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=Cultivation)
fixef(fm1Cult)


boot.fn <- function(data, indices){
 data <- data[indices, ]
 mod <- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=data)
 fixef(mod)
 }

set.seed(12345)
Out <- boot(data=Cultivation, statistic=boot.fn, R=99)
Out

প্রশ্নাবলি

  1. bootপ্যাকেজ সহ মিশ্র মডেলগুলির জন্য কীভাবে প্যারামেট্রিক, সেমিপ্রেমেট্রিক এবং ননপ্যারমেট্রিক বুটস্ট্র্যাপিং করবেন ?
  2. আমার ধারণা আমি আমার কোডে মিশ্র মডেলটির জন্য ননপ্যারমেট্রিক বুটস্ট্র্যাপিং করছি।

আমি এই স্লাইডগুলি পেয়েছি কিন্তু আর প্যাকেজটি পেলাম না merBoot। আমি এই প্যাকেজটি পেতে পারি যে কোনও ধারণা। যে কোনও সাহায্যের প্রশংসা করা হবে। আপনার সাহায্য এবং সময় জন্য অগ্রিম ধন্যবাদ।

উত্তর:


8

মিশ্র রৈখিক মডেলগুলিতে বুটস্ট্র্যাপিং হ'ল রিগ্রেশন-এ বুটস্ট্র্যাপিংয়ের মতো তবে আপনার জটিলতা রয়েছে যে প্রভাবগুলি স্থির এবং এলোমেলোভাবে বিভক্ত। প্যারামেট্রিক বুটস্ট্র্যাপ করার জন্য, আপনি প্যারামেট্রিক মডেলটিকে ডেটাতে ফিট করে, মডেলের অবশিষ্টাংশগুলি গণনা করেন, অবশিষ্টাংশগুলি বুটস্ট্র্যাপ করেন, বুটস্ট্র্যাপের অবশিষ্টাংশগুলি নিয়ে যান এবং তথ্যের জন্য বুটস্ট্র্যাপের নমুনা পেতে ফিটেড মডেলটিতে যুক্ত করে তারপরে মডেলটি ফিট করেন বুটস্ট্র্যাপের ডেটাতে বুটস্ট্র্যাপের নমুনার পরামিতি অনুমানগুলি পেতে। আপনি মূল বাকী অংশগুলি আবার বুটস্ট্র্যাপ করে প্রক্রিয়াটির পুনরাবৃত্তি করেন এবং তারপরে পরামিতিগুলির আরেকটি বুটস্ট্র্যাপ নমুনা অনুমান করার জন্য পদ্ধতিতে অন্যান্য পদক্ষেপগুলি পুনরাবৃত্তি করেন। ননপ্যারমেট্রিক বুটস্ট্র্যাপের জন্য, আপনি প্রতিক্রিয়া এবং কোভারিয়েট মানগুলির ভেক্টর তৈরি করেন এবং বুটস্ট্র্যাপ নমুনার জন্য ভেক্টরগুলির নির্বাচন বুটস্ট্র্যাপ করেন। বুটস্ট্র্যাপ নমুনা থেকে, আপনি পরামিতিগুলি পেতে মডেলটিকে ফিট করে এবং আপনি প্রক্রিয়াটি পুনরাবৃত্তি করেন। প্যারামেট্রিক এবং ননপ্যারামেট্রিক বুটস্ট্র্যাপের মধ্যে কেবলমাত্র পার্থক্য হ'ল আপনি প্যারামেট্রিক বুটস্ট্র্যাপের অবশিষ্টাংশগুলি বুটস্ট্র্যাপ করুন যখন ননপ্যারমেট্রিক বুটস্ট্র্যাপ ভেক্টরগুলিকে বুটস্ট্র্যাপ করে। মিশ্র মডেলের ক্ষেত্রে আপনি কিছু প্রভাবগুলি প্যারামেট্রিকভাবে এবং অন্যগুলি অপরিকল্পিতভাবে চিকিত্সা করে একটি সেমিপ্রেমেট্রিক বুটস্ট্র্যাপ রাখতে পারেন। যদি আপনার কোডটি ভেক্টরগুলি বুটস্ট্র্যাপ করছে তবে আপনি ননপ্যারামেট্রিক বুটস্ট্র্যাপ করছেন। আর-তে এটি করার জন্য আমার কাছে কোনও নির্দিষ্ট সমাধান নেই তবে আপনি যদি এফ্রন এবং তিবশিরানীর বইটি বা রবার্ট লাবুদে সহ আমার বইটি দেখেন তবে আপনি রৈখিক মিশ্র মডেলের অনুরূপ মডেলের জন্য আর কোড দেখতে পাবেন।


আপনার সুন্দর উত্তরের জন্য @ মিশেল ধন্যবাদ। আপনি কার্যকরভাবে তিনটি বুটস্ট্র্যাপিং পদ্ধতির উদাহরণ ভাগ করে নিলে আমি অত্যন্ত প্রশংসা করব R
MYaseen208

আমার মনে হয় না যে আমার কাছে সেমিপ্রেমেট্রিক উদাহরণ রয়েছে। আমি আর আর প্রোগ্রামার নই (এখনও)। রবার্ট লাবুদে আমাদের বইতে সমস্ত তেহ আর প্রোগ্রামিং করেছিলেন। তিনি স্থির অটোরিগ্রেসিভ মডেল (পিপি 120-122) এর জন্য আর-তে কোডেড প্যারামেট্রিক বুটস্ট্র্যাপের একটি উদাহরণ দেন। বইয়ের 10 পৃষ্ঠায় তিনি আর এর সমস্ত বুটস্ট্র্যাপ ফাংশন দেখান যা "> সহায়তা.সন্ধান ('বুটস্ট্র্যাপ') ব্যবহার করে দেখায় ইফ্রন এবং তিবশিরানী তাদের বইতে আর-তে বুটস্ট্র্যাপ প্যাকেজটি নথিভুক্ত করেছে যে তারা" বুটস্ট্র্যাপ "বলে অন্যটিকে ডাকে ডেভিসন এবং হিঙ্কলির কারণে প্যাকেজটি হ'ল তাদের প্যাকেজ "বুট" যা আপনি তাদের বুটস্ট্র্যাপ বইয়ের বিষয়ে আরও জানতে পারেন
মাইকেল আর। চেরনিক

এটি উল্লেখ করা উচিত যে একটি মিশ্র মডেলের কঠোর ননপ্যারামেট্রিক বুটস্ট্র্যাপিং ব্যর্থতার ঝুঁকিপূর্ণ। এটি সম্পূর্ণরূপে একটি এলোমেলো ভেরিয়েবলের একটি স্তর বাদ দিতে পারে, যা প্রক্রিয়াটি থামিয়ে দেয়।
ব্রায়ান

2

আপনি bootMerlme4 এর বিকাশ সংস্করণে ফাংশনটি একবার দেখতে চান ,

install_github("lme4",user="lme4")
library(lme4)

যা মিশ্র মডেলগুলির মডেল-ভিত্তিক (আধা) প্যারামেট্রিক বুটস্ট্র্যাপিং করতে পারে ... কেবল পরীক্ষা করুন ?bootMer

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.