hpprc/llmjp-kaken huggingface dataset:
進化型多目的局所探索の最適なアルゴリズム構造の決定に関する研究

複数の競合する目的を持つ多目的最適化問題を解く場合では,目的間のトレードオフ関係を明確にするために,パレート最適解と呼ばれる解の集合を求めることが重要である.パレート最適解集合を求めることで,目的関数空間内のパレートフロントが得られ,目的間のトレードオフが明確になる.進化型多目的最適化アルゴリズムでは,多点探索という進化計算の特徴を利用することで,アルゴリズムの1回の実行で全てのパレート最適解の探索が行われる.本研究では,進化型多目的最適化アルゴリズムの計算効率を大きく向上させるための局所探索の組み込みに関する研究調査を行い,進化型多目的局所探索の実装に関する重要な知見が得られた.本研究の最終的な目的は,進化型多目的最適化アルゴリズムに局所探索を組み込んだ進化型多目的局所探索アルゴリズムの自動設計である.この最終目的を実現するために,本年度は,4年間の研究期間の初年度として,以下の研究を試みた.(1)変数の数や目的の数が異なる様々な多目的組合せ最適化問題の問題例を作成し,進化型多目的最適化アルゴリズムを適用することで,問題の特徴や多目的最適化の難しさを調べる.(2)研究代表者の長年の経験と知識に基づき,高性能と思われる構造を持つ進化型多目的局所探索アルゴリズムを実装し,数値実験により性能評価を行う.(3)進化型多目的局所探索アルゴリズムの様々なパラメータを記号列とするメタレベルの遺伝的アルゴリズムを実装し,パラメータ設定の最適化を行う.(4)交叉,突然変異,局所探索,世代更新など,進化型多目的局所探索を構成する操作の順序を記号列とするメタレベルの遺伝的アルゴリズムを実装し,最適なアルゴリズム構造の探索を行う.上記(1)に関しては,ナップサック問題や距離最小化問題などの問題例を作成し,NSGA-IIやMOEA/Dなどの進化型多目的最適化アルゴリズムの性能評価を行うことができた.上記(2)と(3)に関しても,局所探索と遺伝的探索の順序に関するアルゴリズム構造の調査,および,局所探索確率や局所探索繰り返し数の最適化などを行うことができ,国際会議に論文を投稿できる結果が得られた.上記(4)に関しては,単純な単一目的最適化問題に対して,個体のランダム生成や交叉操作,突然変異操作の適用順序と探索速度の関係調査など,準備段階の研究を行った.ただし,最適な進化型多目的局所探索アルゴリズムの構造探索は,2年目以降の課題として残った.なお,関連した研究として,局所探索を行う方向の動的調整メカニズムの検討も行った.進化型多目的最適化アルゴリズムに局所探索を組み込んだ進化型多目的局所探索アルゴリズムの自動設計という最終的な長期目標に向けて,本年度は,以下のことを行った.多目的最適化問題の定式化と多目的最適化アルゴリズムの性能調査:目的の数と変数の数を自由に設定することのできる多目的最適化問題として,高次元決定変数空間に複数の点を配置することで多目的距離最小化問題の定式化を行った.また,この問題において,複数の点の配置を規則的に行うことで,線形空間内で個体群の進化の挙動を観察することができることを示した.多目的距離最小化問題に対する数値実験により,変数の数を大きく増加させると個体群の多様性が大きく減少する傾向がみられることを示した.さらに,多目的ナップサック問題に対して,目的間に様々な相関関係や依存関係を持つ例題を生成し,多目的最適化問題の性能を調べた.多数目的最適化問題に対する進化型多目的最適化アルゴリズムの性能調査:目的関数の数が4以上である多数目的最適化は,進化計算分野で大きな注目を集め,急激に進展している研究分野である.この分野に貢献するために,様々な異なる性質を持つ多様な多数目的ナップサック問題に対して,代表的な進化型多目的最適化アルゴリズムを適用し,個々のアルゴリズムの探索挙動を調べた.多目的ファジィ遺伝的機械学習への応用:多目的ファジィ遺伝的機械学習問題に対して,様々な進化型多目的最適化アルゴリズムを適用し,個々のアルゴリズムの探索挙動の調査を行った.また,多目的ファジィ遺伝的機械学習問題に対して獲得された多数の非劣解集合から最適な1個の解を選択する方法の検討も行った.集合最適化問題としての再定式化:多目的最適化問題の解集合に対する複数の評価基準を目的関数とする多目的集合最適化問題の定式化を行った.本研究の最終的な目的は,進化型多目的最適化アルゴリズムに局所探索を組み込んだ進化型多目的局所探索アルゴリズムの自動設計である.この最終目的を実現するためには,交叉や突然変異などの遺伝的操作の選択,選択圧の調整,局所探索の実装,それらの最適な組合せなど,アルゴリズム開発に関連した幅広い研究が必要となる.また,個々の問題に適したアルゴリズム構造が異なるので,幅広い応用分野における多目的最適化問題に関する研究も必要となる.その結果として,以下のような幅広い多様な研究成果が得られている.1.個々の世代で異なるアルゴリズム構造と異なる遺伝的操作を持つことが可能な進化計算アルゴリズムのメタレベルでのオフライン最適化手法を提案した.この研究により,進化計算の序盤,中盤,終盤および最終盤など,探索の異なる曲面において異なるアルゴリズム構造や異なる遺伝的操作の必要性を示した.さらに,異なる問題では,全く異なるアルゴリズム構造や遺伝的操作が必要であることも示した.2.上記1の研究を進化型多目的局所探索

