[R] Reliability metric

[David Winsemius]

On Sep 9, 2011, at 3:14 PM, Halabi, Anan wrote:

>
> Hi,
> Below is the code I write,
> I am trying to create a metric of h and t while the values is out of  
> function R.
> First I have message warning

It is only a warning.

> Second, the metric is not created

What "not created" means has already been question by Weylandt. What  
you seem to have missed in your introductory material is the concept  
of "recycling of arguments". Look at this code and its results:

 > (1:10) * (1:2)
# [1]  1  4  3  8  5 12  7 16  9 20

 > (1:10) * (1:3)
# [1]  1  4  9  4 10 18  7 16 27 10
Warning message:
In (1:10) * (1:3) :
   longer object length is not a multiple of shorter object length

-- 
David.
>
>
>> h <- seq(0.1, 0.9, by=0.1)
>> t <- seq(0,11000, by=100)
>> z <- cbind(t)
>> eta=10000
>> beta=2
>> R <- array (1:1100, dim= c(110,10))
>> R= exp(-(z/eta*(1-h))^(beta*(1-h)))
> Warning messages:
> 1: In z/eta * (1 - h) :
>  longer object length is not a multiple of shorter object length
> 2: In (z/eta * (1 - h))^(beta * (1 - h)) :
>  longer object length is not a multiple of shorter object length
>> R
>               t
>  [1,] 1.0000000
>  [2,] 0.9995586
>  [3,] 0.9974647
>  [4,] 0.9919725
>  [5,] 0.9801987
>  [6,] 0.9572081
>  [7,] 0.9141353
>  [8,] 0.8341665
>  [9,] 0.6833618
> [10,] 0.9892125
> [11,] 0.9825766
> [12,] 0.9727664
> [13,] 0.9583520
> [14,] 0.9370675
> [15,] 0.9051390
> [16,] 0.8559267
> [17,] 0.7769472
> [18,] 0.6423128
> [19,] 0.9629361
> [20,] 0.9521001
> [21,] 0.9382259
> [22,] 0.9201104
> [23,] 0.8958341
> [24,] 0.8622051
> [25,] 0.8136275
> [26,] 0.7395534
> [27,] 0.6175830
> [28,] 0.9246320
> [29,] 0.9127580
> [30,] 0.8982795
> [31,] 0.8800821
> [32,] 0.8564152
> [33,] 0.8244048
> [34,] 0.7790512
> [35,] 0.7109198
> [36,] 0.5996181
> [37,] 0.8767646
> [38,] 0.8671159
> [39,] 0.8550329
> [40,] 0.8394499
> [41,] 0.8187308
> [42,] 0.7902253
> [43,] 0.7493487
> [44,] 0.6874269
> [45,] 0.5854263
> [46,] 0.8215803
> [47,] 0.8170930
> [48,] 0.8098560
> [49,] 0.7988930
> [50,] 0.7827045
> [51,] 0.7588540
> [52,] 0.7231046
> [53,] 0.6673758
> [54,] 0.5736597
> [55,] 0.7611968
> [56,] 0.7642511
> [57,] 0.7637415
> [58,] 0.7588504
> [59,] 0.7482636
> [60,] 0.7297788
> [61,] 0.6994875
> [62,] 0.6498159
> [63,] 0.5635906
> [64,] 0.6975913
> [65,] 0.7098972
> [66,] 0.7174446
> [67,] 0.7196221
> [68,] 0.7153381
> [69,] 0.7026455
> [70,] 0.6779572
> [71,] 0.6341555
> [72,] 0.5547793
> [73,] 0.6325678
> [74,] 0.6551296
> [75,] 0.6715521
> [76,] 0.6814181
> [77,] 0.6838614
> [78,] 0.6771931
> [79,] 0.6581383
> [80,] 0.6199984
> [81,] 0.5469392
> [82,] 0.5677232
> [83,] 0.6008650
> [84,] 0.6265229
> [85,] 0.6443861
> [86,] 0.6537698
> [87,] 0.6532207
> [88,] 0.6397566
> [89,] 0.6070646
> [90,] 0.5398727
> [91,] 0.5044222
> [92,] 0.5478565
> [93,] 0.5827146
> [94,] 0.6086282
> [95,] 0.6250023
> [96,] 0.6305684
> [97,] 0.6226043
> [98,] 0.5951486
> [99,] 0.5334377
> [100,] 0.4437834
> [101,] 0.4967086
> [102,] 0.5404019
> [103,] 0.5742126
> [104,] 0.5975006
> [105,] 0.6091056
> [106,] 0.6065193
> [107,] 0.5840942
> [108,] 0.5275280
> [109,] 0.3866770
> [110,] 0.4478912
> [111,] 0.4997913
>> length(R)
> [1] 111
>
> Anan Halabi
> Reliability Eng, R&D
> HP Scitex
> Tel: 972-9-8924648
> mobil: 972-52-6624231
>
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius, MD
West Hartford, CT