《新加坡聯合早報》12/16 （西西里島唱播）
* 《紐約時報》週二（15日）引述兩名知情人士的話，為確保輝瑞可以在明年上半年為美國生產更多冠病疫苗，川普政府希望動用其權力，要求藥廠釋出更多疫苗原材料。
如果雙方達成協議，這將能局部緩解美國冠病疫苗供應不足的處境。
美國總統川普與輝瑞磋商爭取明年取得更多冠病疫苗 https://www.zaobao.com/realtime/world/story20201216-1109212?utm_source=ZB_iPhone&utm_medium=share
* 美國候任總統拜登的過渡團隊週二（15日）宣佈，拜登提名民主黨總統初選前候選人布蒂吉格出任交通部長。
拜登將設氣候變化辦公室 布蒂吉格出掌交通部 https://www.zaobao.com/realtime/world/story20201216-1109191?utm_source=ZB_iPhone&utm_medium=share
* 英國國會週二提出立法提案，警告Facebook、推特和中國的TikTok如果不能刪除和限制兒童性虐待和恐怖主義等非法內容的傳播，將面臨高達營收10%
英國或對TikTok開出巨額罰款 https://www.zaobao.com/realtime/china/story20201216-1109195?utm_source=ZB_iPhone&utm_medium=share
* 俄羅斯政府週二（15日）發表聲明稱，俄政府已經批准中國石化參與一個大型石化項目。這與莫斯科加強與北京聯繫的計劃一致。
俄羅斯批准中石化參與遠東地區大型聚合物項目 https://www.zaobao.com/realtime/china/story20201216-1109192?utm_source=ZB_iPhone&utm_medium=share
* 英國首相強森將在1月訪問印度，這是他自去年上任以來對外國進行的首個重要訪問。
英首相將在1月訪問印度 https://www.zaobao.com/realtime/world/story20201216-1109197?utm_source=ZB_iPhone&utm_medium=share
* 新加坡和美國簽署諒解備忘錄，深化經濟合作，擴大給予兩國企業的貿易融資和投資援助。
新加坡與美國簽署備忘錄 支持企業經貿投資 https://www.zaobao.com/zfinance/realtime/story20201216-1109199?utm_source=ZB_iPhone&utm_medium=share
* 在中芯國際被納斯達克剔除股指後，該公司高層人士出現變動，聯合CEO梁孟松辭職，74歲蔣尚義回歸任副董事長。
中芯國際高層變動 梁孟松辭職 蔣尚義回歸 https://www.zaobao.com/realtime/china/story20201216-1109200?utm_source=ZB_iPhone&utm_medium=share
* MSCI明晟週二表示，將把10家中資公司的股票從旗下部分指數中剔除，此前美國發佈了限制購買這些公司股票的行政令。不過MSCI明晟稱，該公司將推出幾個保留這些中資公司股票的新版本指數。
明晟將推出保留被剔除中企的新版指數 https://www.zaobao.com/realtime/china/story20201216-1109210?utm_source=ZB_iPhone&utm_medium=share
* 英國首都倫敦週三（16日）將當地防疫限制措施上調至最高的第三級，這對希望趕在聖誕佳節期間吸引更多顧客上門的酒店、餐飲等服務業者來說是一大打擊。
倫敦將防疫限制措施上調至最高 https://www.zaobao.com/realtime/world/story20201216-1109205?utm_source=ZB_iPhone&utm_medium=share
* 前民主黨立法會議員許智峰以公務外訪名義往丹麥棄保潛逃往英國倫敦，之後向宣布表明未必有機會再返香港，但至今尚未辭去中西區區議員職務，被質疑浪費公帑，港府停發薪資。
港府停發許智峰區議員薪酬 https://www.zaobao.com/realtime/china/story20201216-1109209?utm_source=ZB_iPhone&utm_medium=share
* 中國央行金融穩定局局長孫天琦指出，第三方互聯網金融平台銷售銀行存款產品的業務，屬「無照駕駛」的非法金融活動，應納入金融監管範圍；對此類業務，未來需要明確銀行准入資質和標準、研究出台針對高風險銀行吸收存款行為的有關法律法規、完善審慎監管指標和有關規則等。
中國央行：互聯網平台存款業務屬「無照駕駛」 https://www.zaobao.com/realtime/china/story20201216-1109211?utm_source=ZB_iPhone&utm_medium=share
* 美國選舉人團昨天（15日）投票，民主黨籍拜登獲302張選票，跨越當選門檻，確認成為下一任美國總統。台灣外交部發言人歐江安表示，駐美國代表處已轉致蔡英文賀函，指台灣會以適當方式表達對拜登的誠摯祝賀。
拜登勝選 台駐美處已轉發賀函 https://www.zaobao.com/realtime/china/story20201216-1109196?utm_source=ZB_iPhone&utm_medium=share
* 中國月球探測器「嫦娥五號」月面採樣返回任務即將結束，將成為繼美國和前蘇聯之後，第三個取得月球樣本的國家。嫦娥五號預料明天（17日）凌晨回到地球，在內蒙古中部草原的四子王旗著陸場著陸。
嫦娥五號明天回航 將於內蒙古著陸 https://www.zaobao.com/realtime/china/story20201216-1109207?utm_source=ZB_iPhone&utm_medium=share
* 中國大陸解放軍軍機近期連續進入台灣鄰近空域，今天（16日）早上8時許再次進入台灣西南空域。本月累計14天進入台灣空域，台灣空軍廣播驅離紀錄共計33次。
台空軍本月33次驅離大陸軍機 https://www.zaobao.com/realtime/china/story20201216-1109215?utm_source=ZB_iPhone&utm_medium=share
* 新加坡16歲及以上國人和長期居民，包括客工等長期工作准證持有人，預計明年底可免費完成冠病疫苗接種，受訪的年長者對此表示歡迎，並表示願意前往注射。
但不少中年人則對注射疫苗存有疑慮，因此將先觀察一陣子，確保疫苗不會產生副作用才注射。
受訪年長者願接種冠病疫苗 中年人要觀察再決定 https://www.zaobao.com/znews/singapore/story20201216-1109043?utm_source=ZB_iPhone&utm_medium=share
* 【港國安法下銀行賬戶大轉移】《香港國安法》實施快要半年，近日在香港引起了「銀行戶頭大轉移」這個熱門話題。「黃絲」或其支持者要尋求更安全的儲錢方法；港府高官則讓家人開立戶頭代理薪水。
不少所謂「專家」立即在網上支招，除了效法特區行政長官林鄭月娥把現金放在家中，還有七種方法可保留個人的身家，包括在港開立虛擬銀行戶頭、海外投資以及買加密貨幣等。
* 【林鄭被中央要求押後上京述職 分析指連任機會大減】香港特首按慣例每年12月都會到北京向中央領導人述職，但原定上周訪京述職的林鄭月娥，卻首次被中央要求押後上京。有受訪學者認為，個中原因很可能是北京不滿林鄭月娥近來的抗疫表現，意味林鄭連任的機會將大減。
* 【福奇建議拜登和賀錦麗盡快接種冠病疫苗】美國正在進行大規模的冠病疫苗接種活動，美國首席傳染病專家福奇「強烈建議」候任總統拜登及候任副總統賀錦麗盡快接種。
《紐約時報》報道，即將成為拜登的首席醫療顧問的美國國家過敏症和傳染病研究所所長福奇週二（15日）說：「出於安全考慮，我強烈認為，他們得盡快接種，越快越好。我們必須確保他（拜登）在1月上任時已得到充分保護。」
* 【入冬供暖對煤炭需求增 學者分析中國不會全面禁澳煤進口】中澳緊張關係又現新火頭，中國可能禁止澳大利亞煤炭進口的消息再次傳出後，澳大利亞指責中國違反世界貿易組織（WTO）規則與兩國自貿協定，並說會尋求所有途徑捍衛該國企業的權益。
不過，受訪中國學者認為，隨著入冬開始供暖，對煤炭的需求相應增加，中國全面禁止澳洲煤炭進口的可能性不大。
* 【越南中全會即將選出新領導層 權力鬥爭漸升溫】越南執政的共產黨為期一周的中央委員會全體會議昨天開始，此次會議看點是新一屆領導層將出爐，並為未來五年的國家發展定下基調。
越共中央總書記、國家主席阮富仲是越南數十年來權勢最大的領袖之一，他在新領導層中或保留一席重要職位。然而，可能存在變數的是他的健康狀況。76歲的阮富仲近幾個月公開露面時氣色不佳，其健康狀況引起揣測。
* 【新加坡學者：要有效控制病毒傳播 至少八成國人須接種】新加坡政府不會強制所有國人接種疫苗，但是專家預計，須要有至少八成的國人接種疫苗，才能有效控制病毒的傳播。
杜克—新加坡國立大學醫學院新發傳染病重點研究項目副主任黃英勇教授在接受《聯合早報》訪問時說，由於一些人可能對疫苗出現過敏反應，或正在使用遏制免疫系統的藥物，因此不是所有人都適合接種疫苗。對此，病毒的感染傳染數（也稱R數值）可幫助專家判斷需要多少人接種疫苗，才能停止病毒傳播。
據他估計，冠病的R數值為3，意味著每名確診者，會將病毒傳染給另外三人。他說：「根據這樣的數據，我們必須讓至少66％的國人接種疫苗。不過我認為最好是能讓至少80％的國人接種，才能控制病毒的傳播。
* 【新加坡有完整物流設施 運輸和儲存大量疫苗】新加坡具備完整和成熟的物流設施來運輸和儲存大量冠病疫苗，但鑒於安排醫療機構注射疫苗的流程相當複雜，當局仍在研究如何實施疫苗注射計劃，日後公佈更多詳情。
教育部長兼財政部第二部長黃循財昨天在抗疫跨部門工作小組記者會上說：「我明白大家都很關注何時能讓首批國人接種疫苗，但從物流、冷鏈到運送至接種場所，過程相當複雜。從冷鏈取出疫苗後，必須在短時間內接種。」
黃循財透露，政府正同相關單位制定接種計劃，一旦敲定細節，就會分階段在明年落實。
* 【孟加拉成千上萬回教徒 不顧疫情出席葬禮】孟加拉宗教組織「保衛伊斯蘭聯盟」（Hefazat-e Islam）的75歲副首領卡塞米日前感染冠病不治，其葬禮週一在首都達卡舉行，吸引了成千上萬人出席，把街道擠得水洩不通。孟加拉和世界各地的政府一樣，告誡人們不要舉行大型集會，但當天出席葬禮的人不但沒有遵守社交距離，他們大部分也沒有戴口罩，恐怕葬禮後又再掀起另一波疫情。
* 【日本上週末日增3000病例 菅義偉暫停實施振興旅遊方案】日本全國單日確診病例在週末破3000起，達新高。東京、大阪、名古屋和北海道也因疫情反彈，被列入必須繼續加強防疫的重點區。
首相菅義偉昨天決定從12月28日到明年1月11日，全面暫停實施振興旅遊方案「Go To Travel」（去旅行）。
他在記者會上說，全國的確診病例持續居高不下，疫情擴大的地區不斷增加，已經造成醫療機構和前線醫療人員的負擔增加。
據報道，他原本只考慮追加東京和名古屋，但又唯恐確診數據擴大，醫療單位承受不了。因此，決定擴大範圍到全國，要求民眾別出遠門，將新年假期作為「集中抗疫期」。
r 刪除變數 在 Re: [問題] 用R做複線性迴歸方程分析- 看板R_Language 的推薦與評價
※ 引述《sseeaann (屌哥)》之銘言：
: 資料檔我是載入套件faraway 用data(savings)資料檔去分析
: 指令大部分都了解 但可能我迴歸觀念不太好 有些步驟會搞混
: ex:像是離群值要常態假設前刪除還是確定為常態假設再刪除
: 或是刪除變數該用stepwise或是用full model和reduced model的F分配來檢定
: 這些步驟有點不知道用的時機點在哪裡
: 有沒有比較懂迴歸分析的大大可以稍微詳述一下整套迴歸的步驟(利用此資料檔)?
: 另外小弟我晚一點會在把我打的程式碼PO上來 目前正在打當中^^
: → sseeaann:我在統計版也有發文 不過都沒人回-.-
: → sseeaann:我現在比較想知道的是說 當整體模型檢定出現兩個以上變數
: → sseeaann:T值>0.05 不拒絕Ho 那是否可以直接剔除掉?
: → sseeaann:如果直接剔除 是要用哪個方法 偏F檢定還是用stepwise Reg
你想要了解的部分就是regression analysis中
model checking 跟 model selection的部分
我就分這兩個part來回答 一邊附上我的程式做解說
library(faraway)
data(savings)
# Part I. Model checking
# 一般來說，model checking，主要就是error term的常態假設以及等變異假設的確認
# 另外，還有利用leverage or Cook's distance查看有沒有outliers
# 在R裡面有很方面的方式去做這三件事，程式如下：
lm.fit = lm(sr~pop75+pop15+dpi+ddpi, data=savings)
summary(lm.fit)
par(mfrow = c(2, 2), oma = c(0, 0, 2, 0))
plot(lm.fit)
# a. Normality Checking
# 右上角就是確認你的殘差，通常就是直接看圖片大概跟QQ-line (圖中的broken line)
# 接近，大概我們都會認為error term符合常態假設，至於檢定，就我的經驗來說，
# 我實在不認為Goodness of fit檢定有良好的Power.....
# 常見的檢定像是：Shapiro-Wilk test of normality
# R當然有相對的程式如下：
shapiro.test(resid(lm.fit))
# H_0: Normality
# 或是你要自己畫QQ plot，我這提供一個範本
library(lattice)
qqmath(~ res, data = data.frame(res=resid(lm.fit)), distribution = qnorm,
aspect="xy",prepanel = prepanel.qqmathline, panel = function(x, ...) {
panel.qqmathline(x, ...)
panel.qqmath(x, ...)
})
# b. Equal Variance
# 左上角的圖就是用來check這個假設，一般來說，不要有特別的pattern
# 像是漸小，漸大，先小後大等等，或是有二次曲線出現等
# 通常就認為假設是OK的，我一樣不喜歡等變異假設
# 一樣提供R的等變異假設檢定，在lmtest這個package裡面有
library(lmtest)
gqtest(sr~pop75+pop15+dpi+ddpi, data=savings) # 未測試，直接google的
# equal variance圖的其他畫法
xyplot(res~fit, data=data.frame(res=resid(lm.fit), fit=fitted(lm.fit)))
# c. Outliers checking
# 通常就看右下角跟左下角兩張圖
# 左下角那張 通常不要超過sqrt(3)=1.732就可以了
# 右下角那張看leverage跟Cook's distance，但是至今我還是不太會看...
# 我記得我之前學到的是通常不要超過average leverage的0.7倍吧
# 不過R會自動幫你註記他認為有問題的點，機制我就不知道了
# 其他畫法
xyplot(stdres~fit,data=data.frame(fit=fitted(lm.fit),stdres=rstandard(lm.fit)))
# d. others
# 獨立性檢定就很少再用，因為通常都只能檢驗特定期數，很少用
# R code如下：
library(car)
durbinWatsonTest(lm.fit)
# Part II. Model selection
# 這個part就很複雜了，因為還要牽扯到每一個應用層面的關係
# 通常不同的應用層面，有不同的選擇標準以及方式
# 常見的有stepwise, backward, forward
# 標準有R^2, AIC, BIC, p-value等
# 通常先進行這個Part，才進行model checking的部分
# 我這裡是因為它比較複雜先放在後面說
# 常見的code像是：
step(lm.fit)
step(lm.fit, direction="backward")
step(lm.fit, direction="forward")
library(MASS)
stepAIC(lm.fit)
# 共線性的檢驗
vif(lm.fit)
# 如果大於5，就是有共線性的存在
# 便要去確定covariate的correlation matrix
# 通常不建議有超過0.7相關係數的變數存在
# 其他解決方案則有 lasso 跟 ridge等方式
# 至於留存問題，通常根據你對資料的了解而決定
# 但是這個部分，真的不是我在這裡可以說完的，我就不獻醜了....
# 有錯煩請指教，獻醜了，謝謝
剛有google到下面兩個網址，你可以去瞧瞧
https://www.statmethods.net/stats/regression.html
https://www.statmethods.net/stats/rdiagnostics.html
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.164.79.227
※ 編輯: celestialgod 來自: 218.164.79.227 (09/05 22:14)
# confidence interval for fitted value
predict(lm.fit, interval="confidence", level=0.95)
# confidence interval for prediction
pre.data = matrix(c(30,10,700,4,40,5,1000,2),2,,TRUE)
rownames(pre.data) = c("A_region", "B_region")
newdata=data.frame(pre.data)
names(newdata)=names(savings)[2:5]
predict(lm.fit, newdata=newdata ,interval="prediction", level=0.95)
# To ger more information by R>?predict.lm
confint(lm.fit)
# 但是就我所學來說，刪除離群值這件事真的很少做...
# 我是念純統計的，我們老師通常都說資料最大
# 如果不是記錄錯誤的話，通常都會納入考量...
# 或是influence point才會考慮刪除。
我看了一下，那個只能用在category的covariate...
連續型的covariate不能用那個檢驗等變異....
※ 編輯: celestialgod 來自: 218.164.79.227 (09/06 01:15)
... <看更多>