アルゴリズムの設計に拡張した.メタレベルで最適化を行うアルゴリズムの構造としては,個々の世代における遺伝的探索と局所探索の適用順序を用いた.また,局所探索に関する局所探索の適用確率や局所探索で調べられる近傍解の数などのパラメータもメタレベルでの最適化の対象とした.このような設定により,個々の世代で異なるアルゴリズム構造を持つ進化型多目的局所探索アルゴリズムのオフライン最適化を行った.3.代表的な多目的組合せ最適化問題である多目的ナップサック問題を用いて,様々な多目的最適化アルゴリズムの探索挙動の調査を行った.また,ファジィシステムの多目的最適化やネットワークの多目的設計などの応用問題への進化型多目的最適化アルゴリズムの応用に関する研究も行った.本研究の最終的な目的は,進化型多目的最適化アルゴリズムに局所探索を組み込んだ進化型多目的局所探索アルゴリズムの自動設計である.アルゴリズムの完全な自動設計は今後の課題となるが,今年度は,2目的ナップサック問題に対して,アルゴリズム構造の最適化も含めた進化型多目的局所探索アルゴリズムの設計方法の提案を国際会議IEEE CEC 2015で行った.本研究で提案した設計方法の特徴は,以下のようにまとめることができる.(1)局所探索を適用する確率,局所探索を継続する時間,局所探索適用個体を選択する場合での選択圧などのパラメータだけでなく,世代更新,交叉,突然変異,局所探索をどのような順序で行うかというアルゴリズム構造を自動的に決定することができる.(2)進化の異なる段階で異なるアルゴリズム構造と異なるパラメータを使うことができるように,アルゴリズムの自動設計を世代毎に行うことができる.また,世代毎に異なる構造の探索メカニズムを用いるアルゴリズムの有効性を検討するために,幅広く用いられている単一目的関数最適化問題に対して自動設計されたアルゴリズムの一般性すなわち汎化能力の検討を含めた論文をオープンアクセス論文誌であるSpringerPlus誌で発表した.この論文では,ある特定のテスト問題に対して設計されたアルゴリズムは,類似した性質のテスト問題に対して高い性能を持つという予想された結果だけではなく,探索の序盤で幅広くランダムに探索を行うことは,単純な単峰性関数から複雑な多峰性関数まで多様なテスト問題に対して有効であるという結果が得られた.複数の競合する目的を持つ多目的最適化問題を解く場合では,目的間のトレードオフ関係を明確にするために,パレート最適解と呼ばれる解の集合を求めることが重要である.パレート最適解集合を求めることで,目的関数空間内のパレートフロントが得られ,目的間のトレードオフが明確になる.進化型多目的最適化アルゴリズムでは,多点探索という進化計算の特徴を利用することで,アルゴリズムの1回の実行で全てのパレート最適解の探索が行われる.本研究では,進化型多目的最適化アルゴリズムの計算効率を大きく向上させるための局所探索の組み込みに関する研究調査を行い,進化型多目的局所探索の実装に関する重要な知見が得られた.本研究の最終的な目的である進化型多目的最適化アルゴリズムに局所探索を組み込んだ進化型多目的局所探索アルゴリズムの自動設計までは実現できていないが,その準備となるアルゴリズム開発から応用問題の調査まで,幅広い研究を実施できており,おおむね順調に進展していると言える.特に,個々の世代で異なるアルゴリズム構造を持つことができる進化型多目的局所探索
ラグランジュ型非ダルシー流解析モデルの改良と応用～現場適用を目指して～

本研究は,ダルシー則に従わない浸透流(非ダルシー流)の解析が不可欠な現場の問題に取り組むために研究代表者が開発してきたダルシー非ダルシー流解析モデルについて,より一般的な土壌への適用可能性を検証して当該モデルの適用可能性を拡大し汎用性を向上させるとともに,非ダルシー流の解析が不可欠な現場の問題である土砂災害の発生に寄与する山地斜面への降雨浸透流や暗渠周りの浸透・排水流解析へ適用することを目的とする.平成30年度は,前年度に改良したダルシー非ダルシー流解析モデルを斜面土層における降雨浸透流解析への適用可能性について検討した.土層に充填する試料として5種類(砂利大(50%粒径:14mm),砂利小(50%粒径:8.3mm),ガラスビーズ(粒径:5mm,2mm,1mm))を用意した.これらの土層に対して斜面勾配を10%と20%,降雨強度を50mm/hと100mm/hの計4条件の降雨浸透実験を実施するとともに,本解析モデルによる数値解析を実施し,土層実験結果に対する数値解析結果(浸透流場・浸透流量)の再現性を調べ,本モデルの適用可能性を検証した.土層末端での浸透流速からレイノルズ数を算出すると,粒径,斜面勾配,降雨強度が大きいほどレイノルズ数も大きくなり,砂利大土層の斜面勾配20%・降雨強度100mm/hで最大値(およそ10)を取り,これはダルシー流から非ダルシー流に遷移するといわれている範囲の上限値であった.数値解析による浸透流量の再現性については,離散化手法に採用したMPS法の解像度(粒子径)を1cm,0.5cm,0.2cmで比較した結果,粒子径を1cmから0.2cmすなわち1/5倍すると誤差が最大で1cm程度少なくなったが,反対に,計算時間は2倍に増えた.また,浸透流場全体でみると,本モデルでは考慮していない不飽和領域が観察された.したがって,このような問題への適用には,不飽和領域も解析できるよう本モデルの改良が必要である.ダルシー非ダルシー流解析モデルについて,斜面土層における降雨浸透流解析への適用可能性を検討した.不飽和領域の考慮という課題が残ったため,「(3)やや遅れている」と判断した.本研究において目的としている非ダルシー流の解析が不可欠な現場の問題である暗渠周りの浸透・排水流解析へ適用する.対象とする現象を再現するような土層を作成し,土層実験結果をもとに本モデルによる解析結果の再現性の観点から,本モデルの妥当性を検証する.また,本モデルの不飽和領域への拡張についても検討する.本研究は,ダルシー則に従わない浸透流(非ダルシー流)の解析が不可欠な現場の問題に取り組むために研究代表者が開発してきたダルシー非ダルシー流解析モデルについて,不均質な土壌や成層化した土壌といったより一般的な土壌への適用可能性を検証し,当該モデルの適用可能性を拡大し汎用性を向上させるとともに,非ダルシー流の解析が不可欠な現場の問題である土砂災害の発生に寄与する山地斜面への降雨浸透流や暗渠周りの浸透・排水流解析へ適用することを目的とする.平成29年度は,すでに研究代表者が開発したダルシー非ダルシー流解析モデルについて,より一般的な土壌への適用可能性を検証し,当該モデルの汎用性向上に取り組んだ.より一般的な土壌として,不均質な土壌および成層化した土壌を対象とした.不均質な土壌として,粒径5mmと2mmのガラスビーズを一定の割合で混同した土壌,成層化した土壌として,粒径5mmと2mmのガラスビーズを層状に配置した土壌を用い,これらを充填した土層を作成し浸透実験を実施し,土層浸透実験結果に対する数値解析結果の再現性の観点から,当該モデルの妥当性を検証した.その結果,当該モデルは一般的な土壌に対しても適用可能性が高いことが示された.また,当該モデルの運動方程式において採用した非線形抵抗則であるフォルヒハイマー則における係数のモデル化あるいは与件が重要であり,既往係数モデルを用いるだけでなく,フィッティングパラメータとして与件することで当該モデルの再現性が向上することが分かった.ダルシー非ダルシー流解析モデルについて,より一般的な土壌への適用可能性を検証し,当該モデルの適用可能性を拡大し汎用性を向上できたため,「(2)おおむね順調に進展している」と判断した.本研究は,ダルシー則に従わない浸透流(非ダルシー流)の解析が不可欠な現場の問題に取り組むために研究代表者が開発してきたダルシー非ダルシー流解析モデルについて,より一般的な土壌への適用可能性を検証して当該モデルの適用可能性を拡大し汎用性を向上させるとともに,非ダルシー流の解析が不可欠な現場の問題である土砂災害の発生に寄与する山地斜面への降雨浸透流や暗渠周りの浸透・排水流解析へ適用することを目的とする.

