উজ্জ্বল মধ্যে একটি রূপান্তর সতর্কতা অর্থ

আমি আর glmerএর lme4প্যাকেজ থেকে ফাংশনটি ব্যবহার করছি এবং আমি bobyqaঅপটিমাইজারটি ব্যবহার করছি (অর্থাত্ আমার ক্ষেত্রে এটি ডিফল্ট)। আমি একটি সতর্কতা পাচ্ছি, এবং আমি এর অর্থ কী তা জানতে আগ্রহী।

Warning message:
In optwrap(optimizer, devfun, start, rho$lower, control = control,  :
  convergence code 3 from bobyqa: bobyqa -- a trust region step failed to reduce q

আমি অনুসন্ধান করেছি "একটি বিশ্বস্ত অঞ্চল পদক্ষেপ Q হ্রাস করতে ব্যর্থ হয়েছে।" মিনকা প্যাকেজে কিছু তথ্য পাওয়া গেছে , যা "ব্যাখ্যার জন্য পাওলের সাথে পরামর্শ করুন।" আমি (আপনি চাইলে, আপনিও করতে পারেন! নীচে তাদের সাথে উল্লেখগুলি এবং লিঙ্কগুলি দেখুন), তবে আমি বুঝতে ব্যর্থ হয়েছি। আসলে, আমি q হ্রাস সম্পর্কে কিছুই খুঁজে পেতে ব্যর্থ।

এমজেডি পাওয়েল (২০০)) "ডেরিভেটিভস ছাড়াই নিয়ন্ত্রিত মিনিমাইজেশনের জন্য নিউইউইউএর উন্নয়ন", কেমব্রিজ বিশ্ববিদ্যালয়, ফলিত গণিত ও তাত্ত্বিক পদার্থ বিভাগ, সংখ্যা বিশ্লেষণ গ্রুপ, রিপোর্ট এনএ ২০০7 / ০ 05, http://www.damtp.cam.ac.uk/ ব্যবহারকারী / না / এনএ_প্যাপার্স / NA2007_05.pdf ।

এমজেডি পাওয়েল (২০০৯), "ডেরিভেটিভস ছাড়াই আবদ্ধ বাধিত অপ্টিমাইজেশনের জন্য বোবিকিউএ অ্যালগরিদম", রিপোর্ট নং ডিএএমটিপি ২০০৯ / এনএ 06, গণিত বিজ্ঞান কেন্দ্র, যুক্তরাজ্যের কেমব্রিজ বিশ্ববিদ্যালয়। http://www.damtp.cam.ac.uk/user/na/NA_papers/NA2009_06.pdf ।

পিএস আমি জানি আমি অপটিমাইজারটি পরিবর্তন করতে পারি এবং আমি সতর্কতা বা ত্রুটি ছাড়াই আউটপুট পেতে পারি কিনা তা আমি দেখতে চাই। আমিও, গ্রেডিয়েন্ট এবং চট চেক করবে যদি আমি করতে পারেন প্রতি একটি হিসাবে বেন Bolker মন্তব্য / উত্তর । আমি এর glmerমধ্যে dredgeথেকে ব্যবহার করছি MuMInএবং আমি নিশ্চিত নই যে বেনের উত্তর কোনও অতিরিক্ত টিঙ্কারিং ছাড়া কাজ করবে কিনা তবে আমার কম্পিউটার এটি কী করছে তা শেষ করার পরে আমি এটিতে কাজ করব, যাইহোক, আমি ডিগ্রি করি।

হালনাগাদ

ডঃ বলকারের নীচের মতামত অনুসারে, আমি ফরটান কোডটি সন্ধান করতে শুরু করেছি ( এটি দেখার জন্য তবে এটি ডাউনলোড না করতে আগ্রহী প্রত্যেকের জন্য কোড এখানে )। "430" কোডটির bobyqb.f অংশে উপস্থিত হয়। প্রাসঙ্গিক কোডটি সন্ধান করতে কেবল "430" বা "কিউ হ্রাস" অনুসন্ধান করুন।

