Last data update: 2014.03.03
R: Generating random missing values on a data matrix
missing.gen R Documentation
Generating random missing values on a data matrix
Description
The function generates a number of missing values (NA) completely at random on a data matrix. Totally missing rows (i.e., rows with all NA) are avoided.
Usage
missing.gen(mat, nummiss)
Arguments
mat
a matrix of numerical data.
nummiss
number of missing values.
Details
The function generates a number of missing values (NA) completely at random on a data matrix. Totally missing rows (i.e., rows with all NA) are avoided.
Value
The data matrix with missing values (NA).
Author(s)
Nadia Solaro, Alessandro Barbiero, Giancarlo Manzi, Pier Alda Ferrari
See Also
missing.gen0
Examples
sigma <- matrix(0.4, 4, 4)
diag(sigma) <- 1
x0 <- rmvnorm(n=100, mean=rep(0, 4), sigma=sigma)
x0 # complete matrix
x <- missing.gen(x0, 50)
x # matrix with 50 missing values
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(GenForImp)
Loading required package: mvtnorm
Loading required package: sn
Loading required package: stats4
Attaching package: 'sn'
The following object is masked from 'package:stats':
sd
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/GenForImp/missing.gen.Rd_%03d_medium.png", width=480, height=480)
> ### Name: missing.gen
> ### Title: Generating random missing values on a data matrix
> ### Aliases: missing.gen
> ### Keywords: classes multivariate NA
>
> ### ** Examples
>
> sigma <- matrix(0.4, 4, 4)
> diag(sigma) <- 1
> x0 <- rmvnorm(n=100, mean=rep(0, 4), sigma=sigma)
> x0 # complete matrix
[,1] [,2] [,3] [,4]
[1,] 0.65911650 -0.33331110 0.92160986 1.191047218
[2,] -0.97629778 0.37953656 -1.55386992 -1.127791828
[3,] -0.04720913 -0.58606968 0.74847725 0.609088748
[4,] 0.61129437 0.88537322 -0.74702612 -0.637236006
[5,] -0.12842654 0.10971713 -0.52177286 -0.067578832
[6,] 0.75896943 -0.63242417 -0.26258698 0.873436592
[7,] 0.23827133 0.37605678 -0.29748498 -0.260944278
[8,] -0.96512529 -0.76402807 -0.33611555 0.960044948
[9,] -0.09142860 1.41522037 -0.56131672 -0.015840755
[10,] -1.28898534 0.82731274 0.26209545 0.028187917
[11,] 1.08848018 1.08774018 0.44663401 2.277444058
[12,] 0.20239149 0.49657607 0.36386048 0.003071069
[13,] -0.30876106 0.10649953 0.01624501 1.827849608
[14,] 0.45718557 0.73528992 -1.02160523 0.958628557
[15,] -0.95542422 -1.80143390 -0.77911644 -2.562178742
[16,] -0.09625827 -1.08458374 0.31657756 -0.728423514
[17,] -0.72860965 -0.33460039 -0.22134101 -1.554741351
[18,] 0.13467581 -0.28493411 -1.03705698 -0.399347299
[19,] 0.52698943 0.64204252 -0.68263663 -0.392335844
[20,] -1.64894878 -0.84003164 -0.94553324 -0.699718484
[21,] 0.04177614 0.61652527 0.43348377 0.850892610
[22,] -0.65255709 0.59212214 0.31917875 0.157072780
[23,] -1.01898476 0.16907537 -0.31649785 0.940220646
[24,] 0.34221281 -1.78823411 -2.34172916 -0.986721444
[25,] 1.17356245 1.34681014 0.92968951 0.807312161
[26,] 0.56186163 1.49987668 0.37937205 1.609628155
[27,] -0.14419627 0.11237772 -1.08281431 -1.410903237
[28,] 0.43471086 -0.63245863 0.72686655 1.519624992
[29,] 1.31943059 1.40455267 -0.95041039 -0.084509771
[30,] 0.15071306 0.39806638 0.71314322 0.542321667
[31,] 0.06881901 0.68371377 0.41197381 0.345768337
[32,] -0.05729902 -0.08671733 -1.08994372 -0.195221008
[33,] -1.16459589 -0.26390110 -0.19328159 0.281006064
[34,] -1.36450546 -1.44313563 -0.48093300 -0.717994530
[35,] -0.44887600 -1.27143270 0.13370673 -0.013783312
[36,] -0.74449061 -1.33607331 -3.34093724 -2.113443164
[37,] -0.43355700 -1.13811810 0.19892824 0.204447922
[38,] -1.32169809 -0.62177593 -1.22447931 -0.548409911
[39,] 0.46781018 -0.77010747 0.69154779 0.375059818
[40,] 0.65513941 -1.06695438 0.48585316 0.152154514
[41,] 1.08381678 0.09107419 1.07343748 0.186962159
[42,] 0.40851667 1.62078381 1.02851592 0.055539876
[43,] -0.43297823 -1.75042359 -0.38471056 -0.464807194
[44,] -0.34326087 -0.43023342 -0.74181118 -0.304233785
[45,] 0.70029383 0.28042332 -0.87439959 -0.658390653
[46,] 0.29278737 -1.93864821 -1.95976927 -0.887221684
[47,] 0.36799756 0.66822889 -0.64932378 -1.443119858
[48,] 0.88233979 0.60531030 0.55191571 -0.107283182
[49,] 0.03872267 -1.32126762 0.36510610 -0.475244020
[50,] 1.10825886 0.93917666 1.02908110 0.125881526
[51,] 0.69410028 1.62879739 2.55159510 2.301717085
[52,] -0.69304801 1.50884895 1.67004555 -0.396016234
[53,] -0.04634516 1.86270973 1.06911015 0.954604692
[54,] 0.12902751 -0.48704483 -0.39441835 0.989061281
[55,] 0.74342926 0.36045239 0.02224567 0.082466535
[56,] -1.25118645 -1.39102069 -1.32213573 0.199272860
[57,] -1.65483510 0.26637252 0.73320502 0.180635410
[58,] 1.03036194 0.55486505 0.58671952 0.181288943
[59,] -0.92855784 -0.52050313 -1.24358883 -0.096320669
[60,] 0.60923559 0.01224115 -0.44697800 0.145259726
[61,] 0.16573811 0.80678068 0.06454526 0.859236467
[62,] 1.45748566 0.20033965 1.59803393 0.858646144
[63,] -0.13368987 -0.31916554 -2.07602146 -0.997027959
[64,] 0.70709789 -0.50037594 0.24291981 0.094340149
[65,] -0.05734386 -1.18871683 -0.02776583 1.040118010
[66,] -0.45683209 0.32238116 -0.43269009 -1.560045398
[67,] 0.19262867 0.42792662 0.20013635 0.263897327
[68,] -0.51917260 -0.20828194 -0.09649893 0.515952747
[69,] 0.76875167 0.97928099 -0.60842498 0.994259757
[70,] -0.57401844 -0.31979753 -0.39852976 -0.293839355
[71,] 0.44574651 -0.78286228 0.67732802 0.860872959
[72,] 0.01382214 -3.50347866 -0.89190932 0.314178947
[73,] -1.43331504 -1.23976839 -1.26546753 -0.157642437
[74,] 1.21654211 0.84048607 -0.24037903 -0.412349970
[75,] -0.71321677 -0.40589339 -1.06146454 -0.387413826
[76,] -1.06934530 -0.90245569 0.29049159 -0.674047650
[77,] -1.15371028 -0.74592942 0.29974108 -0.114228645
[78,] 2.59039103 -0.56363434 -0.15330687 1.372953562
[79,] -1.39598026 0.05312889 0.14913805 0.255675990
[80,] 1.47127532 0.41061791 1.32310047 1.546090232
[81,] 1.58650954 1.73629616 1.46094935 0.123005021
[82,] -0.34097341 -0.70810805 1.72370527 1.376248648
[83,] 0.01856457 0.49309434 2.03456715 -0.185176661
[84,] 1.74623010 1.70338311 -0.37877862 -0.549465477
[85,] 2.32607793 1.25787557 0.97900143 0.428970853
[86,] 1.48470239 1.43935892 1.71217005 1.220591167
[87,] -1.27177730 0.11399900 -1.33082494 -0.697988296
[88,] 0.67984825 0.99909587 -0.22604763 -0.280530267
[89,] 0.32480806 -0.22489570 -0.64776258 0.202174211
[90,] 0.18534331 0.92080565 -1.45553739 -0.902169327
[91,] 0.38858272 -0.74699128 -0.61418561 0.074640073
[92,] -0.15603877 -0.70586588 0.48774120 -2.101621283
[93,] 0.24500385 -0.56013223 0.20162786 -1.498879972
[94,] 0.89520242 0.73670128 1.28209018 -0.723480115
[95,] -0.65756120 -1.21044766 0.44562943 -0.194352604
[96,] 0.80358721 2.03053445 0.26257551 -0.146347153
[97,] -0.98633734 -0.64823098 -1.83090053 -1.064402321
[98,] 0.98505558 1.63204372 0.28637085 0.572921649
[99,] 0.49526330 -0.47829192 0.42367047 0.609106675
[100,] 0.93513553 -0.88636329 1.27365919 0.375093673
> x <- missing.gen(x0, 50)
> x # matrix with 50 missing values
[,1] [,2] [,3] [,4]
[1,] 0.65911650 -0.33331110 0.92160986 NA
[2,] -0.97629778 0.37953656 -1.55386992 -1.127791828
[3,] -0.04720913 -0.58606968 0.74847725 0.609088748
[4,] 0.61129437 0.88537322 NA -0.637236006
[5,] -0.12842654 0.10971713 -0.52177286 -0.067578832
[6,] 0.75896943 -0.63242417 NA 0.873436592
[7,] 0.23827133 0.37605678 NA -0.260944278
[8,] -0.96512529 -0.76402807 -0.33611555 0.960044948
[9,] -0.09142860 1.41522037 -0.56131672 -0.015840755
[10,] -1.28898534 0.82731274 0.26209545 0.028187917
[11,] 1.08848018 1.08774018 0.44663401 2.277444058
[12,] 0.20239149 0.49657607 0.36386048 0.003071069
[13,] -0.30876106 NA 0.01624501 NA
[14,] 0.45718557 0.73528992 -1.02160523 0.958628557
[15,] -0.95542422 NA -0.77911644 -2.562178742
[16,] -0.09625827 -1.08458374 0.31657756 -0.728423514
[17,] -0.72860965 -0.33460039 -0.22134101 -1.554741351
[18,] 0.13467581 -0.28493411 -1.03705698 -0.399347299
[19,] 0.52698943 0.64204252 -0.68263663 -0.392335844
[20,] -1.64894878 -0.84003164 -0.94553324 NA
[21,] 0.04177614 0.61652527 NA 0.850892610
[22,] -0.65255709 0.59212214 0.31917875 0.157072780
[23,] -1.01898476 0.16907537 -0.31649785 0.940220646
[24,] 0.34221281 -1.78823411 NA -0.986721444
[25,] 1.17356245 1.34681014 NA 0.807312161
[26,] 0.56186163 1.49987668 0.37937205 1.609628155
[27,] -0.14419627 0.11237772 -1.08281431 -1.410903237
[28,] 0.43471086 -0.63245863 0.72686655 1.519624992
[29,] 1.31943059 1.40455267 -0.95041039 -0.084509771
[30,] 0.15071306 NA 0.71314322 0.542321667
[31,] 0.06881901 0.68371377 0.41197381 0.345768337
[32,] NA -0.08671733 -1.08994372 -0.195221008
[33,] -1.16459589 NA -0.19328159 0.281006064
[34,] -1.36450546 -1.44313563 -0.48093300 -0.717994530
[35,] -0.44887600 -1.27143270 NA -0.013783312
[36,] -0.74449061 -1.33607331 -3.34093724 -2.113443164
[37,] -0.43355700 -1.13811810 0.19892824 0.204447922
[38,] -1.32169809 -0.62177593 -1.22447931 -0.548409911
[39,] 0.46781018 -0.77010747 0.69154779 0.375059818
[40,] 0.65513941 -1.06695438 0.48585316 0.152154514
[41,] 1.08381678 0.09107419 1.07343748 0.186962159
[42,] 0.40851667 1.62078381 1.02851592 0.055539876
[43,] -0.43297823 -1.75042359 -0.38471056 -0.464807194
[44,] -0.34326087 -0.43023342 -0.74181118 -0.304233785
[45,] NA 0.28042332 -0.87439959 -0.658390653
[46,] 0.29278737 -1.93864821 -1.95976927 -0.887221684
[47,] NA NA -0.64932378 -1.443119858
[48,] NA 0.60531030 0.55191571 -0.107283182
[49,] 0.03872267 -1.32126762 NA -0.475244020
[50,] 1.10825886 0.93917666 1.02908110 NA
[51,] 0.69410028 1.62879739 2.55159510 2.301717085
[52,] NA 1.50884895 1.67004555 NA
[53,] NA 1.86270973 NA 0.954604692
[54,] 0.12902751 -0.48704483 -0.39441835 0.989061281
[55,] 0.74342926 0.36045239 0.02224567 0.082466535
[56,] -1.25118645 -1.39102069 -1.32213573 0.199272860
[57,] -1.65483510 0.26637252 0.73320502 NA
[58,] 1.03036194 0.55486505 0.58671952 0.181288943
[59,] -0.92855784 -0.52050313 -1.24358883 -0.096320669
[60,] 0.60923559 0.01224115 -0.44697800 0.145259726
[61,] NA 0.80678068 0.06454526 0.859236467
[62,] 1.45748566 0.20033965 1.59803393 0.858646144
[63,] -0.13368987 NA -2.07602146 -0.997027959
[64,] 0.70709789 -0.50037594 0.24291981 NA
[65,] NA NA -0.02776583 1.040118010
[66,] -0.45683209 0.32238116 -0.43269009 -1.560045398
[67,] 0.19262867 0.42792662 NA 0.263897327
[68,] -0.51917260 -0.20828194 -0.09649893 0.515952747
[69,] 0.76875167 0.97928099 -0.60842498 0.994259757
[70,] -0.57401844 -0.31979753 -0.39852976 NA
[71,] 0.44574651 NA 0.67732802 0.860872959
[72,] 0.01382214 -3.50347866 NA 0.314178947
[73,] -1.43331504 -1.23976839 -1.26546753 NA
[74,] NA 0.84048607 -0.24037903 -0.412349970
[75,] -0.71321677 -0.40589339 -1.06146454 -0.387413826
[76,] -1.06934530 -0.90245569 0.29049159 -0.674047650
[77,] -1.15371028 -0.74592942 0.29974108 -0.114228645
[78,] 2.59039103 -0.56363434 -0.15330687 1.372953562
[79,] -1.39598026 0.05312889 0.14913805 0.255675990
[80,] 1.47127532 0.41061791 1.32310047 1.546090232
[81,] 1.58650954 1.73629616 1.46094935 0.123005021
[82,] -0.34097341 -0.70810805 1.72370527 1.376248648
[83,] 0.01856457 NA 2.03456715 NA
[84,] 1.74623010 NA NA -0.549465477
[85,] 2.32607793 1.25787557 0.97900143 0.428970853
[86,] NA 1.43935892 1.71217005 NA
[87,] -1.27177730 0.11399900 -1.33082494 -0.697988296
[88,] NA 0.99909587 -0.22604763 -0.280530267
[89,] 0.32480806 -0.22489570 -0.64776258 0.202174211
[90,] NA 0.92080565 -1.45553739 NA
[91,] 0.38858272 -0.74699128 -0.61418561 NA
[92,] -0.15603877 -0.70586588 0.48774120 -2.101621283
[93,] 0.24500385 -0.56013223 0.20162786 -1.498879972
[94,] 0.89520242 0.73670128 NA -0.723480115
[95,] -0.65756120 -1.21044766 0.44562943 -0.194352604
[96,] 0.80358721 2.03053445 0.26257551 -0.146347153
[97,] -0.98633734 -0.64823098 -1.83090053 -1.064402321
[98,] 0.98505558 1.63204372 0.28637085 0.572921649
[99,] 0.49526330 -0.47829192 NA 0.609106675
[100,] 0.93513553 -0.88636329 NA 0.375093673
>
>
>
>
>
> dev.off()
null device
1
>