平成30年度は,前年度に改良したダルシー非ダルシー流解析モデルを斜面土層における降雨浸透流解析への適用可能性について検討した.土層に充填する試料として5種類(砂利大(50%粒径:14mm),砂利小(50%粒径:8.3mm),ガラスビーズ(粒径:5mm,2mm,1mm))を用意した.これらの土層に対して斜面勾配を10%と20%,降雨強度を50mm/hと100mm/hの計4条件の降雨浸透実験を実施するとともに,本解析モデルによる数値解析を実施し,土層実験結果に対する数値解析結果(浸透流場・浸透流量)の再現性を調べ,本モデルの適用可能性を検証した.土層末端での浸透流速からレイノルズ数を算出すると,粒径,斜面勾配,降雨強度が大きいほどレイノルズ数も大きくなり,砂利大土層の斜面勾配20%・降雨強度100mm/hで最大値(およそ10)を取り,これはダルシー流から非ダルシー流に遷移するといわれている範囲の上限値であった.数値解析による浸透流量の再現性については,離散化手法に採用したMPS法の解像度(粒子径)を1cm,0.5cm,0.2cmで比較した結果,粒子径を1cmから0.2cmすなわち1/5倍すると誤差が最大で1cm程度少なくなったが,反対に,計算時間は2倍に増えた.また,浸透流場全体でみると,本モデルでは考慮していない不飽和領域が観察された.したがって,このような問題への適用には,不飽和領域も解析できるよう本モデルの改良が必要である.ダルシー非ダルシー流解析モデルについて,斜面土層における降雨浸透流解析への適用可能性を検討した.不飽和領域の考慮という課題が残ったため,「(3)やや遅れている」と判断した.本研究において目的としている非ダルシー流の解析が不可欠な現場の問題である土砂災害の発生に寄与する山地斜面への降雨浸透流および暗渠周りの浸透・排水流解析へ適用する.これらの現象を再現するような土層を作成し,土層実験結果をもとに当該モデルによる解析結果の再現性の観点から,当該モデルの妥当性を検証する.本研究において目的としている非ダルシー流の解析が不可欠な現場の問題である暗渠周りの浸透・排水流解析へ適用する.対象とする現象を再現するような土層を作成し,土層実験結果をもとに本モデルによる解析結果の再現性の観点から,本モデルの妥当性を検証する.また,本モデルの不飽和領域への拡張についても検討する.
重症不全心に対する遺伝子導入とRNAiを組み合わせた新しい治療法の開発