এটি ফরট্রান কোডের সাথে আমার প্রথম মুখোমুখি, তবে আমি মনে করি কোডটি বলেছে যে নিম্নলিখিত শর্তগুলি পূরণ করা হলে সতর্কতা তৈরি করুন: NTRITS> 0, VQUAD> = 0, IPRINT> ০. "পূর্ণসংখ্যা এনটিআরআইটিএস" বিশ্বাস অঞ্চল "নাম্বারে সেট করা হয়েছে সর্বশেষ "বিকল্প" পুনরাবৃত্তির পরে ঘটে যাওয়া পুনরাবৃত্তিগুলি। VQUADবেশ কয়েকবার উপস্থিত হয়েছে এবং আমি এর তাত্পর্য সম্পর্কে এখনও স্পষ্ট নই কারণ এর মান বিভিন্ন ভেরিয়েবলের উপর নির্ভরশীল বলে মনে হয়, যার মানগুলি কখনও কখনও অন্যান্য ভেরিয়েবলের উপর নির্ভর করে F bobyqa.f থেকে: "আইপিআরএনটিটির মান হওয়া উচিত 0, 1, 2 বা 3 তে সেট করুন, যা মুদ্রণের পরিমাণ নিয়ন্ত্রণ করে Spec বিশেষত, আইপিআরএনটি = 0 হলে কোনও আউটপুট পাওয়া যায় না এবং কেবলমাত্র আইপিআরএনটি = 1 থাকলে রিটার্নে আউটপুট থাকে ""

সুতরাং, এটি টাস্কটি হ'ল VQUAD> = 0 হওয়ার তাৎপর্য বের করা এবং সম্ভবত কখন / কখন IPRINTপরিণত হয়েছিল তা বোঝা 0. 0 আমাকে একবার কাগজে ফিরে যেতে হবে তবে গণিত, বা এ এটির প্রতীকী প্রকাশটি আমার পক্ষে কিছুটা বাধা। অ্যালগরিদম সম্পর্কে কেউ না জেনে বা এটি সম্পর্কে জানতে আগ্রহী না হওয়া পর্যন্ত আমি মনে করি যে আমি কাগজপত্র, কোড এবং ইন্টারনেটের মধ্যে আগে থেকে পিছনে গিয়ে সতর্কবার্তাটি বোঝার চেষ্টা করব যতক্ষণ না আমি এটি কী বুঝতে পারি উপায়।

— Jota
সূত্র

আমি মনে করি এই প্রশ্নটি সিভি বি / সি-র জন্য বিষয় হতে পারে বলে মনে হচ্ছে ধারণাটি বুঝতে পারার পরিবর্তে প্রতি ডাব্লু / আর-র সাহায্য না করে।

— গুং - মনিকা পুনরায়

আমি নিশ্চিত নই যে কাগজপত্র এবং ফরট্রান কোডের (যা cran.r-project.org/src/contrib/minqa_1src ডিরেক্টরিতে অন্তর্ভুক্ত রয়েছে) এর মাধ্যমে বিট-বিট বিস্তৃত হওয়ার বাইরে আমার এই ক্ষেত্রে অনেক পরামর্শ দেওয়ার দরকার আছে। 2.3.tar.gz ) এবং এই ত্রুটিটি ( কোডটিতে ত্রুটি কোড 430) ট্রিগার হয়ে গেলে কী ঘটছে তা অবিকল দেখছেন ...

— বেন বোলকার

দ্রুত কাগজটির উপরে স্কিমিং করা আমি মনে করি সতর্কতাটি ইঙ্গিত দেয় যে অপ্টিমাইজারটি আপনি যে ফাংশনটি কমিয়ে আনতে চান, তার ক্রমটি চতুর্ভুজ, Q, হ্রাস হ্রাস করতে পারে, এটি হ'ল, অপটিমাইজার এমন এক পর্যায়ে রয়েছে যা সম্ভবত অনুকূল নয় তবে উদ্দেশ্যটির উন্নতি করতে কোন উপায়টি তা জানেন না। সুতরাং, এটি আটকে আছে।

— সেভেন

