|
0
|
1 #!/usr/bin/env python
|
|
|
2 """
|
|
|
3 Validate the metadata file associated with Affymetrix 96 well plate data.
|
|
|
4 """
|
|
|
5 import argparse
|
|
|
6 import datetime
|
|
|
7 import decimal
|
|
|
8 import re
|
|
|
9 import shutil
|
|
|
10 import sys
|
|
|
11
|
|
|
12 parser = argparse.ArgumentParser()
|
|
|
13 parser.add_argument('--input', dest='input', help='Metadata file for Affymetrix 96 well plate data')
|
|
|
14 parser.add_argument('--output', dest='output', help='Output dataset'),
|
|
|
15 args = parser.parse_args()
|
|
|
16
|
|
|
17 EMAIL_MAX_LEN = 255
|
|
|
18 VALID_EMAIL_RE = re.compile("[^@]+@[^@]+\.[^@]+")
|
|
|
19
|
|
|
20
|
|
|
21 def add_error_msg(accumulated_msgs, msg):
|
|
|
22 return "%s\n%s" % (accumulated_msgs, msg)
|
|
|
23
|
|
|
24
|
|
|
25 def empty_value(line_no, label, accumulated_msgs):
|
|
|
26 return add_error_msg(accumulated_msgs, "The required %s value is missing on line %d." % (label, line_no))
|
|
|
27
|
|
|
28
|
|
|
29 def stop_error(msg):
|
|
|
30 sys.exit(msg)
|
|
|
31
|
|
|
32
|
|
9
|
33 def string_as_boolean_string(string):
|
|
|
34 if str(string).lower() in ['true', 'yes', 'on', '1']:
|
|
|
35 return 'True'
|
|
|
36 else:
|
|
|
37 return 'False'
|
|
|
38
|
|
|
39
|
|
0
|
40 def validate_date_string(line_no, date_string, accumulated_msgs):
|
|
9
|
41 if len(date_string) == 0:
|
|
|
42 return accumulated_msgs
|
|
0
|
43 try:
|
|
9
|
44 datetime.datetime.strptime(date_string, '%Y-%m-%d')
|
|
0
|
45 return accumulated_msgs
|
|
|
46 except ValueError:
|
|
9
|
47 return add_error_msg(accumulated_msgs, "Line %d contains an incorrect date format (%s must be YYYY-MM-DD)." % (line_no, date_string))
|
|
0
|
48
|
|
|
49
|
|
|
50 def validate_decimal(line_no, decimal_string, accumulated_msgs):
|
|
|
51 try:
|
|
|
52 decimal.Decimal(decimal_string)
|
|
|
53 return accumulated_msgs
|
|
|
54 except Exception:
|
|
|
55 return add_error_msg(accumulated_msgs, "Line %d contains an incorrect decimal value (%s)." % (line_no, decimal_string))
|
|
|
56
|
|
|
57
|
|
|
58 def validate_email(line_no, email, accumulated_msgs):
|
|
|
59 if not (VALID_EMAIL_RE.match(email)):
|
|
|
60 return add_error_msg(accumulated_msgs, "Line %d contains an invalid email address (%s). " % (line_no, email))
|
|
|
61 elif len(email) > EMAIL_MAX_LEN:
|
|
|
62 return add_error_msg(accumulated_msgs, "Line %d contains an email address (%) that is longer than the maximum length, %d characters." % (line_no, email))
|
|
|
63 return accumulated_msgs
|
|
|
64
|
|
|
65
|
|
|
66 accumulated_msgs = ""
|
|
|
67 # Parse the input file, skipping the header, and validating
|
|
|
68 # that each data line consists of 31 comma-separated items.
|
|
|
69 with open(args.input, "r") as ih:
|
|
1
|
70 for i, line in enumerate(ih):
|
|
0
|
71 if i == 0:
|
|
|
72 # Skip the header.
|
|
|
73 continue
|
|
|
74 line = line.rstrip("\r\n")
|
|
|
75 if i > 97:
|
|
9
|
76 accumulated_msgs = add_error_msg(accumulated_msgs, "The input file contains more than 97 lines (must be 1 header line and no more than 96 data lines).")
|
|
0
|
77 stop_error(accumulated_msgs)
|
|
17
|
78 items = line.split("\t")
|
|
23
|
79 if len(items) != 32:
|
|
24
|
80 accumulated_msgs = add_error_msg(accumulated_msgs, "Line %d contains %s columns, (must be 32)." % (i, len(items)))
|
|
2
|
81 stop_error(accumulated_msgs)
|
|
0
|
82 # Required and validated.
|
|
|
83 # Required.
|
|
18
|
84 user_specimen_id = items[0]
|
|
0
|
85 if len(user_specimen_id) == 0:
|
|
|
86 accumulated_msgs = empty_value(i, "user_specimen_id", accumulated_msgs)
|
|
|
87 # Optional.
|
|
18
|
88 field_call = items[1]
|
|
0
|
89 # Optional.
|
|
18
|
90 bcoral_genet_id = items[2]
|
|
0
|
91 # Optional.
|
|
18
|
92 bsym_genet_id = items[3]
|
|
0
|
93 # Required.
|
|
18
|
94 reef = items[4]
|
|
0
|
95 if len(reef) == 0:
|
|
|
96 accumulated_msgs = empty_value(i, "reef", accumulated_msgs)
|
|
|
97 # Required.
|
|
18
|
98 region = items[5]
|
|
0
|
99 if len(region) == 0:
|
|
|
100 accumulated_msgs = empty_value(i, "region", accumulated_msgs)
|
|
|
101 # Required and validated.
|
|
18
|
102 latitude = items[6]
|
|
0
|
103 accumulated_msgs = validate_decimal(i, latitude, accumulated_msgs)
|
|
|
104 # Required and validated.
|
|
18
|
105 longitude = items[7]
|
|
0
|
106 accumulated_msgs = validate_decimal(i, longitude, accumulated_msgs)
|
|
|
107 # Optional.
|
|
18
|
108 geographic_origin = items[8]
|
|
9
|
109 # Optional.
|
|
18
|
110 sample_location = items[9]
|
|
0
|
111 # Optional.
|
|
18
|
112 latitude_outplant = items[10]
|
|
0
|
113 # Optional.
|
|
18
|
114 longitude_outplant = items[11]
|
|
0
|
115 # Optional.
|
|
18
|
116 depth = items[12]
|
|
0
|
117 # Optional.
|
|
24
|
118 disease_resist = items[13]
|
|
0
|
119 # Optional.
|
|
24
|
120 bleach_resist = items[14]
|
|
0
|
121 # Optional.
|
|
24
|
122 mortality = items[15]
|
|
0
|
123 # Optional.
|
|
24
|
124 tle = items[16]
|
|
0
|
125 # Optional.
|
|
24
|
126 spawning = string_as_boolean_string(items[17])
|
|
9
|
127 # Required.
|
|
24
|
128 collector_last_name = items[18]
|
|
9
|
129 if len(collector_last_name) == 0:
|
|
|
130 accumulated_msgs = empty_value(i, "collector_last_name", accumulated_msgs)
|
|
0
|
131 # Required.
|
|
24
|
132 collector_first_name = items[19]
|
|
9
|
133 if len(collector_first_name) == 0:
|
|
|
134 accumulated_msgs = empty_value(i, "collector_first_name", accumulated_msgs)
|
|
0
|
135 # Required.
|
|
24
|
136 org = items[20]
|
|
0
|
137 if len(org) == 0:
|
|
|
138 accumulated_msgs = empty_value(i, "org", accumulated_msgs)
|
|
|
139 # Required and validated.
|
|
24
|
140 collection_date = items[21]
|
|
18
|
141 accumulated_msgs = validate_date_string(i, collection_date, accumulated_msgs)
|
|
0
|
142 # Required and validated.
|
|
24
|
143 contact_email = items[22]
|
|
5
|
144 accumulated_msgs = validate_email(i, contact_email, accumulated_msgs)
|
|
0
|
145 # Required.
|
|
24
|
146 seq_facility = items[23]
|
|
0
|
147 if len(seq_facility) == 0:
|
|
|
148 accumulated_msgs = empty_value(i, "seq_facility", accumulated_msgs)
|
|
|
149 # Optional.
|
|
24
|
150 array_version = items[24]
|
|
0
|
151 # Optional.
|
|
24
|
152 public = string_as_boolean_string(items[25])
|
|
0
|
153 # Optional.
|
|
24
|
154 public_after_date = items[26]
|
|
13
|
155 accumulated_msga = validate_date_string(i, public_after_date, accumulated_msgs)
|
|
21
|
156 # Required and validated.
|
|
24
|
157 sperm_motility = items[27]
|
|
21
|
158 accumulated_msgs = validate_decimal(i, sperm_motility, accumulated_msgs)
|
|
|
159 # Required and validated.
|
|
24
|
160 healing_time = items[28]
|
|
21
|
161 accumulated_msgs = validate_decimal(i, healing_time, accumulated_msgs)
|
|
23
|
162 # Optional.
|
|
24
|
163 dna_extraction_method = items[29]
|
|
23
|
164 # Optional.
|
|
24
|
165 dna_concentration = items[30]
|
|
23
|
166 # If dna_concentration has a value, then it must be decimal.
|
|
|
167 if len(dna_concentration) > 0:
|
|
|
168 accumulated_msgs = validate_decimal(i, dna_concentration, accumulated_msgs)
|
|
24
|
169 # Optional.
|
|
|
170 registry_id = items[31]
|
|
23
|
171
|
|
0
|
172
|
|
6
|
173 if len(accumulated_msgs) > 0:
|
|
0
|
174 stop_error(accumulated_msgs)
|
|
|
175
|
|
8
|
176 shutil.copyfile(args.input, args.output)
|