我々は、ラット心不全心モデルに,sarcoplasmic reticulum Ca^<2+>ATPase (SERCA2a)遺伝子を導入することにより心機能と生存率が改善されると報告してきた。またSERCA2aの抑制蛋白であるPhospholamban(PLB)をアンチセンス法で抑制することでも、SERCA2a機能を亢進できる。このように我々はPLBとSERCA2aの比率(PLB/SERCA2a ratio)に着目し不全心の治療が可能であると考えている。しかし、従来のアデノウイルスによる遺伝子導入では小動物では心臓全体に対し治療が可能であるが、大動物やヒトの心臓全体に対しては充分な機能回復をもたらすことは困難であり、その方法は未だ確立されていない。最近のトピックスで、線虫でRNAi (RNA interference : RNA干渉)により遺伝子発現が抑制されていることが発見され、哺乳類でも同様のメカニズムがあり、小さなRNAを用いて遺伝子サイレンシングを行う方法が報告され非常に注目を浴びている。そこでRNAiと遺伝子導入を組み合わせて、より高率の良いしかも毒性の少ない不全心治療法を開発することが本研究の目的である。本年度は以前行っていたエックス線透視下の遺伝子導入法を超音波ガイド下で行えるようになった。またsiRNAの作製については、以前、他施設でPLBに有効なsiRNAを作製したが、in vitroでは有効であったが、in vivoでは無効であったという情報をえた。しかし、その後も生体内でPLBに有効なsiRNAの入手はかなわず今後の発展に期待したい。また、重症心不全に対するSRCA2aの有効性についてはさらに検討し報告を重ねた。現在、アデノウイルスを使用しないより安全な方法を検討中であり、臨床応用にむけさらなる進化が期待される。またSERCA2aの抑制蛋白であるPhospholamban (PLB)をアンチセンス法で抑制することでも、SERCA2a機能を亢進できる。このように我々は(PLBとSERCA2aの比率(PLB/SERCA2a ratio)に着目し不全心の治療が可能であると考えている。しかし、従来のアデノウイルスによる遺伝子導入では小動物では心臓全体に対し治療が可能であるが、大動物やヒトの心臓全体に対しては十分な機能回復をもたらすことは困難であり、その方法は未だ確立されていない。そこでRNAiと遺伝子導入を組み合わせて、より効率の良いしかも毒性の少ない不全心治療法を開発することが本研究の目的である。本年度は以前行なっていたエックス線透視下の遺伝子導入法をエコーガイド下で行なえるように現在調整中である。またiRNAの作製については、以前、他施設でPLBに有効なiRNAを作製したが、in vitroでは有効であったが、in vivoでは無効であったという報告を入手した。PLBに有効なiRNAについては現在設計中であるが、その情報もふまえたうえ慎重に検討する必要があると思われる。我々は、ラット心不全心モデルに,sarcoplasmic reticulum Ca^<2+> ATPase (SERCA2a)遺伝子を導入することにより心機能と生存率が改善されると報告してきた。またSERCA2aの抑制蛋白であるPhospholamban (PLB)をアンチセンス法で抑制することでも、SERCA2a機能を充進できる。

このように我々は(PLBとSERCA2aの比率(PLB/SERCA2a ratio)に着目し不全心の治療が可能であると考えている。しかし、従来のアデノウイルスによる遺伝子導入では小動物では心臓全体に対し治療が可能であるが、大動物やヒトの心臓全体に対しては十分な機能回復をもたらすことは困難であり、その方法は未だ確立されていない。最近のトピックスで、線虫でRNAi (RNA interference : RNA干渉)により遺伝子発現が抑制されていることが発見され、哺乳類でも同様のメカニズムがあり、小さなRNAを用いて遺伝子サイレンシングを行う方法が報告され非常に注目を浴びている。そこでRNAiと遺伝子導入を組み合わせて、より効率の良いしかも毒性の少ない不全心治療法を開発することが本研究の目的である。本年度は以前行なっていたエックス線透視下の遺伝子導入法を超音波ガイド下で行なえるようになった。またsiRNAの作製については、以前、他施設でPLBに有効なsiRNAを作製したが、in vitroでは有効であったが、in vivoでは無効であったという情報をえた。しかし、その後も生体内でPLBに有効なsiRNAの入手はかなわず今後の発展に期待したい。また、重症心不全に対するSERCA2aの有効性についてはさらに検討し報告を重ねた。現在、アデノウイルスを使用しないより安全な方法を検討中であり、臨床応用にむけさらなる進化が期待される。
脳画像解析と分子遺伝学的解析による強迫性障害の病態研究

強迫性障害患者(OCD)89名、健常者89名に対して、遺伝子解析を行った。セロトニントランスポーター遺伝子(5HTT -LPR)に関して患者群内の比較でL型と、うつ病合併OCDとの関連が示唆された。OCD群37名、健常群37名のMRI画像を撮像した。患者-健常二群間比較では、患者群で右後帯状回ほかの部位の体積減少が認められた。強迫症状のタイプと特定の脳の部位の体積との相関性も認められ、強迫症状のタイプごとに異なる神経基盤が責任部位である可能性が示唆された。九州大学病院精神科を受診した強迫性障害患者15名に対して血液DNAサンプルの採取をおこない、6名に対してMRI撮像をおこなった。このうち、血液サンプルの採取とMRI撮像の両方を実施した者は4名であった。血液サンプルを採取した15名(男性10名、女性5名)の平均年齢は33.3±9.7歳であった。症状のうちわけは、加害恐怖/確認強迫が4名、不潔恐怖/洗浄強迫が6名、縁起恐怖/儀式行為が1名、加害恐怖/確認強迫+不潔恐怖/洗浄強迫が2名、不潔恐怖/洗浄強迫+溜めこみが1名、不潔恐怖/洗浄強迫+縁起恐怖/儀式行為が1名であった。セロトニン・トランスポーター(SERT)遺伝子、ドーパミン・トランスポーター(DAT)遺伝子、および、ドーパミン受容体D4(DRD4)遺伝子の繰り返し配列多型(VNTR)に関して、強迫性障害の患者15名の遺伝子型を実験により決定した。対照者15名も同様に遺伝子型を決定した。MRI撮像をおこなった6名(男性3名、女性3名)の平均年齢は32.7±6.6歳であった。症状のうちわけは、加害恐怖/確認強迫が1名、不潔恐怖/洗浄強迫が3名、加害恐怖/確認強迫+不潔恐怖/洗浄強迫が1名、不潔恐怖/洗浄強迫+溜めこみが1名であった。今後の画像解析のための前処置として、得られた個人ごとの画像データについてSPM8を用いて、解剖学的標準化、平滑化等をおこない、体積情報を持たせた画像を作成した。今年度に同意を得たが、日程上の都合のため、まだMRI撮像が実施できていない対象者もおり、次年度以降症例数を増やしVBM解析による健常者との体積比較、症状と画像所見の関連等を検討していく予定である。また、次年度は、血液DNAサンプルも増やし、症状内容、治療反応性、発症年齢等と遺伝子型の関連を検討していくとともに、血液サンプル採取とMRI撮像両方の実施終了者を増やし、神経画像所見と遺伝子多型との関連についても検討していくことを計画している。九州大学病院精神科を受診した強迫性障害患者74名(男性28名、女性46名、そのうち平成23年度は59名が追加)に対して血液DNAサンプルの採取を行い、そのうち19名にMRI撮像を行った。また健常者74名(男性28名、女性46名、そのうち平成23年度は59名が追加)に対して血液DNAサンプルの採取を行い、そのうち23名にMRI撮像を行った。セロトニントランスポーター(SERT)遺伝子に関しては2種類(5HTT-LPR、5HTT-VNTR)、ドーパミントランスポーター(DAT)遺伝子(DAT-VNTR)および、ドーパミン受容体D4(DRD4)遺伝子の繰り返し配列多型2種類(D4-120、D4-48)に関して、強迫性障害の患者64名の遺伝子型を実験により決定した。対照者64名も同様に遺伝子型を決定した。年齢、性別、発症年齢、発症年齢が18歳未満か否か、症状タイプ別による分類を行い、遺伝子の多型との相関を調べている。最終年度は引き続き、強迫性障害患者と健常群のサンプル数を増やしていく。実験を更に進め遺伝子型の決定を急ぐと共に、症状タイプ別、治療反応性、発症年齢等と遺伝子型との関連を検討、解析をしていく。また神経画像所見と遺伝子多型との相関についても検討、解析を行う。MRI撮像を行ったのは、強迫性障害患者19名(男性9名、女性10名)及び健常者23名(男性9名、女性14名)であった。画像解析のための前処理として、得られた個人データ各々の画像データについて、SPMを用いて解剖学的標準化・平滑化等を行った。最終年度は、目標症例数である、患者50名・健常群50名まで症例を増やすと同時に、患者群と健常群の二群間においてVBM解析による体積比較、DTI画像の比較、症状や環境因子などのデータと画像所見との関連などを検討していく予定である。

強迫性障害患者(OCD)89名(男性34名、女性55名、平均年齢35.5歳(SD=10.4))、健常者89名(男性33名、女性56名、平均年齢35.9歳(SD=10.8))に対して血液DNAサンプルの採取を行った。セロトニントランスポーター遺伝子に関しては2種類(5HTT-LPRと5HTT-VNTR)、ドーパミントランスポーター遺伝子の繰り返し配列多型(DAT-VNTR)および、ドーパミン受容体D4遺伝子の繰り返し配列多型2種類(D4-120,D4-48)に関して遺伝子型を同定した。それぞれの遺伝子型、アリルに関して、患者-健常二群間比較及び、患者群内での臨床指標との相関解析を行った。いずれの遺伝子型、アリルの比較ともに健常群との有意差は認められなかった。患者群内において、うつ病合併/既往の有無による2群比較でアリルにおいて有意差を認め、L型がOCDや、うつ病合併のOCDとの関連が示唆された。3TesraMRI装置にて、37名のOCD群(男性14名、平均年齢34.4歳(SD=10.5))、37名の健常群(男性14名、平均年齢36.8歳(SD=10.8))のMRI画像を撮像した。SPMソフトウェアにて灰白質画像に前処理を施し、全脳的探索方法を用いて、患者-健常二群間比較及び、患者群内での臨床指標との相関解析を行った。二群間比較では、患者群で右中側頭回・左背側前頭前野・右後帯状回・右眼窩前頭前野・右縁上回・左下前頭回の体積減少が認められた。更に、DY-BOCS (dimensional yele brown obsessive compulsive scale)内の“攻撃性"スコアと右小脳との負の相関、“汚染"スコアと右島皮質との負の相関が認められ、強迫症状のタイプごとに異なる神経基盤が責任部位である可能性があり、OCDの生物学的多様性が示唆された。候補遺伝子の多型の解析、VBM解析を用いたMRIの形態画像解析および神経心理機能を含む臨床評価についておおむね順調に実施することができている。24年度が最終年度であるため、記入しない。今後は更に症例を増やし、脳画像解析と遺伝子の多型解析を進めていく。症状亜型や脳画像との相関についても多変量解析を行っていく必要がある。
アジア圏英語学習者の作文・発話の体系的修訂に基づく大規模校閲コーパスの開発と分析

3年間のプロジェクトでは,英語を母語とする専門校閲家を組織し,アジアの10か国・地域の320人の学習者による640本の作文に対して,1)ESL Composition Profile (Jacobs et al. 1981)に基づく5観点評点(内容・構成・語彙・言語使用・表記)と,2)その完全な校閲バージョンが用意された。プロジェクトで収集したデータは,ICNALE EditedEssaysと称する新たなモジュールにまとめられ,2018年3月までにダウンロード版およびオンライン検索版として全面公開された。この間,収集したデータを用いて30本以上の関連研究が発表・公刊された。アジア圏学習者の作文に対する大規模な修訂データを集める前提作業として,本年度は以下の事業を行った。1)既有データを用いた分析の実施と論文の執筆・講演の実施2)新規データ収集のための実施計画の策定と業者との打ち合わせ3)新規データ公開のためのオンラインシステムの枠組みの策定と業者との打ち合わせこのうち,1)については,関連する内容について,電子情報通信学会:思考と言語研究会(TL)(H27年12月),津田塾大学言語文化研究所「L2ライティングセミナー」(H28年3月)などで講演を行った。Learner Corpus Research 2015(H27年9月)などの海外で研究発表を行った。2)については,既有データの中から,詳細な費用見積もりを取り,データ収集量を300本にする計画を立てた。また,修訂者に同時に評点付けを依頼することで,修訂コーパスが同時に評価データ付き作文コーパスとしても使えるよう計画を立てた。最後に,3)については,並行する科研プロジェクト(基盤B)の枠組みにおいて,今後の修訂データ公開も見据えたオンラインでの検索システムの大幅改修を行い,品詞情報を用いた検索語の入力,品詞タグの入力支援,特定の学習者の習熟度レベル別の頻度データ(例えば,A2,b11,B12,B2レベルの日本人学習者が特定の語彙をどの程度使用しているか,変化の様子が一覧表示できる)と国際学習者の頻度データ(例えば,特定の語を日中韓学習者および母語話者がどの程度使っているかが一覧表示できる)のグラフィカル表示システムを考案・開発し,ICNALE Onlineの新バージョンとして実装した。初年度に計画を立て,2年目にデータ収集を実施し,3年目に分析結果をまとめるという基本計画で作業を進めており,初年度においては,予定された準備内容がほぼすべて完了した。これにより2年目の冒頭からデータ収集にかかれる体制が整った。◎概要:アジア圏国際英語学習者コーパス(ICNALE)の新規モジュールとなる校閲モジュールの開発を進める。◎本年度進捗:2016/4月に,先行研究レビューを踏まえ,本プロジェクト用の作文評価ルーブリックの作成。2016/月に日本・中国・韓国・台湾・タイの学習者作文(約80本*4=320本)の校訂に着手。同6月にJPNモジュールの先行処理データ(40本)の校訂・評価作業が完了。同6月に今後の作業進捗についての調整会議(評価者間分散を統制するための共通作文に対するキャリブレーション評価システムの構築と新規の作業委嘱)。および,JPNモジュール(残40本)の校訂・評価作業が完了。2016/8月に本年度データ収集が完了。2016/12月にICNALE(本体作文モジュール)の拡張データセットとしてアラブ首長国連邦(UAE)でのデータ収集が完了。2017/2月に学習者の作文データの分析の基礎資料となる「大学入試センター試験英語語彙表(19872014年)」の開発・公開が完了。◎アウトプット:口頭発表として,シンポジウム東アジアの言語研究とL2教育研究コーパスが拓く学際研究の領野(大学英語教育学会東アジア英語研究会特別シンポジウム),Grammatical and Contextual Correctness in L2Writing: An Analysis of the ICNALE-Proofread, A Newly Designed Parallel Corpus Including Learners' Original and Edited Essays (全国英語教育学会), 3つの作文評価法:ばらつきが小さく精度が高いのはどれか(統計数理研究所)ほか。著書・論文として『ベーシック応用言語学』(ひつじ書房)ほか。ELF圏6か国・地域のうち,約半数で校閲データ収集が完了。おおむね予定通りの進捗である。なお,収集データは年度末までに整理を終え,2017年4月中に公開予定。グローバル化の渦中にある日本の英語教育の質的向上を目指す上で喫緊の課題は,言語産出におけるモデルと学習者の乖離を客観的に把握し,両者間の距離の最小化を図ることである。この意味において母語話者校閲は有効な手段となる。筆者は日本人を含む各国学習者の英作文・英語発話を収集し,世界最大のアジア圏学習者コーパスICNALEを構築したが,予算的制約から,校閲データは包含できていなかった。本プロジェクトは,この欠落を埋めることを目的として開始された。プロジェクトで収集したデータは,ICNALE Edited

