global変数について教えてください

Question

with open("sequence.fasta.txt") as f:
    
    max_len = 0
    max_protein = ""
    
    min_len = 0
    min_protein = ""
    
    def change(protein, seq):
        global max_len,max_protein,min_len,min_protein
        seq_len = len(seq)
        if seq_len == 0:
            return

        else:
            a = [seq_len]
            print(a)
            
        if seq_len > max_len:
            max_len,max_protein,min_len,min_protein = seq_len,protein,seq_len,protein
            
        else:
            if seq_len < min_len:
                min_len,min_protein = seq_len,protein

ファイルの中から最大値と最小値を検索し続け、見つかり次第、最小値、最大値、その時のタンパク質を更新していくのですが、change関数内のglobalをやめて、関数外で同じ機能を果たせるようにしたいです。
強いて言えば、globalを使用せずに他の関数で現在のコードのままの機能を果たすにはどのような関数を定義したら最適でしょうか？

また、globalは関数外から定義した数値を変更できる。との認識ですが、間違っていますでしょうか？
pythonは通常がlocalなのでできるだけglobalを使用したくないと考えています。

無理な要望ではあるかと思いますが、知っている方教えていただきたいです。
宜しくお願い致します。

Comments

[KojiDoi]

「関数外で同じ機能を果たせる」とは具体的にどういうことでしょうか。最終的に今のプログラムでできていない何をやりたいのかが全く見えません。

[22Go]

今のプログラムの機能で満足しているのですが、できるだけglobalを使わずにこの機能を持つコードを書きたいと思ってます。

Answer 1

returnで、戻り値で呼び出し側に返せば良いでしょう。
#追記
return命令の使い方：

def foo():
    x = 123
    return x

kekka = foo()
print(kekka)

Comments

[22Go]

どのように現コードに組み込んだらいいでしょうか？

[hayataka2049]

提示されているコードにchange関数の呼び出し部分が書かれていないので、具体的な書き換え方法は誰も答えようがないと思います。

[22Go]

globalを消してchange関数を呼び出したらいいのでしょうか？

[hayataka2049]

あなたが質問に示しているコードではchangeは定義されただけで、どう使っているのかがわかりません。たぶん何らかのループで呼び出しているのだろうなぁと想像できる程度です。
なので、具体的な助言が得たいのなら、changeを呼び出して処理を行っている部分のコードを追記した方が良いです。

[otn]

return max_len,max_protein,min_len,min_protein

[22Go]

with open("sequence.fasta.txt") as f:

max_len = 0
max_protein = ""

min_len = 1000000000000
min_protein = ""

def change(protein, seq):
global max_len,max_protein,min_len,min_protein
seq_len = len(seq)
if seq_len == 0:
return

else:
print([seq_len])

if seq_len > max_len:
max_len,max_protein,min_len,min_protein = seq_len,protein,seq_len,protein

if seq_len < min_len:
min_len,min_protein = seq_len,protein

a_line = ""
a_seq = ""

for line in f:
strip_line = line.rstrip()
if strip_line.startswith(">"):
a_line = strip_line
change(a_line, a_seq)
a_seq = ""
print(strip_line+'\n')
else:
#A = [len(a_seq+strip_line)]
a_seq += strip_line
msg1 = "最大アミノ酸エントリ："
msg2 = "長さ："
msg3 = "最小アミノ酸エントリ："
print(msg1,max_protein,msg2,max_len)
print(msg3,min_protein,msg2,min_len)

このように下の方でもchange関数を使っているのですがこれで関数の使いたあってますでしょうか？

[22Go]

return max_len,max_protein,min_len,min_protein
はどの部分に組み込んだらいいでしょうか？

[otn]

return命令の意味を復習しましょう。

[22Go]

関数内で戻り値を使用することはわかったのですが、
global変数を使用しない場合、max_len,max_protein,min_len,min_protein はchange関数に組み込んでいいのでしょうか？

[otn]

はい。関数の中のローカル変数にします。

[22Go]

with open("short.sequ.txt") as f:

max_len = 0
max_protein = ""

min_len = 10000000000000
min_protein = ""

def change(protein,seq,max_len,max_protein,min_len,min_protein):
seq_len = len(seq)
if seq_len == 0:
return

else:
print([seq_len])

if seq_len > max_len:
seq_len = max_len
max_protein = protein
return protein,seq,max_len,max_protein,min_len,min_protein

if seq_len < min_len:
seq_len = min_len
min_protein = protein
return protein,seq,max_len,max_protein,min_len,min_protein

a_line = ""
a_seq = ""

for line in f:
strip_line = line.rstrip()
if strip_line.startswith(">"):
a_line = strip_line
change(a_line, a_seq,max_len,max_protein,min_len,min_protein)
a_seq = ""
print(strip_line+'\n')
else:
A = [len(a_seq+strip_line)]
a_seq += strip_line

msg1 = "最大アミノ酸エントリ："
msg2 = "長さ："
msg3 = "最小アミノ酸エントリ："
print(msg1,max_protein,msg2,max_len)
print(msg3,min_protein,msg2,min_len)

このように組み込んで見たのですが、
最後の部分が
最大アミノ酸エントリ： 長さ： 0
最小アミノ酸エントリ： 長さ： 10000000000000
となってしまったのですが、どのように修正したらいいでしょうか？

[otn]

return命令の使い方を追記しました。

[22Go]

