inst/test_objs/wqmmlesss/superalloy.R
In Auburngrads/SMRD: Statistical Methods For Reliability Data
col.of.ones <-
FALSE
data.ld <-
structure(list(kilocycles = c(5.733, 13.949, 15.616, 56.723,
12.076, 152.68, 43.331, 18.067, 9.75, 156.725, 112.968, 138.114,
122.372, 21.3, 6.705, 112.002, 11.865, 13.181, 8.489, 12.434,
13.03, 57.923, 121.075, 200.027, 211.629, 155), event = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L), .Label = c("Failure", "Right"
), class = "factor"), count = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L), pstress = c(145.9, 85.2, 116.4, 87.2, 100.1, 85.8, 99.8,
113, 120.4, 86.4, 85.6, 86.7, 89.7, 114.8, 144.5, 91.3, 142.5,
100.5, 118.4, 118.6, 118, 80.8, 87.3, 80.6, 80.3, 84.3), lpstress = c(4.98292,
4.445, 4.75703, 4.4682, 4.60617, 4.45202, 4.60317, 4.72739, 4.79082,
4.45899, 4.44969, 4.46245, 4.49647, 4.74319, 4.97328, 4.51415,
4.95934, 4.61016, 4.77407, 4.77576, 4.77068, 4.39198, 4.46935,
4.3895, 4.38577, 4.43438), lpstress2 = c(24.82951, 19.75804,
22.62936, 19.96485, 21.2168, 19.82047, 21.18916, 22.3482, 22.95195,
19.88257, 19.7997, 19.91349, 20.21825, 22.49787, 24.73351, 20.37756,
24.59507, 21.25355, 22.79173, 22.80785, 22.75943, 19.28946, 19.97509,
19.2677, 19.23498, 19.66374)), class = c("regrs3.life.data",
"life.data", "data.frame"), row.names = c(NA, -26L), right.censor.names = "a,alive,c,censor,censored,end,mend,noreport,r,r-censored,right-censored,removed,right,rightcensored,s,survived,survive,suspend,suspended,2", left.censor.names = "l,l-censored,left-censored,left,leftcensored,start,mstart,3", interval.censor.names = "b,bin,i,interval,i-censored,intervalcensored,interval-censored,4", failure.censor.names = "event,exact,d,dead,died,f,fail,failed,failure,report,repair,repaired,replaced,replacement,1", sinterval.censor.names = "s,sinterval,smallinterval,small-interval,5", response.column = c(kilocycles = "kilocycles"), censor.column = c(event = "event"), case.weight.column = c(count = "count"), data.title = "superalloy", time.units = c(kilocycles = "kilocycles"), x.columns = c(pstress = "pstress",
lpstress = "lpstress", lpstress2 = "lpstress2"), data.note = "")
debug1 <-
FALSE
distribution <-
"lognormal"
distribution.number <-
4
dummy <-
c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1,
1, 2, 1, 1, 1, 1)
e <-
c(1e-04, 1e-04, 1e-04, 1e-04, 1e-04)
escale <-
10000
explan.vars <-
1:3
gamthr <-
c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0)
int <-
1
intercept <-
TRUE
intercept.increment <-
1
is.eyring <-
FALSE
ivec <-
c(number.case = 26, nter = 4, nyresp = 1, ntyresp = 0, distribution.number = 4,
lcheck = 0, nparm = 5, int = 1, maxit = 500, kprint = 0, ierfit = 0,
iervcv = 0)
kprint <-
0
likelihood.method <-
"smoothed"
mathsoft.gamthr <-
c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0)
maxit <-
500
ndscrat <-
316
niscrat <-
12
non.pos.resp <-
c(FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE,
FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE,
FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE)
nparm <-
5
nter <-
4L
ntyresp <-
0
number.cases <-
26L
nyresp <-
1L
oldClasses <-
"regrs3.life.data"
param.names <-
c("b0", "pstress", "lpstress", "lpstress2", "sigma")
parameter.fixed <-
c(FALSE, FALSE, FALSE, FALSE, FALSE)
regression <-
TRUE
RegrNames <-
c("pstress", "lpstress", "lpstress2")
relationship <-
NULL
rvec <-
c(0, escale = 10000, log.likelihood = 0)
startna <-
c(TRUE, TRUE, TRUE, TRUE, TRUE)
test <-
7
the.case.weights <-
c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)
the.censor.codes <-
c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1,
1, 2, 1, 1, 1, 1)
the.truncation.codes <-
c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1)
the.xmat <-
structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 145.9, 85.2, 116.4, 87.2, 100.1, 85.8,
99.8, 113, 120.4, 86.4, 85.6, 86.7, 89.7, 114.8, 144.5, 91.3,
142.5, 100.5, 118.4, 118.6, 118, 80.8, 87.3, 80.6, 80.3, 84.3,
4.98292, 4.445, 4.75703, 4.4682, 4.60617, 4.45202, 4.60317, 4.72739,
4.79082, 4.45899, 4.44969, 4.46245, 4.49647, 4.74319, 4.97328,
4.51415, 4.95934, 4.61016, 4.77407, 4.77576, 4.77068, 4.39198,
4.46935, 4.3895, 4.38577, 4.43438, 24.82951, 19.75804, 22.62936,
19.96485, 21.2168, 19.82047, 21.18916, 22.3482, 22.95195, 19.88257,
19.7997, 19.91349, 20.21825, 22.49787, 24.73351, 20.37756, 24.59507,
21.25355, 22.79173, 22.80785, 22.75943, 19.28946, 19.97509, 19.2677,
19.23498, 19.66374), .Dim = c(26L, 4L), .Dimnames = list(NULL,
c("", "pstress", "lpstress", "lpstress2")))
theta.start <-
c(162.477484560072, 0.0749744343916901, -57.7097548399362, 4.67877353472178,
0.306010225415608)
theta.start.comp <-
c(162.477484560072, 0.0749744343916901, -57.7097548399362, 4.67877353472178,
0.306010225415608)
tyresp <-
c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0)
yresp <-
structure(c(5.733, 13.949, 15.616, 56.723, 12.076, 152.68, 43.331,
18.067, 9.75, 156.725, 112.968, 138.114, 122.372, 21.3, 6.705,
112.002, 11.865, 13.181, 8.489, 12.434, 13.03, 57.923, 121.075,
200.027, 211.629, 155), .Dim = c(26L, 1L), .Dimnames = list(NULL,
c(kilocycles = "kilocycles")))
zout <-
list(ivec = c(26L, 4L, 1L, 0L, 4L, 0L, 5L, 1L, 500L, 0L, 0L,
0L), rvec = structure(c(0, 10000, -98.574462890625), Csingle = TRUE),
number.cases = 26L, nparm = 5L, xmat = structure(c(1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 145.899993896484, 85.1999969482422, 116.400001525879,
87.1999969482422, 100.099998474121, 85.8000030517578, 99.8000030517578,
113, 120.400001525879, 86.4000015258789, 85.5999984741211,
86.6999969482422, 89.6999969482422, 114.800003051758, 144.5,
91.3000030517578, 142.5, 100.5, 118.400001525879, 118.599998474121,
118, 80.8000030517578, 87.3000030517578, 80.5999984741211,
80.3000030517578, 84.3000030517578, 4.98292016983032, 4.44500017166138,
4.75703001022339, 4.46820020675659, 4.60617017745972, 4.45202016830444,
4.6031699180603, 4.72738981246948, 4.79082012176514, 4.4589900970459,
4.4496898651123, 4.46245002746582, 4.49646997451782, 4.74318981170654,
4.97327995300293, 4.51415014266968, 4.95934009552002, 4.6101598739624,
4.77406978607178, 4.77576017379761, 4.77067995071411, 4.39198017120361,
4.46934986114502, 4.3895001411438, 4.38576984405518, 4.43438005447388,
24.8295097351074, 19.7580394744873, 22.6293601989746, 19.9648494720459,
21.2168006896973, 19.8204708099365, 21.1891593933105, 22.3481998443604,
22.9519500732422, 19.8825702667236, 19.7996997833252, 19.9134902954102,
20.2182502746582, 22.4978694915771, 24.7335090637207, 20.3775596618652,
24.5950698852539, 21.2535495758057, 22.7917308807373, 22.8078498840332,
22.7594299316406, 19.2894592285156, 19.9750900268555, 19.2677001953125,
19.2349796295166, 19.6637401580811), Csingle = TRUE), yresp = structure(c(1.74623894691467,
2.6354079246521, 2.74829602241516, 4.03817987442017, 2.49121999740601,
5.02834415435791, 3.7688684463501, 2.89408707618713, 2.27726721763611,
5.05449295043945, 4.72710466384888, 4.92807960510254, 4.80706548690796,
3.05870699882507, 1.90285348892212, 4.71851682662964, 2.47359275817871,
2.57877635955811, 2.13877129554749, 2.52043461799622, 2.56725430488586,
4.05911445617676, 4.79641008377075, 5.29845237731934, 5.35483455657959,
5.0434250831604), Csingle = TRUE), structure(c(1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), Csingle = TRUE),
tyresp = structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), Csingle = TRUE),
truncation.codes = structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), Csingle = TRUE),
parameter.fixed = c(FALSE, FALSE, FALSE, FALSE, FALSE), e = structure(c(1.0000000116861e-07,
1.0000000116861e-07, 1.0000000116861e-07, 1.0000000116861e-07,
1.0000000116861e-07), Csingle = TRUE), ndscrat = c(-0.272048591185772,
-0.160918402913318, -0.218139328049096, -0.164594989148301,
-0.188278051766237, -0.162021536622068, -0.187727847007101,
-0.211915673045318, -0.225458190308356, -0.163124531191021,
-0.161653835999595, -0.163675972170028, -0.16918879354705,
-0.215210861916616, -0.269493360117429, -0.17212775749162,
-0.2658425269835, -0.189011632141424, -0.221799164300334,
-0.222165098151209, -0.221067259118756, -0.152824791290892,
-0.164778792048146, -0.152456717267067, -0.151904595999867,
-0.159263483963994, -0.392078134596684, 0.197129761233392,
-0.0853341733347657, 0.180779715775771, 0.0689894617748955,
0.192254751994403, 0.07170298299045, -0.0521043494270362,
-0.125058630638091, 0.187354165732469, 0.193882931149882,
0.184893769819287, 0.159947224428667, -0.0696334866073754,
-0.376873563413248, 0.146392751332287, -0.355252314444521,
0.0653632071241083, -0.105113486609582, -0.107100418352312,
-0.101144364827746, 0.23204832779536, 0.179954480114474,
0.233599983447248, 0.235920978416669, 0.204391698851601,
0.205142733778425, 0.0425047602457583, -0.168853586952413,
0.141339698194682, 0.231276882404968, 0.0793173128050823,
0.237605068314846, -0.0846216384175995, -0.240436031265336,
0.111163233129827, 0.0712602780866516, 0.119846370862436,
0.228276741377881, -0.132081144298253, 0.129508300306438,
0.260613440702354, 0.0285752038754299, 0.226475627638526,
-0.205626619150437, -0.2069814644753, -0.203451452577207,
-0.30836494152614, 0.149388938439295, -0.331806721985138,
-0.366109996115812, -0.0143750611585846, -0.367316925422332,
-0.0929509627726615, 0.222199154889537, -0.0347338013607677,
0.205051042857077, -0.0748839997530497, 0.201926629795687,
0.241732685984165, 0.184940648975126, -0.0573411622579664,
-0.0808754623324668, -0.0487278019941064, 0.0295963665009092,
0.232822206157822, -0.324407427718035, 0.0660564845040099,
-0.265390290610709, 0.209001757802188, 0.205422816110325,
0.203522775013726, 0.209120821740065, -0.243167944321851,
-0.0320225493959938, -0.250717816493719, -0.262192490558188,
-0.121157373117266, 544.028414574992, 12.6245090143564, 0.000464850124155653,
0.099332612193532, 0, 544.028414574992, 12.6245090143564,
0.000464850124155653, 0.099332612193532, 1, -0.00920299684650748,
-0.978997004837596, -0.0427701611757984, -0.199125293925851,
0, 0.0764317462669101, -0.20318165028925, 0.278418094589266,
0.935606097539776, 0, -0.890755402878415, 0.00110767429062365,
0.450367722460029, -0.0610122927510824, 0, -0.447915511290406,
-0.0167587160134649, -0.847244352288314, 0.28507516149007,
0, 0, 0, 0, 0, 1, 159.480357547547, 0.0789829211588798, -56.1714662537231,
4.46722251849938, 0.622712597158995, -17.4609198590718, 9.01039443185972,
-0.0779220237313291, -2.24199644022407, 0.622712597158995,
1.03368896264848e-08, 6.1922222628893e-08, 9.1711808979689e-08,
-2.39137486425147e-07, -8.21100930174995e-07, 1602517.92484304,
-1982.55443140329, -809560.802446477, 109569.397154723, -3.6291089130661,
-1982.55443140329, 2.46524626139566, 1002.32239483082, -135.781674408635,
0.00474394034947916, -809560.802446477, 1002.32239483082,
409022.654099572, -55366.5085171776, 1.85381329469482, 109569.397154723,
-135.781674408635, -55366.5085171776, 7495.78551935397, -0.253737683335568,
-3.6291089130661, 0.00474394034947916, 1.85381329469482,
-0.253737683335568, 0.00900737224479613, 0.410552725195885,
-0.0297180041670799, 0.0127193277264155, 0.0203430641886371,
-0.00431313363323937, -0.0297180041670799, 0.427510774461553,
-0.0228698483819016, -0.0296461430934869, 0.0049693637852819,
0.012719327724465, -0.0228698483820367, 0.436276184808221,
0.0315794473085141, 0.00189799812546457, 0.0203430640067381,
-0.0296461430993986, 0.0315794473085041, 0.415898181828869,
-0.0017395607027195, -0.00431313363323937, 0.00496936378528057,
0.00189799812546457, -0.00173956070271947, 0.00900737224479613,
-16.8559924900869, 8.8938099743442, -0.0778719945042013,
-2.23662067447407, -0.473670187407928, -16.8559924900869,
8.8938099743442, -0.0778719945042013, -2.23662067447407,
-0.473670187407928, 1.23660847332748e-07, 1.23660847332748e-07,
1.23660847332748e-07, 1.23660847332748e-07, 1.0000000116861e-07,
0.0178581557445247, 0.0260819065953998, 0.0154953735439986,
0.0104511197209923, 0.0150021623436904, 0.0692461467476347,
0, 0, 0, 0, 0, 0.047412502947372, 0, 0, 0, 0, 0, 0.07980501210998,
0, 0, 0, 0, 0, 0.118323060718901, 0, 0, 0, 0, 0, 0.0666570584144278,
2.06101375965773e-08, 0, 0, -4.44089209850063e-16, -5.55111512312578e-17,
0, 0, 0, 0, 0, 0, -17.4609198590711, 9.01039443185966, -0.0779220237313291,
-2.24199644022408, 0.622712597158995, 0.00991094275622827,
-0.0112305966533537, -0.00558704764000531, 0.00339356622070928,
0.999866448841777, 0.627473987970219, 0.692392344119565,
-0.0659669863001409, 0.350025131873861, 7.26304827743343e-07,
-0.483781496790864, 0.00422168524415608, -0.383680106630351,
0.78659215331873, 2.91550638402485e-05, -0.473604703544061,
0.465049247472005, 0.744592416724462, 0.0694138761836834,
0.0138430110975356, -0.384489552235385, 0.551523426794226,
-0.542207386845377, -0.503910100013438, 0.00868646497384417,
43.6075339960987, -8.23720678094108e-18, 7.46652276440504e-15,
5.55117381615358e-17, -2.3732086588876e-18, -1.61558713389263e-27,
0.999999948240337, -3.45434235224535e-18, -4.29693667256604e-18,
-1.22222302987982e-13, 7.447897005093e-15, -3.05366768811573e-21,
0.99986420331803, -6.3229577456098e-18, -1.66924811346982e-14,
5.86930017960607e-22, 2.58489384977489e-26, 2.58493941422821e-26,
0.930170986860423, 1.093614288447e-13, -1.10071143586102e-24,
-1.22195023862954e-13, -1.67694589214338e-14, 1.09360796667858e-13,
0.778848937572703, 0.00496936378528195), niscrat = c(1L,
2L, 3L, 4L, 5L, 0L, 1L, 2L, 3L, 4L, 5L, 0L), theta.hat = structure(c(159.480361938477,
0.0789829194545746, -56.1714668273926, 4.46722269058228,
0.6227126121521), Csingle = TRUE), first.derivative = structure(c(1.03368895665312e-08,
6.19222220166193e-08, 9.17118114784898e-08, -2.39137477819895e-07,
-8.21100911707617e-07), Csingle = TRUE), vcv.matrix = structure(c(1602517.875,
-1982.55444335938, -809560.8125, 109569.3984375, -3.62910890579224,
-1982.55444335938, 2.46524620056152, 1002.32238769531, -135.781677246094,
0.00474394019693136, -809560.8125, 1002.32238769531, 409022.65625,
-55366.5078125, 1.85381329059601, 109569.3984375, -135.781677246094,
-55366.5078125, 7495.78564453125, -0.253737688064575, -3.62910890579224,
0.00474394019693136, 1.85381329059601, -0.253737688064575,
0.00900737196207047), Csingle = TRUE), correlation.matrix = structure(c(1,
-0.997455477714539, -0.999940931797028, 0.999722003936768,
-0.0302064381539822, -0.997455477714539, 1, 0.998167932033539,
-0.99885493516922, 0.0318353772163391, -0.999940931797028,
0.998167932033539, 1, -0.999919176101685, 0.0305417031049728,
0.999722003936768, -0.99885493516922, -0.999919176101685,
1, -0.0308800023049116, -0.0302064381539822, 0.0318353772163391,
0.0305417031049728, -0.0308800023049116, 1), Csingle = TRUE),
residuals = structure(c(0.63914692401886, 0.0313747748732567,
1.36396968364716, 0.425873070955276, 0.220956474542618, 1.63276100158691,
1.66140580177307, 1.37586140632629, 0.802886426448822, 1.89532685279846,
0.97111576795578, 1.63007056713104, 2.21718859672546, 2.02721238136292,
0.819184124469757, 2.46718263626099, 2.02681088447571, 0.266159176826477,
0.577037036418915, 1.07738256454468, 1.12178349494934, 0.130222976207733,
1.46457695960999, 0.913447737693787, 0.938291907310486, 1.26856291294098
), Csingle = TRUE), fitted.values.and.deviance = structure(c(7.5759449005127,
120.431838989258, 12.8713827133179, 96.5186309814453, 30.9194355010986,
112.510360717773, 31.5868091583252, 14.8113222122192, 11.1783494949341,
105.249732971191, 115.048767089844, 101.881225585938, 74.5326156616211,
13.7173290252686, 7.59166765213013, 63.825927734375, 7.64207410812378,
30.0552272796631, 11.9552326202393, 11.870080947876, 12.1301279067993,
206.132354736328, 95.4692153930664, 211.627151489258, 220.191513061523,
133.658935546875, -2.29168963432312, -9.07253646850586, -3.24173831939697,
-4.84777641296387, -4.07622051239014, -5.59379625320435,
-4.34299802780151, -3.39026141166687, -2.74663496017456,
-5.70417213439941, -0.670034050941467, -1.16297733783722,
-1.54671537876129, -3.75366067886353, -2.3680112361908, -5.57155895233154,
-3.16840648651123, -3.90008354187012, -2.73520660400391,
-2.96848058700562, -3.01912593841553, -0.020967923104763,
-5.31447505950928, -5.74781847000122, -5.80213117599487,
-5.51698780059814, -0.447620898485184, -3.46175098419189,
0.310399293899536, -0.853613972663879, -1.5097895860672,
0.490272462368011, 0.507664144039154, 0.319080054759979,
-0.219541996717453, 0.639391303062439, 0.779319643974304,
1.13276600837708, 1.3644859790802, 0.706661641597748, -0.199446365237236,
0.903076887130737, 0.706463515758514, -1.32366073131561,
-0.549848794937134, 0.0745345130562782, 0.114919818937778,
0.0510103814303875, 0.381566405296326, -0.0905291438102722,
-0.0636941567063332, 0.237884685397148, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0.630180597305298, 0.729663014411926, 0.775363624095917,
1, 1, 1, 1, 1, 1, 1, 1, 0.106587082147598, 1, 1, 1, 1), Csingle = TRUE))
Auburngrads/SMRD documentation built on Sept. 14, 2020, 2:21 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
col.of.ones <-
FALSE
data.ld <-
structure(list(kilocycles = c(5.733, 13.949, 15.616, 56.723,
12.076, 152.68, 43.331, 18.067, 9.75, 156.725, 112.968, 138.114,
122.372, 21.3, 6.705, 112.002, 11.865, 13.181, 8.489, 12.434,
13.03, 57.923, 121.075, 200.027, 211.629, 155), event = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L), .Label = c("Failure", "Right"
), class = "factor"), count = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L), pstress = c(145.9, 85.2, 116.4, 87.2, 100.1, 85.8, 99.8,
113, 120.4, 86.4, 85.6, 86.7, 89.7, 114.8, 144.5, 91.3, 142.5,
100.5, 118.4, 118.6, 118, 80.8, 87.3, 80.6, 80.3, 84.3), lpstress = c(4.98292,
4.445, 4.75703, 4.4682, 4.60617, 4.45202, 4.60317, 4.72739, 4.79082,
4.45899, 4.44969, 4.46245, 4.49647, 4.74319, 4.97328, 4.51415,
4.95934, 4.61016, 4.77407, 4.77576, 4.77068, 4.39198, 4.46935,
4.3895, 4.38577, 4.43438), lpstress2 = c(24.82951, 19.75804,
22.62936, 19.96485, 21.2168, 19.82047, 21.18916, 22.3482, 22.95195,
19.88257, 19.7997, 19.91349, 20.21825, 22.49787, 24.73351, 20.37756,
24.59507, 21.25355, 22.79173, 22.80785, 22.75943, 19.28946, 19.97509,
19.2677, 19.23498, 19.66374)), class = c("regrs3.life.data",
"life.data", "data.frame"), row.names = c(NA, -26L), right.censor.names = "a,alive,c,censor,censored,end,mend,noreport,r,r-censored,right-censored,removed,right,rightcensored,s,survived,survive,suspend,suspended,2", left.censor.names = "l,l-censored,left-censored,left,leftcensored,start,mstart,3", interval.censor.names = "b,bin,i,interval,i-censored,intervalcensored,interval-censored,4", failure.censor.names = "event,exact,d,dead,died,f,fail,failed,failure,report,repair,repaired,replaced,replacement,1", sinterval.censor.names = "s,sinterval,smallinterval,small-interval,5", response.column = c(kilocycles = "kilocycles"), censor.column = c(event = "event"), case.weight.column = c(count = "count"), data.title = "superalloy", time.units = c(kilocycles = "kilocycles"), x.columns = c(pstress = "pstress",
lpstress = "lpstress", lpstress2 = "lpstress2"), data.note = "")
debug1 <-
FALSE
distribution <-
"lognormal"
distribution.number <-
4
dummy <-
c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1,
1, 2, 1, 1, 1, 1)
e <-
c(1e-04, 1e-04, 1e-04, 1e-04, 1e-04)
escale <-
10000
explan.vars <-
1:3
gamthr <-
c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0)
int <-
1
intercept <-
TRUE
intercept.increment <-
1
is.eyring <-
FALSE
ivec <-
c(number.case = 26, nter = 4, nyresp = 1, ntyresp = 0, distribution.number = 4,
lcheck = 0, nparm = 5, int = 1, maxit = 500, kprint = 0, ierfit = 0,
iervcv = 0)
kprint <-
0
likelihood.method <-
"smoothed"
mathsoft.gamthr <-
c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0)
maxit <-
500
ndscrat <-
316
niscrat <-
12
non.pos.resp <-
c(FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE,
FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE,
FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE)
nparm <-
5
nter <-
4L
ntyresp <-
0
number.cases <-
26L
nyresp <-
1L
oldClasses <-
"regrs3.life.data"
param.names <-
c("b0", "pstress", "lpstress", "lpstress2", "sigma")
parameter.fixed <-
c(FALSE, FALSE, FALSE, FALSE, FALSE)
regression <-
TRUE
RegrNames <-
c("pstress", "lpstress", "lpstress2")
relationship <-
NULL
rvec <-
c(0, escale = 10000, log.likelihood = 0)
startna <-
c(TRUE, TRUE, TRUE, TRUE, TRUE)
test <-
7
the.case.weights <-
c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)
the.censor.codes <-
c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1,
1, 2, 1, 1, 1, 1)
the.truncation.codes <-
c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1)
the.xmat <-
structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 145.9, 85.2, 116.4, 87.2, 100.1, 85.8,
99.8, 113, 120.4, 86.4, 85.6, 86.7, 89.7, 114.8, 144.5, 91.3,
142.5, 100.5, 118.4, 118.6, 118, 80.8, 87.3, 80.6, 80.3, 84.3,
4.98292, 4.445, 4.75703, 4.4682, 4.60617, 4.45202, 4.60317, 4.72739,
4.79082, 4.45899, 4.44969, 4.46245, 4.49647, 4.74319, 4.97328,
4.51415, 4.95934, 4.61016, 4.77407, 4.77576, 4.77068, 4.39198,
4.46935, 4.3895, 4.38577, 4.43438, 24.82951, 19.75804, 22.62936,
19.96485, 21.2168, 19.82047, 21.18916, 22.3482, 22.95195, 19.88257,
19.7997, 19.91349, 20.21825, 22.49787, 24.73351, 20.37756, 24.59507,
21.25355, 22.79173, 22.80785, 22.75943, 19.28946, 19.97509, 19.2677,
19.23498, 19.66374), .Dim = c(26L, 4L), .Dimnames = list(NULL,
c("", "pstress", "lpstress", "lpstress2")))
theta.start <-
c(162.477484560072, 0.0749744343916901, -57.7097548399362, 4.67877353472178,
0.306010225415608)
theta.start.comp <-
c(162.477484560072, 0.0749744343916901, -57.7097548399362, 4.67877353472178,
0.306010225415608)
tyresp <-
c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0)
yresp <-
structure(c(5.733, 13.949, 15.616, 56.723, 12.076, 152.68, 43.331,
18.067, 9.75, 156.725, 112.968, 138.114, 122.372, 21.3, 6.705,
112.002, 11.865, 13.181, 8.489, 12.434, 13.03, 57.923, 121.075,
200.027, 211.629, 155), .Dim = c(26L, 1L), .Dimnames = list(NULL,
c(kilocycles = "kilocycles")))
zout <-
list(ivec = c(26L, 4L, 1L, 0L, 4L, 0L, 5L, 1L, 500L, 0L, 0L,
0L), rvec = structure(c(0, 10000, -98.574462890625), Csingle = TRUE),
number.cases = 26L, nparm = 5L, xmat = structure(c(1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 145.899993896484, 85.1999969482422, 116.400001525879,
87.1999969482422, 100.099998474121, 85.8000030517578, 99.8000030517578,
113, 120.400001525879, 86.4000015258789, 85.5999984741211,
86.6999969482422, 89.6999969482422, 114.800003051758, 144.5,
91.3000030517578, 142.5, 100.5, 118.400001525879, 118.599998474121,
118, 80.8000030517578, 87.3000030517578, 80.5999984741211,
80.3000030517578, 84.3000030517578, 4.98292016983032, 4.44500017166138,
4.75703001022339, 4.46820020675659, 4.60617017745972, 4.45202016830444,
4.6031699180603, 4.72738981246948, 4.79082012176514, 4.4589900970459,
4.4496898651123, 4.46245002746582, 4.49646997451782, 4.74318981170654,
4.97327995300293, 4.51415014266968, 4.95934009552002, 4.6101598739624,
4.77406978607178, 4.77576017379761, 4.77067995071411, 4.39198017120361,
4.46934986114502, 4.3895001411438, 4.38576984405518, 4.43438005447388,
24.8295097351074, 19.7580394744873, 22.6293601989746, 19.9648494720459,
21.2168006896973, 19.8204708099365, 21.1891593933105, 22.3481998443604,
22.9519500732422, 19.8825702667236, 19.7996997833252, 19.9134902954102,
20.2182502746582, 22.4978694915771, 24.7335090637207, 20.3775596618652,
24.5950698852539, 21.2535495758057, 22.7917308807373, 22.8078498840332,
22.7594299316406, 19.2894592285156, 19.9750900268555, 19.2677001953125,
19.2349796295166, 19.6637401580811), Csingle = TRUE), yresp = structure(c(1.74623894691467,
2.6354079246521, 2.74829602241516, 4.03817987442017, 2.49121999740601,
5.02834415435791, 3.7688684463501, 2.89408707618713, 2.27726721763611,
5.05449295043945, 4.72710466384888, 4.92807960510254, 4.80706548690796,
3.05870699882507, 1.90285348892212, 4.71851682662964, 2.47359275817871,
2.57877635955811, 2.13877129554749, 2.52043461799622, 2.56725430488586,
4.05911445617676, 4.79641008377075, 5.29845237731934, 5.35483455657959,
5.0434250831604), Csingle = TRUE), structure(c(1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), Csingle = TRUE),
tyresp = structure(c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), Csingle = TRUE),
truncation.codes = structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), Csingle = TRUE),
parameter.fixed = c(FALSE, FALSE, FALSE, FALSE, FALSE), e = structure(c(1.0000000116861e-07,
1.0000000116861e-07, 1.0000000116861e-07, 1.0000000116861e-07,
1.0000000116861e-07), Csingle = TRUE), ndscrat = c(-0.272048591185772,
-0.160918402913318, -0.218139328049096, -0.164594989148301,
-0.188278051766237, -0.162021536622068, -0.187727847007101,
-0.211915673045318, -0.225458190308356, -0.163124531191021,
-0.161653835999595, -0.163675972170028, -0.16918879354705,
-0.215210861916616, -0.269493360117429, -0.17212775749162,
-0.2658425269835, -0.189011632141424, -0.221799164300334,
-0.222165098151209, -0.221067259118756, -0.152824791290892,
-0.164778792048146, -0.152456717267067, -0.151904595999867,
-0.159263483963994, -0.392078134596684, 0.197129761233392,
-0.0853341733347657, 0.180779715775771, 0.0689894617748955,
0.192254751994403, 0.07170298299045, -0.0521043494270362,
-0.125058630638091, 0.187354165732469, 0.193882931149882,
0.184893769819287, 0.159947224428667, -0.0696334866073754,
-0.376873563413248, 0.146392751332287, -0.355252314444521,
0.0653632071241083, -0.105113486609582, -0.107100418352312,
-0.101144364827746, 0.23204832779536, 0.179954480114474,
0.233599983447248, 0.235920978416669, 0.204391698851601,
0.205142733778425, 0.0425047602457583, -0.168853586952413,
0.141339698194682, 0.231276882404968, 0.0793173128050823,
0.237605068314846, -0.0846216384175995, -0.240436031265336,
0.111163233129827, 0.0712602780866516, 0.119846370862436,
0.228276741377881, -0.132081144298253, 0.129508300306438,
0.260613440702354, 0.0285752038754299, 0.226475627638526,
-0.205626619150437, -0.2069814644753, -0.203451452577207,
-0.30836494152614, 0.149388938439295, -0.331806721985138,
-0.366109996115812, -0.0143750611585846, -0.367316925422332,
-0.0929509627726615, 0.222199154889537, -0.0347338013607677,
0.205051042857077, -0.0748839997530497, 0.201926629795687,
0.241732685984165, 0.184940648975126, -0.0573411622579664,
-0.0808754623324668, -0.0487278019941064, 0.0295963665009092,
0.232822206157822, -0.324407427718035, 0.0660564845040099,
-0.265390290610709, 0.209001757802188, 0.205422816110325,
0.203522775013726, 0.209120821740065, -0.243167944321851,
-0.0320225493959938, -0.250717816493719, -0.262192490558188,
-0.121157373117266, 544.028414574992, 12.6245090143564, 0.000464850124155653,
0.099332612193532, 0, 544.028414574992, 12.6245090143564,
0.000464850124155653, 0.099332612193532, 1, -0.00920299684650748,
-0.978997004837596, -0.0427701611757984, -0.199125293925851,
0, 0.0764317462669101, -0.20318165028925, 0.278418094589266,
0.935606097539776, 0, -0.890755402878415, 0.00110767429062365,
0.450367722460029, -0.0610122927510824, 0, -0.447915511290406,
-0.0167587160134649, -0.847244352288314, 0.28507516149007,
0, 0, 0, 0, 0, 1, 159.480357547547, 0.0789829211588798, -56.1714662537231,
4.46722251849938, 0.622712597158995, -17.4609198590718, 9.01039443185972,
-0.0779220237313291, -2.24199644022407, 0.622712597158995,
1.03368896264848e-08, 6.1922222628893e-08, 9.1711808979689e-08,
-2.39137486425147e-07, -8.21100930174995e-07, 1602517.92484304,
-1982.55443140329, -809560.802446477, 109569.397154723, -3.6291089130661,
-1982.55443140329, 2.46524626139566, 1002.32239483082, -135.781674408635,
0.00474394034947916, -809560.802446477, 1002.32239483082,
409022.654099572, -55366.5085171776, 1.85381329469482, 109569.397154723,
-135.781674408635, -55366.5085171776, 7495.78551935397, -0.253737683335568,
-3.6291089130661, 0.00474394034947916, 1.85381329469482,
-0.253737683335568, 0.00900737224479613, 0.410552725195885,
-0.0297180041670799, 0.0127193277264155, 0.0203430641886371,
-0.00431313363323937, -0.0297180041670799, 0.427510774461553,
-0.0228698483819016, -0.0296461430934869, 0.0049693637852819,
0.012719327724465, -0.0228698483820367, 0.436276184808221,
0.0315794473085141, 0.00189799812546457, 0.0203430640067381,
-0.0296461430993986, 0.0315794473085041, 0.415898181828869,
-0.0017395607027195, -0.00431313363323937, 0.00496936378528057,
0.00189799812546457, -0.00173956070271947, 0.00900737224479613,
-16.8559924900869, 8.8938099743442, -0.0778719945042013,
-2.23662067447407, -0.473670187407928, -16.8559924900869,
8.8938099743442, -0.0778719945042013, -2.23662067447407,
-0.473670187407928, 1.23660847332748e-07, 1.23660847332748e-07,
1.23660847332748e-07, 1.23660847332748e-07, 1.0000000116861e-07,
0.0178581557445247, 0.0260819065953998, 0.0154953735439986,
0.0104511197209923, 0.0150021623436904, 0.0692461467476347,
0, 0, 0, 0, 0, 0.047412502947372, 0, 0, 0, 0, 0, 0.07980501210998,
0, 0, 0, 0, 0, 0.118323060718901, 0, 0, 0, 0, 0, 0.0666570584144278,
2.06101375965773e-08, 0, 0, -4.44089209850063e-16, -5.55111512312578e-17,
0, 0, 0, 0, 0, 0, -17.4609198590711, 9.01039443185966, -0.0779220237313291,
-2.24199644022408, 0.622712597158995, 0.00991094275622827,
-0.0112305966533537, -0.00558704764000531, 0.00339356622070928,
0.999866448841777, 0.627473987970219, 0.692392344119565,
-0.0659669863001409, 0.350025131873861, 7.26304827743343e-07,
-0.483781496790864, 0.00422168524415608, -0.383680106630351,
0.78659215331873, 2.91550638402485e-05, -0.473604703544061,
0.465049247472005, 0.744592416724462, 0.0694138761836834,
0.0138430110975356, -0.384489552235385, 0.551523426794226,
-0.542207386845377, -0.503910100013438, 0.00868646497384417,
43.6075339960987, -8.23720678094108e-18, 7.46652276440504e-15,
5.55117381615358e-17, -2.3732086588876e-18, -1.61558713389263e-27,
0.999999948240337, -3.45434235224535e-18, -4.29693667256604e-18,
-1.22222302987982e-13, 7.447897005093e-15, -3.05366768811573e-21,
0.99986420331803, -6.3229577456098e-18, -1.66924811346982e-14,
5.86930017960607e-22, 2.58489384977489e-26, 2.58493941422821e-26,
0.930170986860423, 1.093614288447e-13, -1.10071143586102e-24,
-1.22195023862954e-13, -1.67694589214338e-14, 1.09360796667858e-13,
0.778848937572703, 0.00496936378528195), niscrat = c(1L,
2L, 3L, 4L, 5L, 0L, 1L, 2L, 3L, 4L, 5L, 0L), theta.hat = structure(c(159.480361938477,
0.0789829194545746, -56.1714668273926, 4.46722269058228,
0.6227126121521), Csingle = TRUE), first.derivative = structure(c(1.03368895665312e-08,
6.19222220166193e-08, 9.17118114784898e-08, -2.39137477819895e-07,
-8.21100911707617e-07), Csingle = TRUE), vcv.matrix = structure(c(1602517.875,
-1982.55444335938, -809560.8125, 109569.3984375, -3.62910890579224,
-1982.55444335938, 2.46524620056152, 1002.32238769531, -135.781677246094,
0.00474394019693136, -809560.8125, 1002.32238769531, 409022.65625,
-55366.5078125, 1.85381329059601, 109569.3984375, -135.781677246094,
-55366.5078125, 7495.78564453125, -0.253737688064575, -3.62910890579224,
0.00474394019693136, 1.85381329059601, -0.253737688064575,
0.00900737196207047), Csingle = TRUE), correlation.matrix = structure(c(1,
-0.997455477714539, -0.999940931797028, 0.999722003936768,
-0.0302064381539822, -0.997455477714539, 1, 0.998167932033539,
-0.99885493516922, 0.0318353772163391, -0.999940931797028,
0.998167932033539, 1, -0.999919176101685, 0.0305417031049728,
0.999722003936768, -0.99885493516922, -0.999919176101685,
1, -0.0308800023049116, -0.0302064381539822, 0.0318353772163391,
0.0305417031049728, -0.0308800023049116, 1), Csingle = TRUE),
residuals = structure(c(0.63914692401886, 0.0313747748732567,
1.36396968364716, 0.425873070955276, 0.220956474542618, 1.63276100158691,
1.66140580177307, 1.37586140632629, 0.802886426448822, 1.89532685279846,
0.97111576795578, 1.63007056713104, 2.21718859672546, 2.02721238136292,
0.819184124469757, 2.46718263626099, 2.02681088447571, 0.266159176826477,
0.577037036418915, 1.07738256454468, 1.12178349494934, 0.130222976207733,
1.46457695960999, 0.913447737693787, 0.938291907310486, 1.26856291294098
), Csingle = TRUE), fitted.values.and.deviance = structure(c(7.5759449005127,
120.431838989258, 12.8713827133179, 96.5186309814453, 30.9194355010986,
112.510360717773, 31.5868091583252, 14.8113222122192, 11.1783494949341,
105.249732971191, 115.048767089844, 101.881225585938, 74.5326156616211,
13.7173290252686, 7.59166765213013, 63.825927734375, 7.64207410812378,
30.0552272796631, 11.9552326202393, 11.870080947876, 12.1301279067993,
206.132354736328, 95.4692153930664, 211.627151489258, 220.191513061523,
133.658935546875, -2.29168963432312, -9.07253646850586, -3.24173831939697,
-4.84777641296387, -4.07622051239014, -5.59379625320435,
-4.34299802780151, -3.39026141166687, -2.74663496017456,
-5.70417213439941, -0.670034050941467, -1.16297733783722,
-1.54671537876129, -3.75366067886353, -2.3680112361908, -5.57155895233154,
-3.16840648651123, -3.90008354187012, -2.73520660400391,
-2.96848058700562, -3.01912593841553, -0.020967923104763,
-5.31447505950928, -5.74781847000122, -5.80213117599487,
-5.51698780059814, -0.447620898485184, -3.46175098419189,
0.310399293899536, -0.853613972663879, -1.5097895860672,
0.490272462368011, 0.507664144039154, 0.319080054759979,
-0.219541996717453, 0.639391303062439, 0.779319643974304,
1.13276600837708, 1.3644859790802, 0.706661641597748, -0.199446365237236,
0.903076887130737, 0.706463515758514, -1.32366073131561,
-0.549848794937134, 0.0745345130562782, 0.114919818937778,
0.0510103814303875, 0.381566405296326, -0.0905291438102722,
-0.0636941567063332, 0.237884685397148, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0.630180597305298, 0.729663014411926, 0.775363624095917,
1, 1, 1, 1, 1, 1, 1, 1, 0.106587082147598, 1, 1, 1, 1), Csingle = TRUE))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.