Essaysと称する新たなモジュールにまとめられ,2018年3月までにダウンロード版およびオンライン検索版として全面公開された。また,この間,収集したデータを用いて10本以上の関連研究が発表・公刊された。それらの研究では,(1)各国学習者ごとの作文の観点評価の比較,(2)観点評価と校閲量(校閲者による追加や削除の回数)の関係のモデル化,(3)評価者の判断の安定性検証などのテーマが扱われた。今回のプロジェクトで公開されたデータ従前において内外にほとんど類例のないもので,学習者研究の新たな基盤資料となるものである。3年間のプロジェクトでは,英語を母語とする専門校閲家を組織し,アジアの10か国・地域の320人の学習者による640本の作文に対して,1)ESL Composition Profile (Jacobs et al. 1981)に基づく5観点評点(内容・構成・語彙・言語使用・表記)と,2)その完全な校閲バージョンが用意された。プロジェクトで収集したデータは,ICNALE EditedEssaysと称する新たなモジュールにまとめられ,2018年3月までにダウンロード版およびオンライン検索版として全面公開された。この間,収集したデータを用いて30本以上の関連研究が発表・公刊された。2年目については,まず,修訂データの体系的取集を行う。すでに業者の選定は済ませており,今後,ライティング評価のマトリックスを示して評価者トレーイングを行った上で,データの収集に入る。およそ34か月程度で初期のデータ収集が完了する予定である。その後納品データの内容検査を行い,コーパス公開分析実施に向けた作業を継続して行う。本年度の予定1)残りの2か国・地域での校閲データ収集2)同じく評価データの収集3)収集データの公開のためのオンラインプラットフォーム開発4)関連研究をまとめた新しい論文の作成応用言語学修訂データの依頼業者が年度内に決まらず,発注が2年目になったため。すでに業者は決定したので,2年目の早期にデータ収集を開始できる予定である。
向社会的意思決定を支える神経基盤に関する文化差の解明