ありがとうございます。
ですが、どうしてもうまく機能しませんでした。
理解できずに申し訳ありません。

[22Go]

with open("short.sequ.txt") as f:

max_len = 0
max_protein = ""

min_len = 10000000000000
min_protein = ""

def change(protein, seq,max_len,max_protein,min_len,min_protein):
seq_len = len(seq)
if seq_len == 0:
return

else:
print([seq_len])

if seq_len > max_len:
seq_len = max_len
protein = max_protein
return seq_len and protein

if seq_len < min_len:
seq_len = min_len
protein = min_protein
return seq_len and protein

kekka = change(protein, seq,max_len,max_protein,min_len,min_protein)

a_line = ""
a_seq = ""

for line in f:
strip_line = line.rstrip()
if strip_line.startswith(">"):
a_line = strip_line
change(a_line, a_seq)
a_seq = ""
print(strip_line+'\n')
else:
A = [len(a_seq+strip_line)]
a_seq += strip_line
msg1 = "最大アミノ酸エントリ："
msg2 = "長さ："
msg3 = "最小アミノ酸エントリ："
print(msg1,max_protein,msg2,max_len)
print(msg3,min_protein,msg2,min_len)

このようにしてみました

[otn]

def add(x,y)のような、引数の2つの数値の和を返す関数を書いて、それを使って1+2の結果を求める
というプログラムを書いてみてください。

[22Go]

def add(x,y):
a = x+y
return a

n = add(1,2)
print(n)

こういったものでしょうか？

[otn]

そうですね。
ではなぜ、
change(a_line, a_seq,max_len,max_protein,min_len,min_protein)
などと書いているのか？

[22Go]

addの中身をifの条件にtrueのときだけ変えるには
change(a_line, a_seq,max_len,max_protein,min_len,min_protein)
にするしかないのではと考えてしまうのですが、、、、

[otn]

changeの返り値を代入しないと意味がありません。
やはり、関数と返り値についての理解が足りていません。

Answer 2

現在のコードのままという制約が厳しすぎて回答難しいのですが…
関数の引数としてリストや辞書を渡して、その要素値を変更する手法はどうででしょうか？

def change(protein, seq, ret):
    # 略
    if seq_len > max_len:
        ret['max_len'],ret['max_protein'],～ = seq_len, protein, ～

ret = {}
change(～, ret)

Comments

[22Go]

できるだけこのコードのままでありたいと思っているのですが手を加えても全然大丈夫です！
リストや辞書も考えたのですがうまく応用ができなくて、、、やり方を教えていただけたらと思います。
宜しくお願い致します。

[22Go]

with open("sequence.fasta.txt") as f:

max_len = 0
max_protein = ""

min_len = 10000000000000
min_protein = ""

def change(protein, seq, ret):
max_len,max_protein,min_len,min_protein
seq_len = len(seq)
if seq_len == 0:
return

else:
a = [seq_len]
print(a)

if seq_len > max_len:
ret["max_len"],ret["max_protein"],ret["min_len"],ret["min_protein"] = seq_len,protein,seq_len,protein
K = {}
change(protein, seq, ret)

else:
if seq_len < min_len:
min_len,min_protein = seq_len,protein

a_line = ""
a_seq = ""

for line in f:
strip_line = line.rstrip()
if strip_line.startswith(">"):
a_line = strip_line
change(a_line, a_seq)
a_seq = ""
print(strip_line+'\n')
else:
#A = [len(a_seq+strip_line)]
a_seq += strip_line
msg1 = "最大アミノ酸エントリ："
msg2 = "長さ："
msg3 = "最小アミノ酸エントリ："
print(msg1,max_protein,msg2,max_len)
print(msg3,min_protein,msg2,min_len)

このような感じで書き換えましたが
change(a_line, a_seq)
TypeError: change() missing 1 required positional argument: 'ret'
のエラーが出てしまいました。
これはどのように対処したらいいでしょうか

[can110]

回答コードと異なるようです。
具体的にはchangeにretを渡していないようです。
回答内容をご理解されてからご自身のコードに組み込むようお願いします。

[22Go]

with open("sequence.fasta.txt") as f:

max_len = 0
max_protein = ""

min_len = 10000000000000
min_protein = ""

def change(protein, seq, ret):
seq_len = len(seq)
if seq_len == 0:
return

else:
print([seq_len])

if seq_len > max_len:
ret["max_len"],ret["max_protein"],ret["min_len"],ret["min_protein"] = seq_len,protein,seq_len,protein

if seq_len < min_len:
min_len,min_protein = seq_len,protein

ret = {}
change(protein, seq, ret)
a_line = ""
a_seq = ""

for line in f:
strip_line = line.rstrip()
if strip_line.startswith(">"):
a_line = strip_line
change(a_line, a_seq)
a_seq = ""
print(strip_line+'\n')
else:
#A = [len(a_seq+strip_line)]
a_seq += strip_line
msg1 = "最大アミノ酸エントリ："
msg2 = "長さ："
msg3 = "最小アミノ酸エントリ："
print(msg1,max_protein,msg2,max_len)
print(msg3,min_protein,msg2,min_len)

[22Go]

これであってますでしょうか？
エラーの内容としては
change(protein, seq, ret)
NameError: name 'protein' is not defined
なのですが、最初のdef changeでproteinは定義されていることにはならないのでしょうか？

Answer 3

全てのデータを関数のパラメータとして渡すようにすればいい……というような話でしょうか。