দুটি পেপারগুলির মধ্যে আপনি কোনটি এড়িয়ে গেছেন এবং আনুমানিক কোথায় আপনি এই তথ্যটি পেয়েছেন? (আমি খুব

— ঝাঁকুনিতে পড়েছি

আমি BOBYQA কাগজ পড়েছি। তারা কী চলছে এবং কি কি হচ্ছে এর একটি বিস্তৃত ধারণা পেতে আমি প্রায় অর্ধেকের মধ্যে প্রথমার্ধে গিয়েছিলাম। সত্যই কোনও নির্দিষ্ট পৃষ্ঠায় নির্দেশ করতে পারে না।

— সেভেন

$f(x)$ $x_k$ $k$

$\Delta_k$
$f(x)$ $x = x_k$ $Q(x)$
$s_k$ $Q_k(x_k + s_k)$ $||s_k|| \leq \Delta_k$
$s_k$ $x_{k+1} = x_k + s_k$
অন্যথায়, আপনার মডেলটিকে পরিমার্জন করুন এবং আবার চেষ্টা করুন

$s_k$

এর মান বুঝতে VQUAD, আমাদের প্রথমে কয়েকটি অন্যান্য ভেরিয়েবল বুঝতে হবে। ভাগ্যক্রমে, ঘোষণার ঠিক নীচে ভাল মন্তব্য রয়েছে SUBROUTINE BOBYQB। প্রধান ভেরিয়েবলগুলি হ'ল:

GOPT, মডেলের গ্রেডিয়েন্ট
HQ, মডেল হেসিয়ান
D $s_k$

410 এর উপরে কয়েকটি লাইন শুরু করা, আপনি দেখতে পাবেন DO 410 J=1,N। এটি একটি ফর-লুপ (এবং নেস্টেড ফর লুপ) শুরু করে যা পরীক্ষামূলক পদক্ষেপ ব্যবহার করে মডেল দ্বারা পূর্বাভাসিত পরিবর্তনের মূল্যায়ন করে D। এটি মধ্যে পূর্বাভাস পরিবর্তন জমে VQUAD। ফর-লুপের প্রথম অংশটি প্রথম-ক্রমের শর্তাদি মূল্যায়ন করে এবং নেস্টেড ফর-লুপ দ্বিতীয়-ক্রমের শর্তাদি মূল্যায়ন করে। লুপগুলি ইন্ডেন্ট করা থাকলে সম্ভবত এটি পড়া সহজ হবে:

    DO 410 J=1,N
        VQUAD=VQUAD+D(J)*GOPT(J)
        DO 410 I=1,J
            IH=IH+1
            TEMP=D(I)*D(J)
            IF (I .EQ. J) TEMP=HALF*TEMP
410         VQUAD=VQUAD+HQ(IH)*TEMP

মডেলের অন্যান্য প্যারামিটারগুলি অন্তর্ভুক্ত করার জন্য এর পরে আরও একটি লুপ রয়েছে। আমাকে স্বীকার করতে হবে, আমি এটি পুরোপুরি বুঝতে পারি না - আমার সর্বোত্তম অনুমান যে তারা মডেলটি কীভাবে তৈরি করেন তা এটি নির্দিষ্ট।

এই সব শেষে, VQUADমডেল দ্বারা পূর্বাভাস উদ্দেশ্যমূলক ফাংশন পরিবর্তন ধারণ করে। সুতরাং যদি VQUADঅ-নেতিবাচক হয় তবে এটি খারাপ। এখন এই বিশেষ সমাধানকারী একটি বিকল্প পদক্ষেপের গণনা (সম্ভবত একটি লাইন অনুসন্ধান) ব্যবহার করতে পারেন, এটিই NTRITSখেলতে আসে। সুতরাং 430-এর যুক্তিটি বলছে, "যদি সর্বশেষ পুনরাবৃত্তিটি বিকল্প পদক্ষেপের গণনা ব্যবহার করে এবং মডেল হ'ল AND IPRINT> 0 হ্রাসের পূর্বাভাস না দেয় তবে সতর্কতা বার্তাটি মুদ্রণ করুন" " নোট করুন যে সলভার এর মান নির্বিশেষে সমাপ্ত হতে চলেছে IPRINT।

কথা বলছি IPRINT, সেই মানটি BOBYQAকলিং ফাংশনটি দিয়ে যায়। এই ক্ষেত্রে, আপনার আর রুটিন হল কলিং ফাংশন। এখানে একটি verboseপ্যারামিটার রয়েছে glmer- আমি ডলারে ডাইম হয়ে যাব যে একই মানটি পাস করা হবে BOBYQA। verbose0 এ সেট করার চেষ্টা করুন এবং আপনি সম্ভবত সতর্কতাটি দেখতে পাবেন না। তবে অবশ্যই হুডের নীচে কী চলছে তা পরিবর্তন হবে না।

— বিল ওউসনার
সূত্র

এটি খুব সহায়ক; আমি মনে করি না আমি আরও ভাল করতে যাচ্ছি, অনুগ্রহ প্রদান করে ...

— বেন বলকার

@BenBolker। সুতরাং আমি এখানে কী ঘটছে সে সম্পর্কে উদ্বিগ্ন হওয়া উচিত , বা এটি কি কোডের মধ্যে কেবল উপদ্রব? (অন্য কথায়, এই সতর্কতার অর্থ কি আমার ফলাফলগুলি বৈধ নয় ??)

— বনাঞ্চল বিশেষজ্ঞ

এক্ষেত্রে সাধারণ পরামর্শটি যা বলা হয়েছে ?lme4::convergence: অপ্টিমাইজেশান পদ্ধতির সম্পূর্ণ / বিশদ মূল্যায়নের সংক্ষিপ্ততা, আপনার সেরা সেরাটি বিভিন্ন অপ্টিমাইজারের থেকে ফলাফলের তুলনা করা।

— বেন বলকার