本研究は、向社会的な意思決定を支える神経基盤の文化差と、そうした文化差を生み出す要因について検討する。意思決定において他者に選択の余地を残そうとする向社会的意思決定は、他者への共感や視点取得に関わる脳部位と関連することが近年明らかにされている。一方で比較文化研究は、こうした共感に基づく向社会的意思決定は北米的パターンであり、東アジア人ではむしろ向社会的にふるまわなかった場合に生じる集団からの排斥といった社会的罰への恐怖から向社会的意思決定が生じると論じている。そこで本研究は、日本人と北米人の間で比較研究を行い、向社会的意思決定を支える神経基盤の文化差と、そうした文化差の説明要因を明らかにする。
マイクロコイル状カーボンファイバーを用いた超小型大変位計の開発

光学顕微鏡を用いた引張り測定実験において、我々は荷重としてカバーグラスを用いた。コイル状サンプルの一端は固定され、他端は荷重に接続された。光学顕微鏡は台と共に傾斜され、傾斜角からサンプルにかかる荷重を求めた。40°と50°傾斜させたときのコイル状カーボンの伸びが顕微鏡写真から測定された。剛性率Gは64nr^4w/d^4δで与えられる。但し、dはファイバーの直径、rはコイル径、nはピッチ数でδはその伸びである。コイルカーボンの測定値から得られた剛性率G=85GPaであった。バルクの電気抵抗率はバルク密度に依存して、1.0-0.2Ωcmの値を示した。一方、測定に用いたカーボンコイルの長さは約1mm、コイル径は8μm、ファイバーの直径は2μmであった。これらの値から計算されたカーボンコイル1本の電気抵抗率は1.09×10^<-2>Ωcmであった。長さの変化を測定するには抵抗変化、インピーダンス変化や光梃子の原理を用いた方法などが通常用いられている。しかし、通常これらの場合サンプルの大きさは数センチメートル以上で、数ミリメートル以下のサンプルに適用することは不可能である。またサンプルの大きさが大きい場合でもクラックのような局所的部分で変位が生ずるような場合、平均的な変位は測定できても局所的変化は測定することが困難であった。我々はこれまでにアセチレンの熱分解法により、マイクロコイル状カーボンファイバーを容易に生成させることに成功した。このファイバーは約3倍までは弾性的変形により伸び縮みが可能であるので、これらを用いた超小型変位計の開発を試みた。相異なる2種類の物質A,Bの境界で引っ張りにより破壊が生ずるような場合、その境界における変位の測定が必要となることがある。このときA,B自身も伸びるので、その両端で通常の方法で変位を測定すると平均的なものしか測定できない。そこで、コイル状カーボンファイバーを用いて、超小型変位計を作製し、1μm数mmの広い範囲の局所的変位を連続的に測定することを可能となった。当該ファイバーは長さは短いが断面積も小さいので、かなり大きな抵抗値を有すると思われるので、測定の際の接触抵抗の影響は少ないと考えられる。このファイバーの抵抗値を通常用いられるホイ-ストン・ブリッジ法を用いて測定することができた。1.コイル状カーボンファイバーの生成とその力学的特性我々はこれまでに微量の不純物活性剤の下で、アセチレンの熱分解法により、コイル径数μm、線径1μm以下、長さ数mmのコイル状カーボンファイバーを容易に生成させることに成功した。このファイバーは、約3倍までは弾性的変形により伸び縮みが可能であった。さらに、このコイル状ファイバーの伸びと荷重との関係を次の方法で測定することによって、剛性率Gを求めることができた。コイル状ファイバーの両端に適当なファイバーを接着する。スライドグラスの上に、その上をすべらすことのできる荷重用としてカバーグラス(約50mg)をのせ、ファイバーの一端をカバーグラス上で接着する。スライドグラスをこの顕微鏡上に静かにのせ、傾斜台をゆっくりと傾斜(最大85′)させる。傾斜角θから次の式よりそのときの荷重Wを求めることができる。W=W_0(sinθ-μcosθ)但し、W_0:カバーグラスの重さ、μ:最大静止摩擦係数。ファイバーの伸びと荷重Wとの関係から、剛性率Gは、但し、n:ピッチ数、r:コイル半径、d:ファイバーの直径、δ:ファイバーの伸び。2.コイル状カーボンファイバーの電気的特性長さ約2mmのコイル状カーボンファイバーの電気抵抗を4端子法により測定した結果、R=11.18KΩであった。このコイル状カーボンファイバーはピッチ間隔がなく、線が接触しあって、外見上はパイプ状となっていた。従って、パイプの外径と内径をSEM観察した結果、外径4.35μm、内径1.45μmであった。測定間のファイバーの長さは、1.35mmであったので、これらの値から電気抵抗率ρ=1.09×10^<-2>Ωcmと求められた。一方、コイルの線間の電気抵抗が大きいと考えるとファイバーの直径1.45μm、長さ12.31mmとなるので、電気抵抗は率ρ=1.50×10^<-4>
海氷・海洋結合系におけるアイスアルベドフィードバック

