<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />目的
此腳本採用重複測量數據表(每個單元具有多行的表)
並將其轉換為每單元 1 行和 Pseudo Failure time 的數據表。
原始數據表必須包含:
1) 標識每個單元的列,
2) 測量時間的變量(可以使用,例如行駛里程),以及 3) a
響應列。
用法
只需通過以下任何一種方法運行此腳本:
編輯 > 運行腳本
控制-R
單擊工具欄中的“運行腳本”按鈕
//*******************一般信息
//此腳本採用重複測量數據表(每單元多行的表)
//並把它變成一個數據表,每個單元有1行和一個偽故障時間。
//原始數據表必須包含:
1) 標識每個單元的列,
2) 測量時間的變量(可以使用,例如行駛里程),以及
3) 一個響應列。
//**********************初始對話框
//初始對話框允許您選擇響應、時間度量
//(可以是時間或使用度量,例如行駛里程)和單位標識符。
//您也可以指定1個分組變量。 如果一個
//group by variable是指定的,用戶可以選擇不同的變換
//對於分組變量的每一層
//左下角的單選框允許用戶指定是否繪圖
//應該顯示回歸線或連接系列中的每個點
//單元。 這會影響顯示圖表,但不會改變內部計算的方式
//執行
//************************選擇轉換窗口
//指定初始變量後,腳本會打開一個窗口,顯示
//時間變量和響應之間的關係。 時間的轉變
//並且可以從側面的下拉列表中選擇響應
//地塊。 轉換時,繪圖將更新
//選擇顯示新的關係。
//如果指定了分組變量,每個變量都會有一個單獨的圖
//分組變量的級別。
//當用戶找到一個導致線性關係的轉換時
//轉換響應vs轉換時間,他/她可以選擇
//“創建偽故障時間數據”
//創建新的數據表。
//有關選擇轉換窗口的更多信息,請選擇幫助
//對話框中的按鈕。
//**********************創建數據表對話框
//此對話框允許您輸入失敗閾值,
//指定額外的協變量(可選),並在數據中指定一列
//控制偽故障時間大小的表。
//生成的數據表每個單元 ID 將有一行。
//細節
//失敗模型的協變量(可選):
// 要包含在失敗中的可選協變量列表
// 時間模型。 這些列從原始數據表複製到
// 偽故障時間表。 它們不影響計算
// 如果每個協變量有多個度量,則
// LAST VALUE 包含在最終表格中。 此時,腳本執行
// 不支持時間相關的協變量
//最大預測時間(可選)
// 包含最大預測時間的可選列。
// 如果線性模型預測的失效時間更大
// 比該列中的值,偽失敗時間為
// 此列的最大預測時間。
//審查列(可選)
// 數據表中的一列,對於測試期間失敗的觀察結果為 0
// 否則為 1。 如果指定了此列,則該列具有 1 的任何單位
// 將被視為“硬故障”,偽故障時間將是
// 對這個變量進行的最後一次測量。
//失敗閾值(必填)
// 單元被認為失敗的響應值。
//***************有關輸入、輸出和內部計算的詳細信息
//輸入:
// * 數據表,如假設中所述
// --將使用活動表或腳本將提示打開
// * 失敗級別:組件被認為“失敗”的損壞級別
// --用戶將輸入到一個帶標籤的框中
// * Censor Time (optional): 如果預測時間超過這個時間,組件
// 應被視為“審查”
// --數據表將有一列包含最大預測時間
// 用戶將允許每個單元
// * 時間和/或響應的轉換
// --用戶將從列表中選擇
// -- 將向用戶顯示轉換值和預測值的圖
// 並允許更改選擇
//輸出:
// * 交互式繪圖。 允許用戶在視覺上適合
// 線性變換
// * 包含以下內容的數據表:
// -- 每個組件的一行
// --“Predicted Time”:包含預測故障時間的數字列
// -- 一個包含原始失敗時間的數字列
// --"Pseudo Failure Time" 包含偽故障時間的數字列。
// 這是線性模型預測的時間,或者,
// 如果指定了審查變量,
// 適當的失敗時間。 這個變量是基於一個公式
// 分配原始的
// 時間或根據用戶指定的預測時間。
// --"Censored" 一個 0,1 變量,指示單位是否被刪失 (1)
// 審查 (0)。 這是應該指定為的變量
// 審查生命分配平台或其他生存中的變量
// 模型。
// --"Estimate Type":一個序數列,指示該行是否被審查
// 以下方式之一
// --測試失敗:0
//(在這種情況下,建模時間設置為最後一次
//單位被測量)。
// --在測試期間沒有失敗,並且預測的時間是有意義的:1
//(偽故障時間=預測時間)
// --predicted censoring time 比用戶指定的時間長:2
//(在這種情況下,建模時間設置為最大時間
// 由用戶指定)
// --該單元只有一行--無法預測:3
//(在這種情況下,建模時間是到最後一次測量單位)。
// --在某些情況下,預測時間的計算會產生一個時間
// 這對於 JSL 來說太大而無法處理。 在這種情況下,偽失敗設置為
// 最大時間或到最後一次測量單位 if
// 沒有最大時間列
// --來自原始數據表的協變量。 有一個對話框允許用戶
// 選擇將在以後的分析中使用的協變量。
// * 不支持與時間相關的協變量。 為每個測量的最後值
// 組件出現在數據表中。
//計算:
// * 輸入值的線性變換
// * 使用線性變換預測故障時間
// --模型為每個單元計算一個單獨的行
// -- 腳本使用 JSL 的矩陣語言執行計算
//假設:
// * 數據表格式如下:
// --1 行,每個項目的每次測量(重複測量)
// -- 包含一個時間列,指定採取每個措施的時間
// -- 包含一個審查列 (0,1) 指示組件是否
// 測試失敗。
// -- 如果單位的最後一個度量有 censoring=0,
// 新數據表將使用表中給出的失敗時間,而不是
// 預測失敗時間。
// * 有一個合適的線性變換適用於
// 所有行或按組內的所有行(即我們沒有情況
// 有些行需要一個
// 對數變換和其他需要冪變換)
//使用的平台:
// * 雙變量平台用於創建允許用戶選擇的圖
// 適當的轉換。 適合每個值用於允許用戶查看
// 如果擬合是近似線性的。
關於根據響應變量查找故障預測時間的方法的問題。
如果有故障前後的數據,如何計算故障時間。 假設該設備在 24 小時時很好,但在 96 小時時很差。
如果故障前只有數據,但設備處於故障軌跡上(它正在漂移)怎麼辦? 故障時間是如何計算的?
謝謝