日韩精品一区二区三区在线视频放-无码中文字幕V?一区二区-成年片免费观看视频-国内少妇人妻丰满av-国产精品中文字幕免费观看-亚洲成人久久一区二区三区-国内少妇偷人精品视频无缓冲-一区二区国产精品日本一区二区三区在线网

ARTICLE DETAIL

資訊詳情

深耕商務(wù)建站與企業(yè)官網(wǎng)運(yùn)營(yíng)的一線實(shí)戰(zhàn)洞察。

開(kāi)放研究工作流實(shí)操指南:從數(shù)據(jù)版本控制到可復(fù)現(xiàn)論文

開(kāi)放研究工作流實(shí)操指南:從數(shù)據(jù)版本控制到可復(fù)現(xiàn)論文 不少朋友看到“OpenResearch”這個(gè)詞第一反應(yīng)是某個(gè)具體的研究項(xiàng)目或者是某個(gè)高校實(shí)驗(yàn)室的代號(hào)。其實(shí)往大了說(shuō)它更像是一種正在席卷全球科研圈的方法論和工具鏈組合用開(kāi)放源代碼的思路去做研究把論文、數(shù)據(jù)、代碼、實(shí)驗(yàn)記錄統(tǒng)統(tǒng)攤開(kāi)來(lái)曬在陽(yáng)光下。往小了說(shuō)它也是一批具體開(kāi)源項(xiàng)目的集合幫助你從文獻(xiàn)管理、實(shí)驗(yàn)記錄到數(shù)據(jù)分析和論文寫作形成一條完全透明、可復(fù)現(xiàn)的流水線。這篇文章不只是解釋概念我打算把它當(dāng)作一份“實(shí)踐地圖”來(lái)寫聊聊開(kāi)放研究到底解決了什么痛點(diǎn)有哪些核心工具能直接上手以及我踩過(guò)幾次坑之后總結(jié)出的避坑要點(diǎn)。無(wú)論你是研究生、剛進(jìn)實(shí)驗(yàn)室的科研助理還是純粹對(duì)科研過(guò)程好奇的技術(shù)愛(ài)好者照著這篇文章的路徑走一遍都能搭起一套屬于你自己的開(kāi)放研究工作流。1. OpenResearch到底在“開(kāi)”什么先聊一個(gè)最基本的問(wèn)題開(kāi)放研究開(kāi)放的核心對(duì)象到底是什么很多人以為開(kāi)放研究就是“論文免費(fèi)下載”那只是冰山一角。開(kāi)放研究真正想撬動(dòng)的是整個(gè)知識(shí)生產(chǎn)鏈條的透明度。鏈條的第一環(huán)是研究問(wèn)題和假設(shè)第二環(huán)是實(shí)驗(yàn)設(shè)計(jì)和數(shù)據(jù)采集第三環(huán)是數(shù)據(jù)分析和代碼第四環(huán)是論文寫作和同行評(píng)審最后一環(huán)才是出版和傳播。傳統(tǒng)模式下除了最后一環(huán)的最終論文前四環(huán)基本都鎖在抽屜里審稿人看不到、同行看不到、社會(huì)公眾更看不到。OpenResearch的思路是把前四環(huán)也一并打開(kāi)。這帶來(lái)的直接好處有三層。第一層是可復(fù)現(xiàn)性。劍橋大學(xué)曾經(jīng)做過(guò)一項(xiàng)統(tǒng)計(jì)心理學(xué)和經(jīng)濟(jì)學(xué)領(lǐng)域有相當(dāng)比例發(fā)表過(guò)的論文在別人手里無(wú)法復(fù)現(xiàn)出同樣的結(jié)果。不是因?yàn)樵旒俣且驗(yàn)樵紨?shù)據(jù)、分析腳本、參數(shù)設(shè)置根本沒(méi)公開(kāi)別人無(wú)從下手。第二層是協(xié)作效率。當(dāng)數(shù)據(jù)和代碼托管在公開(kāi)平臺(tái)上任何人發(fā)現(xiàn)bug或者想到更好的分析方法可以直接提Pull Request像給開(kāi)源軟件貢獻(xiàn)代碼一樣給論文貢獻(xiàn)分析這種協(xié)作模式在傳統(tǒng)科研里幾乎不可想象。第三層是學(xué)術(shù)公平。發(fā)展中國(guó)家的小實(shí)驗(yàn)室訂閱不起頂刊數(shù)據(jù)庫(kù)但開(kāi)放研究讓核心數(shù)據(jù)和預(yù)印本免費(fèi)可得等于把知識(shí)壁壘拉低了一截。用一個(gè)生活化的比喻來(lái)說(shuō)傳統(tǒng)科研像私房菜大廚把菜端上來(lái)你只能看到成品色香味但用了什么火候、什么調(diào)料、中間有沒(méi)有翻車全憑大廚一面之詞開(kāi)放研究像“云廚房直播做菜”從洗菜切菜到出鍋裝盤全程可見(jiàn)你想學(xué)還能直接拿同一份菜譜在自己廚房里復(fù)刻一遍。這樣“做菜過(guò)程”本身也成了知識(shí)資產(chǎn)而不是只有“那盤菜”才算成果。所以當(dāng)你聽(tīng)到OpenResearch這個(gè)概念時(shí)不要只把它理解成一堆開(kāi)放獲取期刊的名單它本質(zhì)上是一種科研生產(chǎn)關(guān)系的重構(gòu)。它不要求每個(gè)研究者都變成程序員但它要求每個(gè)研究者具備“把過(guò)程記錄下來(lái)并共享出去”的意識(shí)和能力。1.1 從Open Access到Open Research的進(jìn)化邏輯要理解OpenResearch得先分清兩個(gè)經(jīng)常被混為一談的詞Open Access開(kāi)放獲取和Open Research開(kāi)放研究。開(kāi)放獲取是開(kāi)放研究的前置階段它關(guān)注的核心是“論文本身能不能免費(fèi)讀”。最有名的行動(dòng)是柏林宣言和布達(dá)佩斯開(kāi)放獲取倡議它們推動(dòng)了數(shù)千種期刊轉(zhuǎn)向開(kāi)放獲取模式。國(guó)內(nèi)讀者最熟悉的例子是Sci-Hub雖然它的法律爭(zhēng)議很大但它把“論文應(yīng)該免費(fèi)讀”這個(gè)訴求推到了大眾面前。開(kāi)放獲取解決的是知識(shí)消費(fèi)端的公平問(wèn)題也就是“讀得起”。開(kāi)放研究則往前走了一大步它關(guān)心“研究過(guò)程能不能公開(kāi)查”。光讀論文還不夠還得能拿到數(shù)據(jù)、腳本、實(shí)驗(yàn)材料清單、倫理審批記錄、同行評(píng)審意見(jiàn)甚至包括“被斃掉”的失敗實(shí)驗(yàn)。最典型的表現(xiàn)形式就是注冊(cè)報(bào)告Registered Report研究者在做實(shí)驗(yàn)之前先把研究問(wèn)題和分析方法寫成預(yù)注冊(cè)文檔提交給期刊期刊先審方案方案通過(guò)就原則上接受論文不管實(shí)驗(yàn)結(jié)果是什么。這種做法直接打擊了所謂的“文件抽屜效應(yīng)”——很多陰性結(jié)果因?yàn)椴弧捌痢倍绘i在抽屜里導(dǎo)致整個(gè)領(lǐng)域的知識(shí)體系被陽(yáng)性結(jié)果污染。從開(kāi)放獲取到開(kāi)放研究本質(zhì)上是把科研的評(píng)價(jià)單元從“最終論文”前移到“整個(gè)流程”。這種進(jìn)化不是某一個(gè)機(jī)構(gòu)推動(dòng)的而是一系列基礎(chǔ)設(shè)施成熟后的必然結(jié)果Git和GitHub讓代碼共享毫無(wú)門檻Zenodo和Figshare讓數(shù)據(jù)和附件有了穩(wěn)定的DOI編號(hào)Jupyter Notebook和R Markdown讓“分析代碼圖表解釋文字”能揉在同一個(gè)文檔里Overleaf讓LaTeX寫作協(xié)作變得像Google Docs一樣順滑。我自己的體會(huì)是很多人一聽(tīng)“開(kāi)放”兩個(gè)字就緊張覺(jué)得這是道德綁架好像不把全部家底掏出來(lái)就不配做科研。其實(shí)真不是這樣。開(kāi)放是有梯度的你可以在論文投稿時(shí)附帶數(shù)據(jù)也可以把代碼托管到GitHub還可以更進(jìn)一步托管整個(gè)實(shí)驗(yàn)環(huán)境比如Docker鏡像一步一步來(lái)只要比“啥都不開(kāi)”前進(jìn)一小步整個(gè)科研生態(tài)都會(huì)因此受益。1.2 開(kāi)放研究適合哪些人和團(tuán)隊(duì)回到實(shí)踐層面我接觸到的情況是會(huì)認(rèn)真搭建開(kāi)放研究工作流的人通常集中在這么幾類第一類是博士和青年教師。他們面臨最直接的考核壓力論文數(shù)量、引用量都是硬指標(biāo)。開(kāi)放研究對(duì)他們而言不僅是“道德選擇”更是“傳播策略”把數(shù)據(jù)和代碼公開(kāi)的論文引用率普遍比同類封閉論文高出不少這在多個(gè)實(shí)證研究里都被驗(yàn)證過(guò)。第二類是數(shù)據(jù)密集型的計(jì)算科學(xué)團(tuán)隊(duì)。生物信息、地理信息、計(jì)算化學(xué)、機(jī)器學(xué)習(xí)這些領(lǐng)域研究成果本身就是代碼和數(shù)據(jù)閉著反而會(huì)拖累團(tuán)隊(duì)內(nèi)部協(xié)作。第三類是真正想讓研究產(chǎn)生社會(huì)影響的人。比如做環(huán)境監(jiān)測(cè)、公共衛(wèi)生、教育政策的研究者他們的成果要被地方政府或公益組織采用首要前提就是對(duì)方能看懂并驗(yàn)證你的數(shù)據(jù)和分析過(guò)程封閉論文基本做不到這一點(diǎn)。反過(guò)來(lái)也有一類人目前不適合強(qiáng)行擁抱開(kāi)放研究涉及商業(yè)機(jī)密的企業(yè)實(shí)驗(yàn)室涉及個(gè)人隱私的醫(yī)療數(shù)據(jù)研究以及涉及國(guó)家安全等敏感領(lǐng)域的研究。這類研究的開(kāi)放邊界需要非常謹(jǐn)慎不是簡(jiǎn)單把所有文件丟到GitHub就完事。后文我會(huì)專門講一講在合規(guī)框架下怎樣做“有限度開(kāi)放”。2. 搭建一套OpenResearch工作流的頂層設(shè)計(jì)說(shuō)了這么多理念接下來(lái)進(jìn)入實(shí)操層面。別急著下載工具先花十分鐘把整體架構(gòu)想清楚不然很容易裝了一堆軟件卻不知道怎么配合。我給自己的團(tuán)隊(duì)設(shè)計(jì)過(guò)一套工作流核心原則是每一個(gè)環(huán)節(jié)都產(chǎn)出公開(kāi)可追蹤的“證據(jù)”。具體來(lái)說(shuō)整個(gè)科研生命周期被我切成五個(gè)階段每個(gè)階段都有對(duì)應(yīng)的開(kāi)放工具和交付物。第一階段是選題與預(yù)注冊(cè)。研究問(wèn)題確定后先寫一份預(yù)注冊(cè)文檔說(shuō)明研究假設(shè)、樣本量計(jì)算、主要和次要結(jié)局指標(biāo)、分析方法。交付物是一份帶有時(shí)間戳的版本化文檔可以放到OSFOpen Science Framework或自己搭建的Git倉(cāng)庫(kù)里。第二階段是數(shù)據(jù)采集與清洗。這個(gè)階段的關(guān)鍵是原始數(shù)據(jù)絕對(duì)不能動(dòng)所有的清洗和轉(zhuǎn)換操作要用腳本完成。交付物是原始數(shù)據(jù)目錄只讀屬性和清洗腳本數(shù)據(jù)文件建議用CSV或Parquet等開(kāi)放格式避免用SPSS的.sav或MATLAB的.mat這種專有格式替未來(lái)的自己和合作者省點(diǎn)力氣。第三階段是分析與建模。每一個(gè)分析步驟都要有可重跑的腳本或Notebook對(duì)關(guān)鍵結(jié)果保留運(yùn)行日志。交付物是Jupyter Notebook或R Markdown文件、依賴環(huán)境說(shuō)明、隨機(jī)種子記錄。第四階段是論文寫作。用可版本控制的格式寫作Markdown或LaTeX都比Word合適配合Git能看清每一處修改是誰(shuí)在什么時(shí)候做的審稿意見(jiàn)也能逐條對(duì)應(yīng)到修改記錄。第五階段是發(fā)布與歸檔。論文預(yù)印本放到arXiv、bioRxiv或SocArXiv數(shù)據(jù)和代碼歸檔到Zenodo獲取DOI必要時(shí)把整個(gè)軟件環(huán)境打包成Docker鏡像。這一步做完任何人理論上都能從零開(kāi)始把你的研究完整復(fù)現(xiàn)一遍。聽(tīng)起來(lái)復(fù)雜但真正落地時(shí)你只需要把精力集中在少數(shù)幾個(gè)“樞紐工具”上。下面我把每一個(gè)環(huán)節(jié)的工具選型和配置方式拆開(kāi)講。2.1 文獻(xiàn)管理Zotero依然是我最推薦的開(kāi)源方案文獻(xiàn)管理是開(kāi)放研究工作流里最容易被低估的環(huán)節(jié)。很多人用的是EndNote或NoteExpress但它們?cè)凇伴_(kāi)放”這件事上有個(gè)天然短板文獻(xiàn)庫(kù)文件是閉源格式換軟件等于換一套生態(tài)想在團(tuán)隊(duì)里協(xié)作同步也很麻煩。Zotero是我從研究生階段一直用到現(xiàn)在的工具大概用了七八年期間試過(guò)Mendeley、ReadCube Papers和EndNote最終都換回了Zotero。原因很簡(jiǎn)單第一它完全開(kāi)源免費(fèi)數(shù)據(jù)存在本地SQLite數(shù)據(jù)庫(kù)里任何軟件都能直接讀取不存在“綁架”問(wèn)題第二插件生態(tài)極其豐富能自動(dòng)抓取網(wǎng)頁(yè)元數(shù)據(jù)、和Overleaf聯(lián)動(dòng)引用文獻(xiàn)、把文獻(xiàn)庫(kù)同步到WebDAV或自建服務(wù)器第三它在處理中文文獻(xiàn)時(shí)細(xì)節(jié)做得好CNKI條目也能正確抓取作者和期刊信息。安裝Zotero之后的第一個(gè)動(dòng)作我建議設(shè)置“附件存儲(chǔ)目錄”和“同步策略”。打開(kāi)設(shè)置面板把數(shù)據(jù)目錄放到一個(gè)堅(jiān)果云或WebDAV同步的文件夾里這樣換電腦時(shí)文獻(xiàn)庫(kù)能一鍵恢復(fù)。第二個(gè)動(dòng)作是安裝兩個(gè)必裝插件Better BibTeX用來(lái)生成穩(wěn)定的引用鍵ZotFile用來(lái)把PDF附件重命名并移動(dòng)到統(tǒng)一目錄。如果說(shuō)Zotero是一座文獻(xiàn)圖書館這兩個(gè)插件就是圖書館的編目系統(tǒng)和文件傳輸帶缺一個(gè)都會(huì)感覺(jué)別扭。最進(jìn)階的用法是把Zotero的文獻(xiàn)庫(kù)和Quarto或LaTeX寫作鏈打通。Better BibTeX插件會(huì)自動(dòng)生成一份.bib文件你在Markdown里像寫代碼一樣引用文獻(xiàn)時(shí)引用鍵直接匹配Zotero條目編譯生成的參考文獻(xiàn)表一字不差。整個(gè)流程從“查文獻(xiàn)”到“寫論文”不需要手動(dòng)復(fù)制粘貼一條引用信息這種體驗(yàn)一旦習(xí)慣就再也回不到Word手輸引用的時(shí)代了。2.2 數(shù)據(jù)管理用DataLad或DVC管理數(shù)據(jù)和代碼的版本科研數(shù)據(jù)和軟件代碼有一個(gè)本質(zhì)區(qū)別代碼是文本文件用Git管理天經(jīng)地義數(shù)據(jù)往往是大文件可能有幾千個(gè)GB直接塞進(jìn)Git倉(cāng)庫(kù)會(huì)把倉(cāng)庫(kù)撐爆。所以數(shù)據(jù)管理需要專門的工具我試過(guò)的方案里有三個(gè)值得認(rèn)真考慮DataLad、DVC和Git LFS。Git LFSLarge File Storage)是GitHub官方推出的擴(kuò)展方案思路是讓Git在提交大文件時(shí)只記錄文件指針實(shí)際內(nèi)容上傳到LFS服務(wù)器。它在單個(gè)文件不超過(guò)幾百M(fèi)B時(shí)很好用但每月的流量有配額限制免費(fèi)額度用完就得付費(fèi)對(duì)動(dòng)輒幾十GB的視頻或影像數(shù)據(jù)并不友好。DVCData Version Control)本質(zhì)上是“Git的副駕駛”。它不會(huì)替代Git而是把數(shù)據(jù)文件路徑記錄在Git倉(cāng)庫(kù)里數(shù)據(jù)本體存放在本地或云端的存儲(chǔ)桶里。當(dāng)你切換Git分支時(shí)DVC會(huì)根據(jù)記錄自動(dòng)檢出對(duì)應(yīng)的數(shù)據(jù)版本。對(duì)機(jī)器學(xué)習(xí)項(xiàng)目特別友好模型文件、數(shù)據(jù)集、訓(xùn)練參數(shù)都能精確追溯。DataLad是我目前最偏愛(ài)的選擇因?yàn)樗菍iT為神經(jīng)科學(xué)和心理學(xué)研究設(shè)計(jì)的原生支持“數(shù)據(jù)即Git子模塊”的概念。它不僅能版本化數(shù)據(jù)還能從遠(yuǎn)程數(shù)據(jù)集獲取數(shù)據(jù)子集適合大型公開(kāi)數(shù)據(jù)集的多團(tuán)隊(duì)協(xié)作。它的學(xué)習(xí)曲線比DVC陡峭一些但一旦理解“數(shù)據(jù)集是可以嵌套的”這個(gè)思想能力上限會(huì)高很多。走一條最穩(wěn)妥的入門路徑我建議從DVC入手。安裝很簡(jiǎn)單用pip安裝然后在項(xiàng)目根目錄初始化把原始數(shù)據(jù)目錄加入管理每次改動(dòng)數(shù)據(jù)后跑一遍dvc add和dvc push。我自己的習(xí)慣是把DVC和Git配合使用代碼走Git數(shù)據(jù)走DVC兩者靠.dvc文件關(guān)聯(lián)。這樣做的好處是別人克隆你的Git倉(cāng)庫(kù)后只需要一條dvc pull命令就能把你全部的數(shù)據(jù)資產(chǎn)拉到本地整個(gè)復(fù)現(xiàn)過(guò)程非常順暢。2.3 電子實(shí)驗(yàn)記錄讓實(shí)驗(yàn)過(guò)程像代碼一樣可追溯實(shí)驗(yàn)記錄本在傳統(tǒng)科研里是紙質(zhì)手寫本重要實(shí)驗(yàn)要蓋章、簽字、存檔。但在計(jì)算驅(qū)動(dòng)的研究里紙質(zhì)實(shí)驗(yàn)本根本記錄不了腳本的一步步演化過(guò)程。所以現(xiàn)在越來(lái)越多的團(tuán)隊(duì)開(kāi)始用電子實(shí)驗(yàn)記錄本這里要區(qū)分兩個(gè)層面ELNElectronic Lab Notebook和計(jì)算工作流版本控制。正式的ELN產(chǎn)品有開(kāi)源方案比如開(kāi)源界的eLabFTW和RSpace。它們提供帶時(shí)間戳的記錄頁(yè)面、實(shí)驗(yàn)?zāi)0?、?quán)限管理、電子簽名適合化學(xué)、生物這類濕實(shí)驗(yàn)為主的實(shí)驗(yàn)室。我所在的領(lǐng)域是計(jì)算方向的更常用的其實(shí)是“Markdown Git”組合每天新建一個(gè)Markdown文件記錄當(dāng)天做了什么實(shí)驗(yàn)、跑了什么參數(shù)、得到什么結(jié)果、遇見(jiàn)了什么報(bào)錯(cuò)文件名按日期命名每周提交一次到Git倉(cāng)庫(kù)。很多程序員朋友聽(tīng)到這里可能會(huì)笑這不就是寫開(kāi)發(fā)日志嗎確實(shí)本質(zhì)上是同一回事。但科研和開(kāi)發(fā)有個(gè)重要區(qū)別科研的記錄必須保留“失敗分支”。做軟件開(kāi)發(fā)時(shí)你通常只關(guān)注最終能跑的代碼路徑做科研時(shí)那些“跑不通的參數(shù)組合”往往是下一個(gè)正確決策的重要依據(jù)。所以我給自己定了一條紀(jì)律無(wú)論實(shí)驗(yàn)成功還是失敗都在當(dāng)天記錄里寫清楚“今天嘗試了什么”“為什么這么嘗試”“結(jié)果是什么”“下一步改什么”。這套方法堅(jiān)持半年后回頭翻記錄你會(huì)發(fā)現(xiàn)以前模棱兩可的決策來(lái)源其實(shí)都清清楚楚寫在某一天的記錄里。2.4 寫作與發(fā)布Quarto是一條連接分析和寫作的捷徑論文寫作是整個(gè)流程的“出口”也是最容易讓開(kāi)放研究卡殼的環(huán)節(jié)。很多研究者的痛點(diǎn)在于分析結(jié)果在Notebook里論文在Word里圖表需要手動(dòng)導(dǎo)出后一張張粘貼一旦數(shù)據(jù)更新所有圖表要重新做一遍。這種工作方式不僅低效還容易出錯(cuò)——貼錯(cuò)圖、版本對(duì)不上是家常便飯。我推薦用Quarto來(lái)打通這個(gè)環(huán)節(jié)。Quarto是RStudio團(tuán)隊(duì)推出的開(kāi)源科學(xué)出版系統(tǒng)支持Python、R、Julia和Observable JS等多種語(yǔ)言你可以用純文本格式寫論文在文中直接嵌入分析代碼塊編譯時(shí)自動(dòng)執(zhí)行代碼并生成圖表最后輸出為HTML、PDF或Word格式。簡(jiǎn)單說(shuō)論文里每一個(gè)數(shù)字、每一張圖都是從數(shù)據(jù)實(shí)時(shí)計(jì)算出來(lái)的不存在“貼錯(cuò)版本”的問(wèn)題。對(duì)于還沒(méi)用過(guò)Quarto的朋友我提供一個(gè)兩周上手路徑。第一周先別碰論文拿一篇自己的舊文章練手把文字內(nèi)容轉(zhuǎn)成Markdown格式圖表用Quarto的代碼塊重畫一遍不追求完美只求流程走得通。第二周選一篇計(jì)劃投稿的新論文嘗試用Quarto從零開(kāi)始寫重點(diǎn)練習(xí)“參數(shù)化報(bào)告”的寫法把樣本選擇標(biāo)準(zhǔn)、統(tǒng)計(jì)閾值、文件名這些要素定義成YAML參數(shù)這樣同一份Quarto文檔改一下參數(shù)就能生成不同分析條件下的不同版本。這套流程熟練后你會(huì)感覺(jué)寫論文像在寫一份會(huì)自我更新的分析報(bào)告效率提升非常明顯。發(fā)布環(huán)節(jié)國(guó)內(nèi)作者優(yōu)先考慮兩條路徑一是把預(yù)印本放到arXiv、bioRxiv等平臺(tái)上獲得搶跑的公開(kāi)時(shí)間戳二是把最終數(shù)據(jù)和代碼歸檔到Zenodo獲取DOI后寫入論文的Data Availability Statement。如果論文是雙盲評(píng)審記得預(yù)印本發(fā)布和投稿時(shí)間要錯(cuò)開(kāi)具體策略我放在常見(jiàn)問(wèn)題章節(jié)里詳細(xì)說(shuō)。3. 一個(gè)能直接抄作業(yè)的實(shí)操案例水文觀測(cè)數(shù)據(jù)的開(kāi)放研究全流程理論講太多沒(méi)用我拿一個(gè)簡(jiǎn)化但完整的具體案例把上面提到的工具串起來(lái)跑一遍。案例背景是某區(qū)域水文觀測(cè)站積累了10年的日尺度降雨和徑流數(shù)據(jù)我想研究“極端降雨事件對(duì)徑流響應(yīng)時(shí)間的影響”并把這個(gè)研究做成一個(gè)任何人都能復(fù)現(xiàn)的開(kāi)放項(xiàng)目。3.1 第一步建立項(xiàng)目目錄結(jié)構(gòu)目錄結(jié)構(gòu)是整個(gè)流程的骨架我的習(xí)慣是一開(kāi)始就按功能分好目錄避免后期一團(tuán)亂麻。下面是一個(gè)經(jīng)過(guò)多次迭代后我比較滿意的目錄結(jié)構(gòu)extreme-rainfall-runoff/ ├── README.md ├── LICENSE ├── data/ │ ├── raw/ # 原始數(shù)據(jù)只讀權(quán)限 │ ├── processed/ # 清洗后的數(shù)據(jù) │ └── metadata/ ├── code/ │ ├── 01_download_data.py │ ├── 02_clean_data.py │ ├── 03_analysis.py │ └── 04_make_figures.py ├── notebooks/ │ ├── 01_exploration.ipynb │ └── 02_model_diagnostics.ipynb ├── docs/ │ ├── protocol.md # 預(yù)注冊(cè)文檔 │ └── data_dictionary.md # 數(shù)據(jù)字典 ├── paper/ │ ├── manuscript.qmd # 主論文 │ ├── references.bib # 文獻(xiàn)庫(kù) └── environment.yml # 環(huán)境依賴允許我多解釋兩句data/raw目錄一旦放入原始數(shù)據(jù)就設(shè)置成只讀整個(gè)分析過(guò)程只在processed目錄里操作這樣能保證原始數(shù)據(jù)不被污染code目錄里的腳本用數(shù)字前綴標(biāo)明執(zhí)行順序別人克隆項(xiàng)目后照著順序跑一遍就能復(fù)現(xiàn)docs/protocol.md是預(yù)注冊(cè)文檔要記錄研究問(wèn)題和分析方法這個(gè)文件的時(shí)間戳就是“我不作弊”的證據(jù)environment.yml是conda環(huán)境導(dǎo)出文件鎖定所有依賴庫(kù)的版本。3.2 第二步用DVC管理數(shù)據(jù)版本我在項(xiàng)目根目錄執(zhí)行以下命令dvc init dvc add data/raw dvc remote add myremote s3://my-bucket/rainfall-runoff-dvc dvc push這里說(shuō)明一下每個(gè)命令的用途。dvc init初始化DVC環(huán)境dvc add data/raw會(huì)把原始數(shù)據(jù)目錄的元信息寫入一個(gè).dvc文件同時(shí)把數(shù)據(jù)內(nèi)容加入DVC的緩存dvc remote add添加一個(gè)遠(yuǎn)程存儲(chǔ)位置可以是S3、阿里云OSS或者本地其他磁盤dvc push把數(shù)據(jù)推送到遠(yuǎn)程。此后每次有新的原始數(shù)據(jù)進(jìn)來(lái)重復(fù)一遍dvc add和dvc pushGit提交時(shí)連同.dvc文件一起提交。如果哪天發(fā)現(xiàn)某個(gè)分析結(jié)果對(duì)不上你可以用dvc checkout切換到任意歷史版本的數(shù)據(jù)重新跑分析定位問(wèn)題到底出在數(shù)據(jù)更新還是代碼改動(dòng)。這個(gè)能力在傳統(tǒng)工作流里幾乎沒(méi)法實(shí)現(xiàn)。3.3 第三步以提交記錄為時(shí)間軸的實(shí)驗(yàn)管理分析階段我用的是“Notebook用于探索Python腳本用于生產(chǎn)”的模式。探索階段我會(huì)開(kāi)著Jupyter Notebook快速驗(yàn)證各種統(tǒng)計(jì)模型的適用性一旦確認(rèn)某個(gè)分析思路可行就把Notebook里的核心代碼重構(gòu)到code/04_make_figures.py這樣的腳本里保證最終圖表是可重復(fù)生成的。同時(shí)每次跑完一輪完整的分析我會(huì)在Git里打一個(gè)tag命名格式類似analysis_20241013_v1然后在實(shí)驗(yàn)記錄文檔里寫下這個(gè)tag對(duì)應(yīng)什么階段、用了哪套參數(shù)、得到了什么結(jié)論。很多人在這一步會(huì)偷懶覺(jué)得“反正代碼都在Git里歷史記錄不會(huì)丟”。但Git歷史只記錄了代碼變更不能直觀告訴你“哪次提交對(duì)應(yīng)論文的哪個(gè)版本”。打tag加實(shí)驗(yàn)記錄文檔等于給時(shí)間線加上了可讀的注釋三個(gè)月后再看也清楚。3.4 第四步用Quarto生成可復(fù)現(xiàn)論文論文主體用Quarto的qmd格式寫正文是Markdown圖表用代碼塊直接生成--- title: 極端降雨對(duì)徑流響應(yīng)時(shí)間的影響 author: 你的名字 format: pdf bibliography: references.bib --- ## 引言 正文省略 ## 結(jié)果 我們分析了{(lán){ param site_id }}站點(diǎn)近10年的降雨徑流數(shù)據(jù)。 {python} #| label: fig-response #| fig-cap: 降雨事件與徑流響應(yīng)時(shí)間關(guān)系 import pandas as pd import matplotlib.pyplot as plt df pd.read_csv(data/processed/event_response.csv) fig, ax plt.subplots() ax.scatter(df[rainfall_intensity], df[response_time]) ax.set_xlabel(降雨強(qiáng)度 (mm/h)) ax.set_ylabel(徑流響應(yīng)時(shí)間 (h)) plt.savefig(figures/response_time.png, dpi300)編譯時(shí)執(zhí)行quarto renderQuarto會(huì)自動(dòng)運(yùn)行代碼塊、生成圖表、引用文獻(xiàn)、排版成PDF。如果審稿人要求換一種圖的形式比如把散點(diǎn)圖改成箱線圖我只需要改一下代碼塊的繪圖邏輯然后重新編譯所有相關(guān)圖表和描述文字會(huì)自動(dòng)更新不會(huì)出現(xiàn)正文和圖表對(duì)不上的尷尬。 ### 3.5 第五步發(fā)布?xì)w檔 最后一步也是很多人最容易忽略的一步。論文寫完并不代表項(xiàng)目結(jié)束一個(gè)真正開(kāi)放的項(xiàng)目要讓人“接力”得起來(lái)。 我會(huì)做四個(gè)動(dòng)作一是把論文預(yù)印本上傳到arXiv或SocArXiv獲得公開(kāi)時(shí)間戳和預(yù)印本編號(hào)二是把數(shù)據(jù)和代碼打包上傳到Zenodo通過(guò)GitHub和Zenodo的集成自動(dòng)生成版本記錄并分配DOI三是在README里寫清楚“如何復(fù)現(xiàn)”克隆倉(cāng)庫(kù)、創(chuàng)建conda環(huán)境、運(yùn)行特定腳本、編譯論文四是把這個(gè)DOI和倉(cāng)庫(kù)鏈接寫進(jìn)論文的Data Availability Statement里。 關(guān)于第四步期刊投稿時(shí)如果要求雙盲評(píng)審我的策略是先在投稿系統(tǒng)里勾選“數(shù)據(jù)可用但評(píng)審期間保密”把Git倉(cāng)庫(kù)設(shè)為匿名訪問(wèn)或者把代碼托管到像Anonymous GitHub這樣的匿名平臺(tái)等論文正式接收后再公開(kāi)。這個(gè)操作往往比選任何“開(kāi)放數(shù)據(jù)”選項(xiàng)都關(guān)鍵。 ## 4. 開(kāi)放研究實(shí)踐中的常見(jiàn)坑與排查手冊(cè) 方案再好實(shí)操中總會(huì)踩坑。下面這幾個(gè)問(wèn)題是幾乎每一個(gè)剛開(kāi)始做開(kāi)放研究的人都會(huì)遇到的我把它們整理成了一個(gè)速查表也附帶了我自己的排查思路。 | 問(wèn)題現(xiàn)象 | 常見(jiàn)原因 | 排查與解決思路 | | --- | --- | --- | | 數(shù)據(jù)文件太大Git提交時(shí)卡死 | 沒(méi)有用DVC或Git LFS管理大文件 | 把data目錄從Git中移除改用DVC確認(rèn).dvc文件是文本文件 | | Notebook在別人電腦上跑不通 | 依賴庫(kù)版本不一致 | 用conda導(dǎo)出environment.yml并鎖定版本號(hào)用requirements.txt固定pip依賴 | | 圖表的數(shù)字和論文正文對(duì)不上 | 圖表是手動(dòng)粘貼的沒(méi)有實(shí)時(shí)生成 | 改用Quarto或R Markdown嵌入代碼塊實(shí)現(xiàn)“一處改動(dòng)全文更新” | | 審查結(jié)束后數(shù)據(jù)集被刪除 | 沒(méi)有歸檔到長(zhǎng)久的公共平臺(tái) | 上傳到Zenodo或Figshare這些平臺(tái)承諾長(zhǎng)期保存并分配DOI | | 預(yù)印本時(shí)間早于期刊投稿導(dǎo)致期刊拒稿 | 部分期刊不支持提前公開(kāi)預(yù)印本 | 投稿前查一下期刊的預(yù)印本政策雙盲評(píng)審期刊要特別注意 | | 數(shù)據(jù)里包含個(gè)人隱私信息 | 沒(méi)有做去標(biāo)識(shí)化處理 | 刪除姓名、證件號(hào)、精確坐標(biāo)等直接標(biāo)識(shí)地理數(shù)據(jù)做空間模糊化涉及隱私的問(wèn)卷數(shù)據(jù)只發(fā)布聚合統(tǒng)計(jì) | 這里展開(kāi)說(shuō)幾個(gè)最典型的坑。 ### 4.1 坑一數(shù)據(jù)太大Git倉(cāng)庫(kù)直接罷工 我見(jiàn)過(guò)不少朋友用GitHub管理科研項(xiàng)目一開(kāi)始特別順利代碼提交速度飛快。直到某一天他們把一份幾GB的原始影像數(shù)據(jù)塞進(jìn)了Git倉(cāng)庫(kù)然后一切都變了每次git push要等半小時(shí)GitHub提示倉(cāng)庫(kù)超過(guò)建議大小甚至直接拒絕推送。 解決這個(gè)問(wèn)題最干凈的方案是“從一開(kāi)始就分區(qū)管理”代碼走Git數(shù)據(jù)走DVC兩者之間只通過(guò)文本化的.dvc指針文件關(guān)聯(lián)。如果已經(jīng)踩了坑處理起來(lái)會(huì)有一些麻煩需要用git filter-repo把歷史提交中的大文件徹底抹掉然后重新推送倉(cāng)庫(kù)。注意只是刪除當(dāng)前文件是不夠的Git會(huì)保留歷史中的所有版本倉(cāng)庫(kù)體積不會(huì)變小。 ### 4.2 坑二環(huán)境“能跑”但“復(fù)現(xiàn)不了” 科研項(xiàng)目最尷尬的瞬間是審稿人回信說(shuō)“按照你的README我在我電腦上運(yùn)行到第三步就報(bào)錯(cuò)了”。這往往不是代碼寫錯(cuò)而是環(huán)境沒(méi)對(duì)齊Python版本不同、依賴庫(kù)的版本不同、系統(tǒng)庫(kù)缺失都有可能讓“在我機(jī)器上能跑”變成“在你機(jī)器上跑不了”。 最可靠的方案是使用容器技術(shù)把整個(gè)環(huán)境打包。Docker和Singularity是科研領(lǐng)域最常用的兩個(gè)Singularity在HPC集群上更常見(jiàn)因?yàn)樗恍枰猺oot權(quán)限且與集群調(diào)度器集成更好。我一般會(huì)在項(xiàng)目的environment.yml里定義conda環(huán)境然后寫一個(gè)Dockerfile把這個(gè)環(huán)境完整地打包成鏡像。這樣做最大的好處是復(fù)現(xiàn)方只需要裝一個(gè)Docker然后run一條命令環(huán)境就一模一樣地跑起來(lái)了。雖然前期學(xué)習(xí)成本高一點(diǎn)但相比“審稿人跑不通”帶來(lái)的損失這點(diǎn)成本完全值得。 ### 4.3 坑三開(kāi)源協(xié)議沒(méi)選對(duì)數(shù)據(jù)被用了還不能聲張 開(kāi)放研究里的“開(kāi)放”不等于“沒(méi)有權(quán)限約束”。代碼、數(shù)據(jù)、文本這三種資產(chǎn)最好分別選擇不同的許可證。代碼通常用MIT、BSD、Apache 2.0或GPL數(shù)據(jù)通常用CC0或CC BY論文文本通常用CC BY。選許可證不是隨便找個(gè)模板復(fù)制要明確“別人能用你的東西做什么、必須滿足什么條件”。 讓我舉一個(gè)真實(shí)的例子我在一次項(xiàng)目里拿到了某機(jī)構(gòu)提供的水文數(shù)據(jù)對(duì)方在公共門戶上標(biāo)了“開(kāi)放使用”但實(shí)際數(shù)據(jù)包里又夾著一份PDF使用協(xié)議規(guī)定“禁止用于商業(yè)目的”。這就是典型的許可證不清晰。如果不注意這一層貿(mào)然把數(shù)據(jù)整合進(jìn)自己的開(kāi)放項(xiàng)目里后期可能面臨法律糾紛。 在做數(shù)據(jù)發(fā)布前我的自查清單是數(shù)據(jù)來(lái)源有沒(méi)有明確的許可或授權(quán)聲明如果沒(méi)有就只發(fā)布自己采集或明確擁有權(quán)利的數(shù)據(jù)數(shù)據(jù)里有沒(méi)有第三方版權(quán)的內(nèi)容有的話需要先獲得授權(quán)選用的許可協(xié)議是否和期刊、資助方的要求沖突各個(gè)機(jī)構(gòu)對(duì)開(kāi)放許可的要求不完全一樣以資助方的口徑為準(zhǔn)。 ### 4.4 坑四把“公開(kāi)”當(dāng)成“可復(fù)現(xiàn)” 有些項(xiàng)目雖然把所有文件都公開(kāi)了但第三方根本復(fù)現(xiàn)不了。比如公開(kāi)了一份Excel表但沒(méi)說(shuō)明每一列的計(jì)算邏輯公開(kāi)了一份訓(xùn)練好的模型權(quán)重但沒(méi)有提供訓(xùn)練代碼和超參數(shù)公開(kāi)了一個(gè)Notebook但中間狀態(tài)已經(jīng)被手動(dòng)修改過(guò)執(zhí)行順序?qū)Σ簧稀?“公開(kāi)”是手段“可復(fù)現(xiàn)”才是目的。判斷一個(gè)項(xiàng)目是否可復(fù)現(xiàn)我常用一個(gè)“陌生人測(cè)試”讓自己完全不參與這個(gè)項(xiàng)目的同事只憑README和公開(kāi)文件從頭跑一遍完整流程看能不能得到論文里的核心結(jié)果。如果這一步走通了這個(gè)項(xiàng)目才能真正被稱為開(kāi)放研究項(xiàng)目而不是單純把文件掛在網(wǎng)上。 ## 5. 我沉淀下來(lái)的幾條選型和建議 工具選型是一個(gè)動(dòng)態(tài)調(diào)整的過(guò)程沒(méi)有哪一套方案是放之四海而皆準(zhǔn)的。根據(jù)自己的項(xiàng)目類型和團(tuán)隊(duì)背景我對(duì)不同工具有一些親身的體會(huì)和建議。 ### 5.1 關(guān)于“先選工具”還是“先定流程” 我的經(jīng)驗(yàn)是先定流程再選工具。流程是骨架工具是血肉。如果你連“數(shù)據(jù)怎么組織”“分析怎么記錄”“產(chǎn)出怎么發(fā)布”都沒(méi)想清楚裝再多的開(kāi)源工具也只會(huì)讓項(xiàng)目更亂。反過(guò)來(lái)流程一旦確定工具選型就清晰很多需要版本控制就選Git需要管理大文件就選DVC或DataLad需要一體化寫作就選Quarto或R Markdown。 在團(tuán)隊(duì)協(xié)作場(chǎng)景下流程比工具更關(guān)鍵。我見(jiàn)過(guò)兩個(gè)團(tuán)隊(duì)用完全一樣的工具套件一個(gè)協(xié)作順暢一個(gè)天天沖突差別就是前者在項(xiàng)目啟動(dòng)第一天就約定了分支策略、命名規(guī)范和數(shù)據(jù)管理流程而后者完全靠個(gè)人自覺(jué)。所以如果要組隊(duì)做開(kāi)放研究建議第一次組會(huì)就拿出一份書面的“項(xiàng)目操作手冊(cè)”把流程約定下來(lái)。 ### 5.2 關(guān)于開(kāi)放研究的“度”與合規(guī)邊界 我始終強(qiáng)調(diào)開(kāi)放不是目的創(chuàng)造價(jià)值才是。有些數(shù)據(jù)天然不適合完全公開(kāi)比如涉及個(gè)人隱私的醫(yī)療健康數(shù)據(jù)、涉及商業(yè)機(jī)密的行業(yè)調(diào)研數(shù)據(jù)。這類研究的開(kāi)放策略是“有限度開(kāi)放”數(shù)據(jù)本身不公開(kāi)但分析代碼、處理流程、結(jié)果報(bào)告可以公開(kāi)原始數(shù)據(jù)留在受控環(huán)境里需要申請(qǐng)和審批才能獲取。 具體做法上可以參考一些大型生物醫(yī)學(xué)項(xiàng)目的做法數(shù)據(jù)分級(jí)管理一級(jí)數(shù)據(jù)完全公開(kāi)二級(jí)數(shù)據(jù)需要在平臺(tái)注冊(cè)并提交使用申請(qǐng)三級(jí)數(shù)據(jù)則只能在指定的安全計(jì)算環(huán)境內(nèi)分析連下載都不允許。這種模式既保障了開(kāi)放協(xié)作又守住了安全和隱私的底線。 對(duì)于涉及人類被試的研究還有一個(gè)容易被忽視的點(diǎn)知情同意書里有沒(méi)有寫“數(shù)據(jù)將公開(kāi)共享”。如果當(dāng)初收集數(shù)據(jù)時(shí)只告知被試“數(shù)據(jù)僅用于本研究”那后續(xù)是不能直接把原始數(shù)據(jù)公開(kāi)的。這種情況下可以考慮發(fā)布聚合統(tǒng)計(jì)結(jié)果或脫敏后的衍生數(shù)據(jù)并說(shuō)明數(shù)據(jù)共享的邊界。 ### 5.3 關(guān)于長(zhǎng)期維護(hù)與“人走項(xiàng)目涼”的問(wèn)題 學(xué)術(shù)圈人員流動(dòng)快一個(gè)博士生畢業(yè)后他負(fù)責(zé)的項(xiàng)目很可能就沒(méi)人維護(hù)了。開(kāi)放研究項(xiàng)目一旦發(fā)布理論上就屬于公共知識(shí)資產(chǎn)但“發(fā)布”不等于“有人維護(hù)”。如果項(xiàng)目倉(cāng)庫(kù)里的代碼依賴了某個(gè)已經(jīng)停止維護(hù)的庫(kù)或者數(shù)據(jù)更新到一半就斷了后續(xù)使用者的體驗(yàn)會(huì)很快惡化。 我的應(yīng)對(duì)策略有三條一是盡量在論文接收后的一兩個(gè)月內(nèi)把代碼、數(shù)據(jù)、文檔一次性完整歸檔到Zenodo生成快照版本哪怕以后倉(cāng)庫(kù)不再更新快照版的DOI仍然可用這意味著研究結(jié)論至少還能被找到二是盡量用長(zhǎng)期存在的公共平臺(tái)而非個(gè)人服務(wù)器做托管用好GitHub和Zenodo這些老牌平臺(tái)三是在README里明確“維護(hù)狀態(tài)”寫上最后更新時(shí)間和維護(hù)人聯(lián)系方式后人接手時(shí)至少有線索可循。 說(shuō)到底開(kāi)放研究的終極目標(biāo)不是讓每一個(gè)項(xiàng)目都萬(wàn)古長(zhǎng)青而是讓每一個(gè)研究過(guò)程都經(jīng)得起追溯。哪怕人走了證據(jù)還在后來(lái)者可以接續(xù)工作這就是它比封閉科研更有生命力的地方。 最后再分享一個(gè)小體會(huì)搭建這套工作流前期的學(xué)習(xí)曲線確實(shí)會(huì)陡一點(diǎn)尤其是剛開(kāi)始接觸Git和DVC時(shí)會(huì)有那么一兩個(gè)星期覺(jué)得“這比寫論文還累”。但熬過(guò)這個(gè)階段當(dāng)你三個(gè)月后回頭翻自己的Git提交記錄發(fā)現(xiàn)每一個(gè)結(jié)論都能找到對(duì)應(yīng)的代碼和數(shù)據(jù)當(dāng)你把項(xiàng)目鏈接發(fā)給合作者對(duì)方一個(gè)下午就復(fù)現(xiàn)了你的全部結(jié)果——那種踏實(shí)感是傳統(tǒng)科研方式給不了的。我的建議是不用等“準(zhǔn)備好”再開(kāi)始就從下一篇論文甚至從這篇論文的“數(shù)據(jù)分析”階段開(kāi)始把至少一個(gè)環(huán)節(jié)接入開(kāi)放工作流。邁出這一步你就已經(jīng)走在“OpenResearch”這條路上了。
返回列表
PREV
查看更多資訊
NEXT
返回資訊列表
超碰超碰95| 久久久久久久久久久久黄色| 麻豆国产视频精品观看| 国产又粗又长又大的视频| 日韩精品.久久精品.AV女优.天美传媒| 久久激情亚洲精品无码?V| 综合网亚洲1| 91青视频| 日韩另类色图| 另类小说五月天| 秋霞免费无码视频日韩A片| 欧美伦乱爱| 亚洲人妻av| 国产又粗又大硬免费色网视频| 性色av一区二区| 亚洲综合五月天| 国产一区二区a毛片| 亚洲中文字幕av | 五月婷婷爱六月丁香色| av无码av无码专区| 十八禁av无码免费网站APP| 天天综合青苹果| 俺也射| 东京热男人的天堂网| 亚洲 欧美 手机在线观看| 日韩三级伊人| 一级性爱网| 欧美 亚洲 在线| 成人一二三区| 91综合无码| 久久久国产护士丝袜美腿一| 91精品久久久久久77777| 久久久99久9| 97爱免费插| 精品一区二区三区蜜桃臀赵总 | 天美av在线观看| 亚洲综合网91| 中文字幕亚洲永久精品| 99久久久99久久91熟女| 亚洲偷拍自拍在线视频| 久久精品国产亚洲AV片多多 | 国产1024在线播放| 高清国产无码av| 精品免费视频国产一区| 强奸乱伦大香蕉| 青娱乐欧美激情一区二区 | 1024亚洲中文字幕久在线看片你懂的 | 十八禁啪啪视频| 美女裸体无遮挡永久免费观看网站| 久久人妻一区二区三区高清| 国产高清在线自在拍69| 久久久久深夜无码| 97免费视频在线| 中国熟妇| 久久91| 91 欧美| 澳门黄片一香蕉视频| 亚洲国产欧美一区二区潘金莲| 大香蕉免费3| 黄色成人网久久久久久| 激情五月天丁香社区| 欧美日韩国内不卡| 天天天天干| 日韩人妻大香蕉| 蜜臀精品1区2区| 天天干人妻| 欧美BT 亚洲色图| 超碰美女97| 东京热大香蕉| 亚洲欧洲精品视频发布| 国产刺激视频| 黑人精品久久97| 日本熟妇人妻一区二区三区| 久9re热视频这里只有精品| 亚洲图片激情综合另类| 国产乱婷婷精品二区三区| 吉川爱美亚洲二区在线| 97亚洲欧美| 日韩在线性爱免费视频| 色婷婷丁香五月| 97综合在线观看| 四虎视频在线观看| 伊人骚琪琪亚洲天堂网站| A 天堂| 久久精品国产亚洲AV无码电影| 97干在线| 加勒比伊人综合| 97久久超碰国产精品| dy888午夜老子影视达达兔| 九九综合久久中文字幕| 国产熟女少妇一区| 久久香蕉国产线看观看亚洲女人 | 久久精品成人一区二区三区蜜臀| 久久久96| 7月婷婷综合| 人妻大香蕉| 91久久久亚洲| 天天干天天爽| 18禁无码永久免费无限制| 懂色中文一区二区三区| av天天在线观看| 青青草久久| 欧美AAAA黄片| 又大又长又爽| 密臀在线一区尤物| 午夜欧美女人操逼| 久久影视二区三区行押| 黄色电影观看久久9| 五月婷婷综合在线| 色九九久九九| 久久 国产精品 一区| 日韩精品国模| 国产www色在线观看| 亚洲 图片 欧美 色图| 欧美综合91| 国产精品熟女丝袜一区二区| 欧美97在线欧| 亚洲激情综合| 国产精品九九九| 久99久视频精选| 综合情欲网| 91网站18+| 国产东北女人在线视频| 色好看av| 国产精品亚洲高清在线| 日韩9999| 亚洲最大无码中文字幕网站| 九9精品| 国产成人亚洲精品无码古代早漏男| 欧美日日网| 国产99999| 超碰97欧美在线| 一起草视频在线| 欧美99热| 婷婷性网| 免费看美国人人爽,人人操| 一本久久久精品| 欧美国产日韩高清在线| 亚洲不卡一| 亚洲成人美女无吗| 国产乱不卡| 90后性网国产欧美| 欧美色网络| 很很干很很操| 国产亚洲精品无码三区| 欧美综合自拍亚洲综合图| 亚洲日韩成人性爱视频| 久久精品导航| 国产精品分类在线观看| AAAA欧美日韩| 艹精品| 视频二区熟女人妻| 超碰人人乐97| 欧美精品成人一区二区在线观看| 尤物视频偷拍免费| 长长久久88视频| 中文久久| 人妻 欧美亚洲| 99无码视频| 女人综合网| 宗合情欲网| 美中日韩无码| 精品78| 另类图片亚洲加勒比另类图片亚洲加勒比另类图片亚洲加勒比 | 日韩激情中文字幕有码| 少妇高潮一区二区三区在线| 97视频网站在线观看| 岛国黄片网站| 久日综合网| 婷婷在线播放| 操逼大黄片| 久久久久久99999国产精品| 综合激情五月丁香| 精品九九国产无码| 91殴美| 极品粉嫩少妇视频| 男人的天堂2018.| 欧美日韩岛国大片在线观看| 国产偷人伦激情在线观看| 在线小说视频一区| 欧美性色综合网| 日韩久久激情精品| 黑丝自慰喷水网站| 成年人网站在线免费观看| 青青操日韩| 成人久久精品| 96爱综合| 97色97好| 97人人夜| 日韩午夜国产| 97超碰天天爱天天爱| 啪啪啪综合网| 激情五月综合网| 日美免费黄片| 国产AV激情无码久久无码| 国产在线综合福利网站| 久久精品国产欧美日韩亚洲欧美日韩中文久久国产一区 | 亚洲色婷婷久久久综合日本| 久久禁| 久湿久久| 丝袜天堂网| 粉嫩AV一区夜夜嗨| 蜜臀人妻少妇久久在线观看| 色噜噜人妻av 中文字幕| 91草草草| 做爱福利视频一区二区| 97亚洲欧美日韩| 蜜桃久久久久久久| 国产性感骚丝袜在线| 亚洲国产第一页综合视频| 国产热RE99久久6国产精品首 | 精品视频97| 中文字幕一区 二区三四五 区日 日骚| 欧美成人一级免费电影| 91三级理论片播放器| 日本狂喷奶水在线播放212| 日韩免费福利在线观看| 欧美亚洲第一页| 啪啪视频亚洲第一| 黑操B| 婷婷激情四射| 久久久久久久 九九九九九九九| 成人在线视频二区| 99青青草国产视频| 亚洲网污污污污| 欧美亚洲高清不卡| 2019亚洲男人天堂| 无码人妻精品酒店| 9久在线视频只有精品| 骚货操死你| 偷偷人人精品女女久久| 五月天亚洲色图| 五月天亚洲色图| 97超碰站| 欧美日韩大黄片| 激情黄色片在线观看| 26uuu最新| 亚洲人妻日日日| 中文乱码字幕观看| 一二三四视频中文字幕在线看| 亚洲色图亚洲| 亚洲图片婷婷五月天| 国内毛片婷婷六月色| 人人看人人摸人人色| AV九九| 校园春色综合| 风骚少妇视频中文字幕| 99久久久| 99国产精品免费| 九九国产| 99热伊人| 丝袜剧情| 久久久com| 美女主播色欲91抠b在线播放| 污污污8888| 亚洲日精品| 午夜性刺激视频免费观看| 天天干天天狼在线视频| 婷婷尹人大香蕉免费| 素人美腿视频网站| 欧美女同在线| 欧美传媒| 青娱乐亚洲热| 国产嫩草精品A88AV| 97日韩欧美| 日本性爱不卡视频| 91丝袜美腿网站| 亚洲精品97久久中文字幕| 狠狠爱综合网| 高清国产av无码| 9久久9综合| 草B在线| 2023天天操夜夜操| 97操97色| 一二三区视频在线观看| 欧美性生活免费网| WWW4虎| 91国产丝袜美女| 亚洲h片在线免费观看| 亚洲高清在线| 91 国产丝袜在线播放-百度| 久久精品午夜国产亚洲AV无码| 亚洲精品日韩国产欧美| 玖玖人人爱| 久久久精品九| 曰本人妻人人澡人人夹| 97视频620| 免费操逼91| 操美女高潮抽搐白浆| 蜜臀99999| 超碰免费人妻在线| 中文字暮97| 国产精品伦理| 99性爱| 日本天堂在线播放| 国产丝袜欧美在线视频| 91麻豆天美传媒HD| 欧美性爱第1 页| 亚洲天堂区| 精久久久| A 天堂| 一色网男人的天堂| 国产精品一区二区a| 五月丁香影院| 色原狠狠天天天| 美国精品国产精品| 日日夜夜青青草母狗| 日韩精品中文字幕人妻| 久久精品国产99国产精品亚洲| 一区二区中文| 色臀aV| 亚洲精品无码成人久久久99| 人人操人人精品影片| 日本一区二区三区四区免费观看| 亚洲情色第一页| 91亚洲综合| 国模限制级电影| 欧美AB在线| 蜜臀在线免费观看在线免费观看| 三男一女不戴套的A片| 中文字幕二区日韩天堂| 精品一级毛片在线观看| 精品99999久久久久久| 性爱视频免费网址| yy少妇精品久久| 日日噜噜夜夜狠狠视频无| 亚洲少妇在线影音| 99热| 成人一道本免费视频| 婷婷激情一区二区三区俺也去| 日韩黄色片子| 久久女女| 国产性感在线观看| 中文字幕一区日韩精| 上床啊啊啊| 欧美黄色图片| AV天堂男人的天堂| 日本肏逼视频在线观看| 91在线视频国产网站| 东北女人的毛片| 欧美婷婷五月天| 狠狠躁伊人中文字幕| 91大胆欧美| 操逼片国产| 99操逼| 激情文学欧美| 九九热精品| 宅男影院久久久,99| 99最新日韩偷拍视频| 久久精品视-一级做a爰片性色毛片16美国-中国女与老外在线精品 | 日本十八禁免费看污网站| 97在线欧| 人人超碰在线观看黄| 婷婷亚洲五月***久久| 起碰97| 色男人色天堂东京热| av中文在线| 激情抓乳插进去啪啪啪日韩| 2020中文字幕在线观看| 呦女网站| 九九九精品一区二区无码| 亚洲综合小说另类图欧美视频激情小说色五月天| 天天看少妇| 91中文字幕制服丝袜免费视频| 久久男人精品| 精品国产三级av韩国在线| m欧洲一级午老| 亚洲熟妇无码一区二区三区| 中文字幕熟女人妻丝袜| 超碰人妻97| 男人的天堂日韩| 91粉芽高清在线一区二区| 一起草日韩| aaa淫乱视频| 凸凹视频在线观看| 中文字幕亚洲热播人妻| 欧美草草高清日韩视频| 成全动漫视频观看免费下载| 日本成熟少妇A∨网站| 国产一级舔足在线观看| 综合色欧美| 蜜臀久久精品久久久久视频| 亚洲丝袜B诱惑| 粉嫩国产精品久久粉嫩| 蜜臀少妇一区二区| 日本有码久久| 91人人爽人人爽人人人,gav福利视频导航,日韩欧美亚洲国产字幕四区 | 亚洲男人天堂网久久| 亚洲成人av色网| 欧美激情综合| 欧美亚洲日韩人妻在线观看| 精品人妻美妇91job| 色哟哟 日韩精品| 超碰午夜| 日韩有码 一区二区三区| 人妻丝袜日本| 亚洲色图91| 日韩日本欧美在线观看| 99精品久久久久久久婷婷| 久久中文字幕不卡人妻| 免费少妇一区二区| 精品无码一区二区三区| 最新岛国大片| 在线可观看的黄色网址| 精品性爱无码在线播放| www.色婷婷| 99热免费| 天天看特黄的免费网站| 国产亚洲日本| 东京男人天堂| 欧美精品成人一区二区在线观看 | 精品人妻av在线播放| av天堂精品久久| 隔壁邻居波多野结衣中文字幕| 亚洲最新中文字幕免费| 久久XX| 狠狠躁AV| 国产十八禁视频| 噜噜噜亚洲精| 亚欧无码线免费观看视频| 久久久精品中文字幕爱豆| 久久精品国产AV一区二区三区| 国产福利在线视频网站| 亚洲18禁| 97久久国产亚洲精品超碰热| 欧美AB在线| 岛国片在线视频网站| 91操人| 大香蕉久久| 日韩精品人妻中文字幕有码午| 综合天天网| 人妖欧美一区二区| 久草尤物| 亚洲操逼无码| 91在线免费观看处女| 人人妻人人色| 伊人9| 亚洲成人精品在线一区| 入口操逼网站| 观看免费区二区三区二| 亚洲天天精品| 日韩欧美成人综合在线| 女优视频第10页| 婷婷五月天无码| 天堂网亚洲区手机版| 国产亚洲日韩在线三区黑人| 久久婷婷国产一区二区色| 青青草乱入乱欲视频在线观看| 大茄子熟女AV导航| 人人操人人大香蕉| 和协影院中文字幕三区| 91精片| 亚洲性爱无码乱伦av| 97超碰亚洲| 天美91| 人人扣人人操| 后入福利视频| 91精品丝袜在线观看| 五月天久久婷婷亚洲| 欧中美三级一区二区三区| 搞中出久久| 丰满的三级少妇欧美久久久| 亚洲做性| 中文字幕av亚洲在线| 亚洲色性情三级| 欧美色图片欧美色图| 97在线观看播放视频| 国产精品乱人伊人网| 亚洲男人在线观看天堂| 午夜激情成人在线观看| 热久日综合| 国产精品午夜福利亚洲综合网| 日韩人妻精品久久久久| 18精品一二区| 啊啊啊啊嗯嗯嗯用力好爽 | 91老司机精品| 综合一区中亚洲国产成人综合精品| 青青操视频在线| 97在线视频观看网站| 丝袜视频一区二区在线播放国产中文 | 无码一区二区三区四区五区六区七区八区九区十区视频 | 综合欧美亚洲| 四虎影库国产精品免费| 精品999999| 天美av在线观看| 国产成人久久精品蜜臀| 青青草视频久久久久| a片自拍直播视频| 一区二区久久天天干狠狠| 国产精品噜噜噜日日日| 中文字幕av久久爽Av| 亚洲成人精品在线一区| 97视频620| 亚洲色 国产 欧美 日韩| 蜜臀久久99精品久久久久电影| 国语对白在线播放视频| 九九九九九九综合| 天天干少妇| 一个人免费视频观看在线WWW| 情色五月天久久久| 中文字幕日韩人妻视频一区二区三区| 日本天天干天天日一区| 国产蜜臀精品一区免费尤物| 精品999日本| 1769精品一区二区三区| 欧美精品精品一区二区| 国产白丝精品在线观看| 久久久成人免费av电影| 日本天天色| 日韩精品中文字幕人妻| av久日| 国产在线精品偷| 男人天堂2030| 日日干夜夜欢| 按摩中文字幕| 色色香蕉| 91网站18+| 9久综合网| ji熟女.com| 天天碰久久入| 狠狠亚洲| 亚洲色吧网| 搡老女人老熟女91| 人妻少妇三级| 国内91熟女人妻丝袜天天精品视频在线 | 蜜桃臀一区二区三区久久| 日韩成人精品中文字幕| 婷婷去俺也去六月色| 人人妻人人澡人人爽久久av| 精品无码久久久久| 91neishe| 99久久99九九99九九九| 久久五月天婷婷丁香中文字幕| 噜噜噜噜天天狠狠| 日本最新免费韩国1区2区视频播放| 久久久久婷婷| 多毛小伙内射老太婆| 91综合网站| 熟妇激情| 日韩精品9区| 亚洲国产综合久久久性感熟妇| 久久精品国产亚洲AV片多多| 人人摸.人人色| 五十路熟女人妻一区二区在线观看| 老女人老91妇女老热女| 欧美日韩222| 久久久艹艹艹| 激情综合色| 国产吞精a级片激情电影| 黄片免费视频2019| 亚洲一区日韩精品中文字幕| 夜夜嗨一区二区三区直播内容| 黑人综合网| 情色图区| 拍拍拍拍大尺度黄色三级片拍拍拍拍拍照| 欧美一二在线| 97国伦国色| 色第一页| 欧美性暴力猛交XXXX| 成片免费播放| 99操逼| 国产精品老师| 99视频内射三四| 丁香五月天堂| 爱射综合| 精品久热| 禁片 高清 在线观看视频网站| 99热在线播放| 婷婷五月天小说| 啪啪啪综合网| 九九热免费视频| 日韩草久视频| 亚洲AV无码天美传媒一区| 五月天婷婷久久| 一级免费啪啪片| 成人短视频在线观看| 91久操| 99热伊人| 欧洲亚洲人妻无码中字久久三区四区| 久插综合| 男人的亚洲天堂| 91熟女视频网| 好吊色综合| 日本东京热久久久电影| 日韩人妻一区二区| 日韩欧美tv一区二区在线观看| 日韩久久.一级黄色片| 亚洲精品官网在线观看| JULIA一区二区三区在线播放| 淫荡少妇免费| 91色婷婷综合久久中文字幕二区| 人妻丝袜无 码视频专区| 91蜜臀熟女| 久久 亚洲 日韩 人妻| 蘋果手機免費看成人Av| 小电影欧美91| 牛牛操视频逼| 成人性爱高清视频免费看| 26uuu国产免费观看| 国产CHASE男男GAYGA 毛多色婷婷| www.zbzhongsen.com| 国产超碰国产97| 天天综合网在线| 色97欧美| 不卡中文字幕aⅴ在线| 老熟女乱子伦中文字幕一区二区| 久久性爱城| 97玖玖人妻| 无码操逼天堂| 97综合久久| 人人爱人人操人人性| 天天日少妇逼AV| 亚洲国产精品久久久男人的天堂| 1024精品在线| 中国一级操逼视频| 99久久婷婷国产综合精品草原| 一区二区三区蜜桃成人撸久久东京热 | 99热超碰| 天堂亚洲精品| WWW操逼| 懂色AV一区二区三区| 2017大香蕉国产精品久久| 亚洲色图A| 国产精品久久久无码AV网站| 亚洲情色婷婷五月天| 欧美黑人极品高潮喷吹熟女黑人性暴力日韩在线欧美极品一区二区老师 | 精品成人无码| 91在线色综合| 亚洲精品久久久久久| 人人摸人人入| 殴美,日韩国产伦精品| 97超碰这里只有精品| www久久国产精品| 亚洲欧美精品91| 欧美亚洲清纯| 日韩性爱1级片视频| 欧美熟妇乱码在线一区| 久久久青草青青国产亚洲免观精品高清完整版_97久久综合区小说区图片区,国精品 | 少妇贴图| aⅴ日韩成人电影av在线免费看av大全| 亚洲男人天堂Av| 91天堂丝袜美腿| 无码聚合| 中文字幕第95页| 岛国毛片在线观看免费| 日韩美女高潮喷水视频| 先锋影音av先锋一区| 6080yy午夜理论三级一区二区三区无码| 久久精品噜噜噜成人看免欧美大片| 99这里只有精品国产| 色色热| 久操免费电影| 99热导航| 大学生美女口爆| 麻豆视频test| 亚洲综合影院| AV电影在线播放| 欧美亚洲特P| 超碰av在线| 亚洲图片欧美偷拍| 爱射综合| 偷拍 欧美 日韩| 亚洲官网在线| 中文字幕一区二区日韩网| 亚洲欧美另类少妇精品| 国内偷自视频区视频综合| 蜜桃视频一区二区三区在线观看| 日韩性爱小视频| 日本天天干天天搞一区| 精品亚洲国产成人精品| 精品美女人人干| 91狠狠色丁香婷婷综合久久精品| 久久久啊啊啊| 国产欧美精品日韩区二区麻豆天美| 人人模人人看| 日韩99神马视频片| 人人操人人干网页| 色偷偷2020免费视频播放| av绯色| 久久宗合97| 亚洲最大AV网| 亚洲欲色| 亚洲狠狠入| 男人天堂无码| 2017,超碰| 天干天干天干天天做| 国产精彩女在线观看视频| 欧美色一二三| 磁力99AV| 亚洲综合九| 欧美色综合图片| 亚洲蜜臀视频精品久久| 台湾佬中文娱乐自偷自拍| 91女人的网站| 国产精品久久9| 婷婷丁香一区二区三区| 欧美日韩少妇色情| 欧美亚洲成人在线一区二区三区| 91色s| 330dv亚洲成年视频网| 超碰1997| 亚洲a色| 91爱啪| 色综合加勒比| 一级性爱啪啪视频| 极品尤物自安慰| 亚洲 欧美 手机在线观看| 久久综合激情| 亚欧精品久久久久久久久久久| 久久欧美按摩999| 97 国产一区| 超碰天天操你比| 亚洲人综合| 蜜臀99久久国产| 97日韩欧美| 久综合国内精品自在自线| 欧美东京热精品A∨| 日本国产成人亚洲精品无码| 色五月亚洲| 91夜夜蜜桃臀1区2区3区| 日本精品久久久久久久| 超碰97久| 日韩乱伦AⅤ| 大香蕉一线视频| 一区二区偷拍拍视频| 蜜臀久久久国产| 日韩BBN| 亚洲天天操| 日韩国产乱子伦App| 五月天色色色| 欧美亚性天堂| 成熟熟女国产精品一区二区| 69麻豆天美| 强奸乱伦亚洲第一页| 99热最新| 欧美日韩亚洲五月天婷婷| 久久在肏| 日韩国产不卡在线视频| 殴美牲| 白丝在线一区| 亚洲综合有码| 视频国产精品未满十八禁止在线观看| 天天射,天天操,天天爽-国内精品一区二区三区-成人AV | 老师充足的奶水小说| 国产91丝袜在线播放蜜月| 男女性无套 免费九一| 人干人人人操人人摸| 五月天激情小说| 欧美日韩亚洲国产中文永久天天看 | 婷婷久久综合| 亚洲AV不卡在线观看| 日韩肏逼视频| 久久精品美女一区| 日日夜夜模| 视频不卡中文字幕| 丁香激情网| 人妻内射一区二区在线视频| 久久久人妻| 99热日本| 国偷自 一区二区| 中文字幕精品免费一区二区| 91在线欧色| 欧美偷偷网| 免费福利视频中文字幕| 美女啊啊啊啊啊啊啊| 色综合 加勒比| 久久av成人无码免费| 欧美大香蕉在线观看| 国产在线激情| 色婷婷导航| 97就爱干| 四虎国产精品永久地址入口| 五月天综合网| 日韩性爱一级片| 亚洲九月丁香| 欧美在线第五页| 日韩国产乱子伦App| 裸体1区| 免费1级a做爰片观看| 性欧美体内射精| 国产狂喷潮在线精品| 亚洲91少妇| 中文乱码字幕观看视频| 色婷婷基地| 久久婷婷一区| 九九热精品| 鸥美中出| 亚洲久热| 综合第一页| 操少妇很爽av| 3PAV乱伦视频| 男人的天堂不卡一区二区| 天堂亚洲欧美| 国产一级黄色片在线观看| 青青草日韩免费观看高清在线| 啪啪啪东京| 欧美97日韩| 国产精品盗摄 偷窥盗摄| 禁十八久久| 日韩亚洲中文字幕在线| 午夜后入| 久精品无码av一区二免费国产在线观看| 国产中午字一暮区| 亚洲国产欧美中文永久| 欧美懂色综合网| 91色艳| 久久中文字幕不卡人妻| 美女诱惑在线一区| 高跟丝袜AV专区国产| 综合五月天| 长长久久免费视频| 成人国产二区三区在线,男女精品。| 骚货人妻偷情自拍在线视频| 五月亭亭六月丁香| 精品无人区麻豆乱码久久久| 少妇一区二区三区精选| 99福利社| AV无码久久久精品| 亚洲高潮影院| 91国产丝袜足交精品视频| 亚洲成人性| 日本男人插女人的逼黄色| 日本欧美亚洲高清在线看| 97天天在线| 五月综合色| 人人妻碰人人免费| 午夜视频久久久| 日本韩国一本产品小视频日本韩国一本产品久久久产品小视频日本韩国一本产品久 | 五月婷网站| 91丝袜人妻| 泰国AV在线观看| 国产精品成人蜜臀AV在线| 狼人狠干| 欧美日韩国产人人| 亚欧精品久久久久久久久久久| 日本伦乱九九九综合| 加勒比久久av| 久久9精品视频| 97啪啪| 岛国1区2区3区在线观看| 久操影视| 一区二区精品更新提醒| 欧美在线观看综合国产| 91视频伊人| 美女好片色日本| 日韩精彩免费| 五月丁香综合网| 18禁免费视频| 国产兽交视频在线播放| 91熟女视频网| 最新精品久久蜜桃 | 成人免费福利网站国产| A片大香蕉在线| 日日爽熟女| 性性久久| 黑丝日韩av丝袜av| 丁香五月婷婷五月| 极品色www影院| 人人看人人插| 日韩视频中文字幕| 欧美色老汉| 和协影院中文字幕三区| 黑丝91视频| 99久热| 熟女突然公开看18禁影片| 少妇一级无码精品| 国产乱青青草久久| 天天看片麻豆| 免費人妻夜夜爽天天爽爽一区| 天天日夜夜| 国产视频一区二区免费| 日本一区二区三区四区免费观看| 少妇被c 黄 免费观看| 天天久久久久久| 精品无码一区二区人妻久久蜜桃 | 午夜精品久久久久久久99蜜桃一| 久久久国产av美女私房| 一中国女人毛片水真多| 野狼福利社区| 亚洲成人在线乱码色午夜| 骚日日av| 果冻传媒A片一二三区| 色阁阁AV综合网| 91九九九小逼| 国产乱码精品久久久久久| 男人下部插入女人下部| 嗯嗯嗯好爽| 久久大香蕉手机高清| 久久,精品一二三| 天天日少妇逼AV| 五月天社区| 午夜寂寞欧美| 久久久久亚洲Aⅴ无码| 欧美成人一级免费电影| 日韩免费看黄片| 天堂无码精品国产久| 亚洲一区二区三区四区视频| 色情乱伦AV| 成人自拍三级在线观看| 伊人精品久久网站| 日韩性爱1级片视频| av无码av无码专区| 国产黄色 A 片免费看| 99久久久无码国产精品性男| 九九探花视频在线观看| 国产又色又爽又舒服的三级视频| 国内毛片免费h片在线| 天天色播亚洲综合网站| 久久综合九九| 亚洲精品丝袜-不卡成人免费……| 色999人与兽| 一区麻豆 高清中文字幕| 欧美色日本| 亚洲色图欧美色图制服诱惑| 日韩精品人妻中文字幕久久久| 91国内外在线| 色天天野狼综合社区| 亚洲欧美不卡线| 69国产对白刺激| 无码丰满熟妇一区二区浪潮AV| 日韩丰满熟妇| 97人人射| 91色图片| 久久久熟妇熟女国产| 精品欧美日韩在线观看| 中文字幕乱妇免费视频| 小说区 图片区色 综合区| 琪琪精品免费一区二区三区| 亚洲最大无码中文字幕网站| 五月丁香六月| 欧美乱妇狂野欧美在线视频| 久9爱经典视频| 伊人加勒比| 国产a片操逼| 亚洲精品人妻在线| 亚洲不卡AV在线| 粉嫩少妇自慰在线| 东京热,男人的天堂| 国产日本一区二区三区蜜臀在线观看| 欧美的性爱网站免费| 午夜寂寞欧美| 国内精品久久久久影院亚洲| 日本精品一区二区中文字幕| 神马福利久草| 伦理弟一页| 玖玖久久久| 激情综合网五月婷婷五月天| 99操| 91色s| 丝袜AV一区二区三区| 国产高清免费不卡av| 九九九九久久久| 亚洲欧美一区二区不卡视频播放| 91欧美经典| 精品人体无圣光凹凸| 人妻天天操天天爽视频免费| 九九九九九九九九九九九九九九九女| 性色A∨91| 可以在线观看的黄色网址| 精品妇女一区二区三区| 亚州欧美一区| 欧美成人A√在线一区二区| JULIA一区二区三区在线播放| 99九九久久| 国产精品一区二区三| 日韩视频啪啪| 欧美亚洲20p| 91天堂| 免费中文在线| 日韩精品 资源| 国产精品午夜福利| 超碰色男人操熟女| 亚洲国产剧情少妇激情| 天美传媒精品久久视频| 亚洲AV秘 精品久久老牛影视| 蜜桃精品一区二区三区ww| 东京热,男人的天堂| 国产欧美一级在线观看| 欧美日动态视频| 中文字幕三四区| 欧美拳交在线播放| 蜜桃网熟妇| 欧美一级黄色18片免费看| 九九夜精品九九在线| 国产午夜精品理论片一二三区区| 大屁股xxxxx| 国产日韩中文字幕欧美| 午夜一区二区三区国产| 午夜精品久久久久| 久9re热视频这里只有精品| 国产不卡片| 97chaopengongkai| 亚洲天堂,男人| 97操B| 99爱在线视频| 91操熟女视频 | 天天搞在线综合网| 欧美色日本| 伊人久久综合精品欧美| 操逼精品视频| 三级色综合| 蜜桃精品视频一区| 精品国产丝袜一区二区三区乱码| 好吊色一区| 26UUU欧美激情一区二区| 91伊人大香蕉| 伊人久久大香线蕉无码| 欧日韩在线观看| 1禁看欧美黄片免费看| 91无码人妻| 99日精品欧美国产| 日韩成人人妻网站| www.天天干| 97这里有精品| 国产精品视频电影| 激情五月天中文字幕色| 人妻无码后入| 青娱乐欧美激情一区二区| 日产成人久久| 午夜福利视频在线一区| 日本在线伊人啪啪| 98一区二区精品| 成人AV超碰免费在线| 97超碰亚洲| 一级A啪啪啪啪| 国产麻豆91欧美一区二区久久婷婷国产精品| 婷婷丁香五月综合| 婷婷五月天激情四射| 青青草在线视频欧美| 久久首页| 伊人网在线点播| 99热这里只有精品9| 久久久久久久强迫| 久干网| av一区二区三区 中文| 女优大全 - 91n| 天天享受天天看| 国产辣妈在线视频福利| 丝袜狂射91| 十八禁成人网站在线观看| 91亚洲丝袜熟女| 97人妻色| 亚洲 日本 国产 综合| 久久只有精品一区二区三区| 久草国产在线视频| 清纯唯美亚洲综合| 少妇同性| 色啪网| 成人a v在线播放免费| 麻豆a'v电影| 97超碰色五月| 亚洲超碰AV| 91视频综合| 俺去啦自拍| 婷婷情色综合网| 久操视频资源站公开| 88在线一区二区三区| 久久有码| 久久久久免费少妇| 国产在线视频二区| 蜜臀久久99精品久久久久久婷婷 | 天天躁日日躁AAAXX| 少妇被c 黄 免费观看| 天天搞欧美| 黄色二级片网站| 丁香色婷婷| 最新亚洲人成网站在线影院| 国产精品播放| 欧美三级不卡| 国产精品毛片| 嗯嗯啊啊操死我| 在线v中文字幕一区二区三区 | 一本一道久久综合久久| av情色影音| 91操人视频| 亚洲成熟国产精品美女| 国产午夜精品在线观看| 一区二区三区一亚洲中文字幕、综合区灬 | 大屁股人妻女教师撅着屁股| 易易A毛视频| 欧美色图91p| 江都AV在线| 国产福利在线视频网站| 国产深夜福利| 精品人妻一区春色| 中文字幕艹艹| 亚洲精品亚洲人成人网| 在线a v| 日本精品网站在线中文| 女人综合网| 国产一区在线看| 午夜精品久久久久久久| 中文字幕女同在线| 午夜福利久久久噜久噜久久综合| 亚欧精品久久久久久久久久久| 久操热线| 男人的天堂一区三区| 国语精品内射在线观看| 国产精品盗摄 偷窥盗摄| 色综合一本| 97免费在线观看| 夜夜欢天天干| 北野未奈加勒比av| 一区二区三区成人| 五月开心久久AV官网| 九九热超碰97亚洲最新香蕉| 精品v1区| 日韩情色视频| 欧美色图亚洲色图成人在在线| 综合网,亚洲,欧美| 久久人妻办公室视频| 男人的天堂不卡一区二区| 人妻中文在线| av天堂影视中文在字幕在线中文 | 91精品微拍福利| 日本性感人妻91| 亚洲天堂五月天国产| 清纯唯美综合| 97色综合中文网| 久久男人天堂| 成人羞羞视频国产| 国内亚洲高清无码| 啊啊啊啊啊啊啊在线| 日本99一区二区| 91超碰碰在线| 日本孕妇孕交| 激情人妻另类| 中文熟女五十乱码在线| 亚洲啪啪啪啪视香蕉| 欧美精品不卡一二三四在线91| 美女让帅哥通她小鸡鸡| 91热热色| 国产色产精品在线观看| 草草草视频在线免费看| 天天日天天搞天天干| 1人人看人人摸人人操| 丝袜剧情| 亚洲成人在线播放| 激情开心五月天| 泰国AV在线观看| 九九九久久久久| 熟女熟妇伦久久影院毛片一区二区| 大色网久久| 99久久99久久免费精品蜜臀| 亚洲无码国产精品久久| 波多野结衣被操50分钟免费视频| 日韩不卡毛片Av免费高清| 乱伦a片视频| 日本熟女中文| 三级片大波波| 日日夜夜草草草| 日夜尻逼网| 亚洲男人的天堂va亚洲男人社| 久草精品国产99| 日本精品不卡一二三区| 99热自拍| 极品内射| 大香蕉2017| 色色网91| 日日狠狠久久偷偷色综合免费|