まず、今までに各国の砕氷船によって収集された南極海全域の海氷・海洋上層データと過去20年間のマイクロ波放射計による海氷密接度データから、海氷は開水面に入った熱によって主に融解し、海洋上層の熱と塩は概ねローカルバランスが成り立つことが示唆された。また、最も短波放射が大きい時期(12月)での海氷のアノマリーが、アルベドフィードバック効果によって次の結氷期にまで大きな影響を及ぼすことも示された。これらの知見から、開水面から入る熱によってのみ海氷が融解するとする簡略な海氷-海洋結合モデルを作成した。このモデルは、観測された海氷密接度と海洋混合層の関係、及び海氷融解の季節進行を非常によく説明できる。観測データとモデルのフィッティングから、海洋-海氷結合系のキーパラメーターである海洋海氷間熱交換係数も決めることができる。バルクでみたときこの値は全南極域で季節を問わずほぼ同様な値となることがわかった。さらにこの結合モデルは、海氷融解から結氷に至るまでの海氷・海洋場の気候値・年々変動を再現することができ、海氷融解期及び結氷期の年々変動は海氷-海洋結合系のアイスナルベドフィードバック効果によって決まる部分が大きいことが示唆された。すなわち、風の場が海氷の場に対して発散的であると、開水面より短波放射をより吸収し、海氷融解を促進し開水面を拡げ、ますます短波放射を吸収する、というポジティヴフィードバックが重要となる。風が収束的存場合は逆になる。しかし一方、大気も結合させた系で考えないと説明できない場合があることも示された。大気-海氷-海洋結合系士して、アイスアルベドフィードバック効果を考えることが今後の課題として提示された。本年度は南極海の融解・結氷間におけるアイスアルベドフィードバック効果の解析と南極ロス海での観測の解析を行った。過去20年間のマイクロ波放射計による海氷密接度データから、南極海の融解最盛期である12月の海氷のアノマリーは次の結氷期である4・5月の海氷のアノマリーと極めて高い相関があることが示される。このような融解期と結氷期の関係に対し、熱収支解析により以下のメカニズムを提出した。12月は最も短波放射が大きい時期である一方、海氷密接度のアノマリーも最も大きい時期である。その結果、海洋上層中に与えられる熱量が大きなアノマリーを持ち、そのアノマリーが海氷が消滅した後にも記憶され、次の結氷の早晩を生む。このメカニズムは一種のアイスアルベドフィードバックと言うことができ、海氷の年々変動を支配する一つの機構と考えられる。1999年1月においてアラスカ大学と共同で行った夏期の南極ロス海での海氷・海洋観測の結果を解析した。ビデオによる海氷密接度データとインテイクによる海洋上層の水温・塩分の関係を詳しく調べ、Ohshima et al.,(1998)による昭和基地沖の結果と同様な結果が得られた。すなわち、海洋内部領域では、海氷が開水面に入った熱によって主に融解し、熱と塩は概ねローカルバランスが成り立つことが示唆された。また、ポリニヤ域氷縁では日射により昇温した氷縁の外側の開水面に海氷が移流され効率的に海氷融解が生じていることが示唆された。さらに、これらの観測結果は簡単な海氷・海洋結合モデルでも表現できることが示された。本年度は、まず南極海の融解期の海氷・海洋上層データの解析を行い、簡略な海氷・海洋結合モデルとの比較・検討を行った。

南極海East Antarctica沖で広範囲にわたり、オーストラリアの砕氷船R/V Aurora Australisから観測された海洋上層の水温データと、人工衛星のマイクロ波放射計(SSM/I)による日毎の海氷密接度データから、融解期の海洋混合層温度と海氷密接度との関係を解析した。その結果、海域に関わらず海氷野内部では30km平均程度の空間スケールで見ると、海氷密接度は海洋混合層温度と明瞭な負の相関関係があり、海洋表面での大気からの熱のインプットが増加するに従い、負の相関関係の傾きが大きくなることがわかった。この関係は、開水面から入る熱によってのみ海氷が融解すると仮定した海氷-海洋結合モデルからうまく説明でき、またこの関係からバルク的に視たときの海洋-海氷間の熱交換係数も見積もることもできた。次に、上記の知見などに基づいて、海氷-海洋結合モデルの重要なパラメーターの値を決め、南極海全域を対象として、融解から結氷に至るまでの海氷・海洋場の気候値・年々変動の再現を試みた。その結果、開水面から入る熱によってのみ海氷が融解するとする簡略なモデルで、かなりの部分気候値や年々変動が再現できた。また、海氷・海洋結合系のアイスアルベドフィードバック効果が年々変動を決める重要な要因であることが示唆された一方、大気も結合させた系で考えることが年々変動を考える場合非常に重要であることが示された。本年度は、今までに収集された南極海全域での海氷・海洋上層データを統合して解析を行った。データは日本の「しらせ」、オーストラリアの「オーロラ」による東南極全域とアメリカの「パルマー」によるロス海で取得されたものである。これらのデータから、海氷密接度と海洋混合層の関係、及び海氷融解の季節進行は、開水面から入る熱によってのみ海氷が融解するとする簡略な海氷-海洋結合モデルで統一的に説明できることが示された。観測データとモデルのフィッティングから、海洋-海氷結合系のキーパラメーターである海洋海氷間熱交換係数も決めることができる。バルクでみたときこの値は全南極域で季節を問わずほぼ同様な値となることがわかった。この結合モデルは、海氷融解から結氷に至るまでの海氷・海洋場の気候値・年々変動を再現することができ、海氷融解期及び結氷期の年々変動は海氷-海洋結合系のアイスアルベドフィードバック効果によって決まる部分が大きいことが示唆された。すなわち、風の場が海氷の場に対して発散的であると、開水面より短波放射をより吸収し、海氷融解を促進し開水面を拡げ、ますます短波放射を吸収する、というポジティヴフィードバックが重要となる。風が収束的な場合は逆になる。しかし一方、大気も結合させた系で考えないと説明できない場合があることも示された。大気-海氷-海洋結合系として、アイスアルベドフィードバック効果を考えることが今後の課題として提示された。
抗etrone sulfataseモノクローナル抗体の作製とその特徴

子宮内膜癌組織では局所的な高エストロゲン状態が見られ,その原因はエストロンサルファターゼ(ES)やアロマターゼの異常と報告されている。しかしながらESの組織局在性は不明であり,その理由は良いモノクローナル抗体が存在しないためである。本研究では,抗ESモノクローナル抗体(MAb)を作成し,ESの組織局在性を明らかにすることを目的とした。正常胎盤10コからカラムクロマトグラフィーでモノピークを示し,かつES活性の高い蛋白分角を抽出した。BALB/Cマウスに抽出蛋白を計3回免疫し,hybridoma cellを作成。cloning後,ELISAとWestern blotでスクリーニングを行い12種の抗体を得た。正常者と先天性魚鮮癬患者の胎盤とで比較免疫組織化学染色を行った。なお胎盤使用には,患者の同意を得た。Aryl sulfatase C遺伝子を乳癌細胞(MCF-7,T470細胞株)に移入し,ESの酵素活性をみたところ,空VECTORのみの遺伝子移入細胞に比し,ESとDehydro epiandro steronesulfatase(DEAS)の両活性が高値を示した。これら細胞から抽出した蛋白と,MAbとは反応し,胎盤抽出抗原と反応が一致した。この結果は,作製したモノクローナル抗体がESとDEASの両酵素に反応するものと考えられた。正常子宮内膜と内膜癌組織においては,抗ESMAbは,上皮細胞のみに反応し,G_1,G_2,G_3,全ての内膜腺癌組織において腺癌細胞のみが反応し,間質には反応しなかった。以上の結果から,免疫組織化学に有用な抗ESMAbの作製に成功した。[方法]正常胎盤2コから鈴木らの方法に準じカラムクロマトグラフィーでモノピークを示す蛋白分画を抽出した。ついでBALA/cマウス2匹を抽出蛋白を計3回免疫し、electric fusion法でhybridom cellを作製。cloning後、ELISAとWestern blotでスクリーニングを行った。正常社と先天性魚鱗癬患者の胎盤とで比較免疫組織化学染色を行った。なお胎盤使用には、患者の同意を得た。[成績]胎盤より抽出した蛋白は最終カラムの酢酸エチルで抽出するとES活性、DHAS活性の各々高活性ピークが2峯性に分離した。第1ピークはDHAS高活性、第2ピークはES高活性であった。第2ピークのみを免疫し、ELISA,Western blotでスクリーニングし最終的に11種のclone cellを得た。凍結切片を用い11種のES MoAbにより免疫組織染色を行ったところ11種全部が魚鱗癬胎盤には染色されなかった。一方正常胎盤には11種中9種の抗体が染色可能であり絨毛表面の栄養膜細胞に組織局在性を認めた。ES活性とDHAS活性制御作用がMoAbに存在するか検討したが全11種と両酵素の活性抑制はみられなかった。現在Arylsufatase C遺伝子移入TD47細胞ができ、さらに抗体との反応性を検討中である。[結論]凍結切片に使用可能な抗ES MoAbが得られた。[方法]子宮内膜癌組織の凍結切片を作製。間接酵素抗体法により免疫組織染色を行う。同時に抗アロマターゼ抗体を用いて、両酵素の組織局在性を検討する。対照として正常胎盤を用いる。[成績]胎盤は、強くエストロンサルファターゼ抗体とアロマターゼ抗体により染色された。内膜癌では、腺癌細胞自体がよく染色されたが、間質の染色性はみられなかった。一方アロマターゼ抗体では、主に間質が染色され、G_3腺癌のみ上皮成分も染色された。

またエストロンサルファターゼ抗体は、正常内膜においても腺細胞自体がよく染色され、間質には染色性を認めなかった。[結論]子宮内膜癌では、エストロンサルファターゼは、主に腺上皮にアロマターゼは、主に間質に発現し、組織局在性に差を認めた。
無限グラフのハミルトン閉路とt-理想グラフ

本研究は無限グラフのハミルトン閉路とt-理想グラフを研究対象とした。t-理想グラフに関しては、比較的良い性質を持つクローフリーグラフのクラスに限定して研究を進めた。t-理想グラフに関する性質の多くは誘導路とよばれる構造と深く関係する。またハミルトン閉路に関する多くの問題では、グラフの切断集合が様々な必要条件を与える。そこでクローフリーグラフの誘導路と切断集合の関係を調べることが自然な流れとなった。前年度の研究では、クローフリーグラフの2頂点を結ぶ誘導路の本数を、クリーク切断集合の言葉で記述することに成功した。そこで今年度は互いに素な2つの頂点集合間の誘導路の本数を考えた。前年度の研究は2つの集合がともに位数1である場合に相当するので、本年度の研究はその拡張にあたる。一方2つの頂点集合が与えられたとき、それらをそれぞれ1頂点に縮約すると、前年度の研究成果が適用できる可能性が生じる。ただし初めのグラフがクローフリーグラフであっても、頂点を縮約することにより得られるグラフはクローフリーグラフであるとは限らない。ここに大きな困難があった。そこで我々はまず、前年度の研究成果をクローフリーグラフよりも大きいクラスに拡張することから始めた。具体的には、対象となる2頂点がクローの中心となっていても、そのクローフリー条件の欠如を補う条件を導入することで、前年度と同じ帰結を得ることに成功した。この新しい結果はクローフリーグラフの2つの頂点集合を縮約したグラフに適用できることが判明し、これを用いて2つの集合間を結ぶ誘導路の本数をクリーク切断集合の言葉で記述することに成功した。本年度は、クローフリーグラフにおいて指定された2頂点を通る誘導サイクルが存在するための条件を調べた。指定された非隣接2頂点x,yに対して、それらを分離するクリークが存在すれば、求める誘導サイクルが存在しないことは自明に分かる。しかし我々はそれだけが唯一の障害であり、そのようなクリークが存在しなければ常にx,yを通る誘導サイクルがあるのではないかと予想した。そこで調べてみると、xとyの距離が2の場合に反例が見つかった。一方xとyの距離が3以上あれば、予想は正しいことが分かった。また距離が2の場合にも、求める誘導サイクルの存在を妨げる条件を完全に決定することもできた。一般のグラフでは、指定された2頂点を分離するクリーク有無を答える多項式時間アルゴリズムが存在し、一方指定された2頂点を通る誘導サイクルの有無を問う問題はNP-完全である。本結果はクローフリーグラフでは後者がクラスPに属していることを示しており、意義深い。このように本年度は当初の目的を達成することができた。