作者查詢 / AndrewShi
作者 AndrewShi 在 PTT 全部看板的留言(推文), 共340則
限定看板：全部
看板排序：
2F→: l大~我知道，但是就是不知道該怎麼寫XD11/08 01:52
9F→: l大~是的，0代表晚上12點，6代表早上6點以此類推，前11/08 13:01
10F→: 面則是代表該小時颱風所在的經緯度位置。11/08 13:01
11F→: 能找一整天是最好，但我不確定要加上次日0時也有的條11/08 13:07
12F→: 件會不會使得程式很複雜，所以想說先從找一天當中四個11/08 13:08
13F→: 時段（滿18hr）皆有的去理解，之後再進階。11/08 13:08
14F→: 謝謝b大，我研究一下~11/08 13:10
16F→: 謝謝t大，我好好研究一下。11/10 23:22
17F→: T大~不好意思關於從mutate(Hour_lag =lead(Hour),...)11/20 18:34
18F→: 到最後的程式碼我不是完全很懂每一行是在做什麼（已有11/20 18:35
19F→: 查指令用途)，能否請你稍微說明一下嗎？？謝謝。11/20 18:35
21F→: 非常感謝T大，這樣我就懂了，不過還是想請問一下row_nu11/25 02:31
22F→: mber( )==1是將符合條件的時間再做重新排序（rank)的意11/25 02:31
23F→: 思嗎？！11/25 02:31
4F→: https://imgur.com/jAB83f310/31 23:05
5F→: andrew大是指改成這樣嗎？？但改成這樣會有錯誤訊息10/31 23:07
6F→: https://imgur.com/H8d5rKT10/31 23:08
7F→: 我用了c大的寫法就可以了，感謝andrew大和c大~10/31 23:31
8F→: 另外想再請教一個問題，如何我想在台灣的中心點（給定10/31 23:39
9F→: 一個經緯度位置）畫一個半徑300km的空心圓，我可以朝10/31 23:39
10F→: 哪個參數設定去想呢？？10/31 23:39
15F→: 謝謝l大和c大，我這幾天研究看看11/02 22:52
16F→: andrew大~了解，再次謝謝你11/02 22:53
17F→: https://imgur.com/BHtfQvm11/06 15:50
2F→: https://imgur.com/JUWsZdb05/24 14:51
3F→: andrew大~我如果把範圍改成>=100就可以畫的出來（下圖05/24 15:15
4F→: ），我在想可能是值過多的問題，不知道有沒有什麼方法05/24 15:15
5F→: 可以解決05/24 15:15
6F→: https://imgur.com/NMMbhc605/24 15:16
9F→: 感謝andrew大，不過想請教你用wtd.Ecdf（）如果像這樣05/24 18:08
10F→: 不分配權重給它是不是畫出來就等於原本用ecdf（）畫出05/24 18:08
11F→: 來是一樣的，只是把所有的值做分組再畫而已呢？！05/24 18:08
14F→: 什麼意思？？是要分配權重才會等於用ecdf畫出來的一樣05/24 22:44
15F→: 嗎？！05/24 22:44
19F→: 了解，權重的分配就是依據看A.B.C這三類的人數多寡去05/24 23:21
20F→: 分配對吧？！05/24 23:21
21F→: 所以andrew大所寫的程式碼所分配的權重就是依據每個區05/24 23:28
22F→: 間（10為一區間）所發生的次數去給它囉？！05/24 23:28
24F→: 了解，非常感謝andrew大，又讓我學到了不少05/24 23:47
5F→: 原來還可以先轉成矩陣，再用apply系列的指令，再次謝05/17 17:17
6F→: 謝andrew大。05/17 17:17
4F→: andrew大~對，但應該不多，analog.indices那個表格最05/09 20:41
5F→: 左邊那一列是原始資料的天數，右邊這兩列則是所有資料05/09 20:41
6F→: 裡和這天最相近的2天。05/09 20:41
7F→: 權重的表格是和analog.indices的表格相對應，也就是兩05/09 20:49
8F→: 者的［1,V1]做相乘，兩者的[1,V2]做相乘，再相加，只是05/09 20:49
9F→: 做相乘的時候analog.indices的[1,V1]不是取那個數字，05/09 20:49
10F→: 而是取那天的降雨資料出來和權重的[1,V1]這個值做相乘05/09 20:49
11F→: ，[1,V2]也是一樣。05/09 20:49
13F→: 一樣呀，第10天就是拿第10天的降雨值乘上權重[10,V1](005/09 22:15
14F→: .764485)，加上第4407天的降雨值乘上權重[10,V2］(0.2005/09 22:15
15F→: 3293)相加而得。05/09 22:15
17F→: 那是第10天被選了2次，意思是原始資料的第10天和第1105/10 01:20
18F→: 天很像，所以原始資料在第10天選了自己本身，在第11天05/10 01:20
19F→: 時也選了第10天，右邊這兩列的數字是選出跟原本資料（05/10 01:21
20F→: 最左列）那天最相近的2天出來，所以基本上會挑出自己05/10 01:21
21F→: 本身（第一列和第二列數字幾乎一樣）；然而原始資料的05/10 01:21
22F→: 第10天（最左邊的列）是沒有重複的。05/10 01:21
28F→: 表格的列號正常的話也只有9862列，後面是我新創的，忘05/10 11:11
29F→: 記刪掉，抱歉造成你理解上的誤會，再次感謝andrew大，c05/10 11:11
30F→: ode的部分我會再好好的研究。05/10 11:11
32F→: 恩恩，andrew大，想再請教你如果我每天的analog.indice05/10 13:20
33F→: s和weights都有30個的話，col.names和for迴圈的部分除05/10 13:20
34F→: 了自己補足剩下的，該怎麼修改迴圈中$i1.$i2.$w1.$w205/10 13:20
35F→: 的部分讓它自動加總30天呢？？05/10 13:20
40F→: 我想到的是第一種（雙層迴圈），不知道改成這樣對不對05/10 15:16
41F→: ，還請andrew大指點。05/10 15:16
42F→: https://imgur.com/ci54KFs05/10 15:16
45F→: 因為不知道怎麼把$i1.$i2.$w1.$w2寫成迴圈，所以把i1.i05/10 15:19
46F→: 2.w1.w2全改成i和w。05/10 15:19
48F→: 原來是要改成這樣，非常感謝andrew大，最後有2個小問題05/10 18:04
49F→: 想請教你，一個是.［1:9862]最前面的那個.是代表前面05/10 18:04
50F→: 讀進來的那個csv檔嗎？！另一個是col.names如果i有1~3005/10 18:04
51F→: 的話是只能一個一個打（命名）嗎？？05/10 18:04
54F→: 了解，再次感謝andrew大的解答和每次的幫助~05/10 23:38
12F→: 了解，不過andrew大，我改成你說的那樣反而出現error。01/25 00:54
4F→: andrew大是指obs和rnnmm被多次定義嗎?!可是我每算一次01/26 01:23
13F→: 出來的值就會被丟到矩陣裡應該沒問題吧?!未被定義的是?01/26 01:28
15F→: 跑是可以跑，只是print(rnnmm.mean)出來只會有一個值。01/26 13:59
17F→: andrew大~我最多就只想到改成上面那樣，還煩請你指點了01/26 17:11
23F→: https://imgur.com/0ur3ojG01/27 01:19
24F→: andrew大~我會出現這個error，我是用另一個觀測資料，01/27 01:21
25F→: 所以變數名稱不同，且有重跑。01/27 01:22
28F→: 了解，我有成功跑出來了，非常謝謝andrew大^^01/27 13:22
29F→: andrew大~不好意思想再請教你，我今天下午的確有跑成功01/27 21:55
30F→: 但當我改變數名稱之後再去跑就會出現下圖中最上面和最01/27 21:56
31F→: 下面的這兩個錯誤，我有上網查錯誤訊息，好像for迴圈常01/27 21:59
32F→: 會跑出下面這個錯誤訊息，還是是我改的程式碼還是有問01/27 22:00
33F→: 題呢??01/27 22:01
34F→: https://imgur.com/JiLVdxk01/27 22:02
35F→: 我有重跑和重開R跑過，但都還是會出現上面那兩個error01/27 22:06
40F→: 不瞞andrew大，老實說我有2本R的書(輕鬆學習R語言.R錦01/28 01:01
41F→: 囊妙計)，前面這本我幾乎整本看完，裡面介紹的指令我也01/28 01:03
42F→: 都有try過，但我覺得書比如介紹矩陣的部分都太粗淺，我01/28 01:06
45F→: 們真正在寫.應用的程式碼都複雜很多，不過我也承認我基01/28 01:12
46F→: 我知道，i.j不就是我資料經.緯度的維度嗎?!01/28 01:15
54F→: 了解，我好像懂了，我說一下我的理解不確定對不對。01/28 16:49
55F→: https://imgur.com/DK1E2eH01/28 16:49
56F→: 一開始先創一個跟資料總長度一樣大的矩陣(都塞NA值)，01/28 16:52
57F→: 之後再算選定(i.j)網格點的值，算完之後再把值填到原先01/28 16:54
58F→: 創的矩陣裡。01/28 16:55
62F→: 了解，謝謝andrew大01/28 21:03
2F→: 對，或許之後不一定會用到這麼長時間，但20~30年跑不11/22 22:47
3F→: 掉。11/22 22:47
6F→: 4010，這是指我的記憶體空間極限是4010Mb嗎?!11/23 01:51
9F→: 了解，所以只有換電腦記憶體的選擇一途囉?!我是有想到11/23 02:54
10F→: 可以把nc檔先轉成txt檔再轉成csv檔，這樣不知道會不會11/23 02:54
11F→: 比較不占記憶體的空間?!11/23 02:54
16F→: 了解，謝謝andrew大的解答。11/23 12:16
6F→: 我想讓365天按照日期讓R來判斷這是幾月的資料，是一月11/13 19:45
7F→: 就加總前31天的資料，是二月的就加總第32~59天的資料，11/13 19:45
8F→: 最後加總成總共12個月。11/13 19:45
9F→: 好的，我先試試看。11/13 19:45
14F→: andrew大可以給我一點提示嗎??我目前只想到，names(a)11/14 01:51
15F→: <- b，我不確定這是不是把a按b分組的意思。11/14 01:51
17F→: andrew大~你給我的這段程式碼是把365天按照2018年的月11/14 12:34
18F→: 份和日數去分組，最後再從1+2+3+...365除上每個月的天11/14 12:34
19F→: 數的意思嗎?!11/14 12:34
23F→: andrew大~我稍微改了一下上面的程式碼，但還是不行，11/14 15:52
24F→: 我也覺得還是怪怪的，但不知道該怎麼改，tmonth有36511/14 15:52
25F→: 個數字(1~12)，我不知道該怎麼把都是1的那31天的降雨11/14 15:53
26F→: 值放入同一格裡。11/14 15:53
27F→: https://i.imgur.com/MQFzcpY.jpg11/14 15:54
33F→: andrew大~打架是指??我的最後2個迴圈是在做降雨值的加11/14 23:42
34F→: 總，能否請教andrew大我的迴圈有缺或錯哪些地方嗎??如11/14 23:43
35F→: 果沒有錯的很離譜，我想看自己能不能修改並完成它(畢11/14 23:43
36F→: 竟是用自己的邏輯寫的)，至於你給我的那段程式碼，我也11/14 23:43
37F→: 只是勉強看的懂，但我還真的不知道也想不到要那樣寫。11/14 23:43
48F→: 好的，那我明天再來好好的想，謝謝andrew大。11/15 01:45
52F→: 了解，非常謝謝andrew大，我的想法的確是你後面傳給我11/22 12:31
53F→: 的這個寫法，再看了你的程式碼之後我有想出來要怎麼寫11/22 12:31
54F→: 了。11/22 12:31
55F→: andrew大~想請問你確認一下這個迴圈概念套用在兩年或是11/22 17:15
56F→: 多年的資料上應該也可以吧?!就是算兩年(或多年)一月份11/22 17:15
57F→: 的平均降雨值11/22 17:15
4F→: andrew大~好像畫不出來，畫出來是下圖這樣子，且有等了11/08 18:35
5F→: 超過10分鐘，圖還是這樣(空白)。11/08 18:35
6F→: https://imgur.com/pzSSKip11/08 18:36
8F→: 阿...我知道了，lon和lat不能刪掉，xlim和ylim加在後11/08 18:41
9F→: 面就好，謝謝andrew大^^11/08 18:41
1F推: 非常感謝andrew大，願意花這麼多時間幫我解答，我今天10/23 11:32
2F→: 也會好好研究程式碼，如果對於你寫的程式碼有疑惑的話10/23 11:32
3F→: 方便能再請教你嗎XD??10/23 11:32
4F推: abdrew大~我想請教apply裡的c(1,2)是指把每列.每行(每10/23 15:24
5F→: 個經度.緯度)的降雨值都帶到function裡面的意思對嗎?!10/23 15:24
6F→: 另外想請問你是怎麼把時間(月)放到precip.ave.monthly10/23 15:24
7F→: 的第一個維度裡的呢??10/23 15:24
8F推: 感謝andrew大詳細的解答，上面的敘述我能理解，不過我10/24 14:47
9F→: 很好奇為什麼在單用apply算降雨年平均的時候時間是在10/24 14:47
10F→: 降雨的第三個維度，而在算月平均的時候，時間變成降雨(10/24 14:47
11F→: precip.ave.monthly)的第一個維度，而經.緯度則變成第10/24 14:47
12F→: 二.三個維度，我看不出來哪一段程式碼是在做這樣的處理10/24 14:47
13F→: 。10/24 14:47
18F推: 了解，真的非常感謝andrew大，讓我學到很多~10/24 22:25
19F推: andrew大~想再請教你如果想用迴圈的概念來寫的話，我想10/25 16:22
20F→: 到的寫法是:for(i in (1:1440)){10/25 16:22
21F→: for(j in (1:400)){10/25 16:22
22F→: mean(precip[i,j,time=(1:31)])}}10/25 16:22
23F→: 但跑出來也只有一個值，所以想請教你我的迴圈寫法(概10/25 16:22
24F→: 念)是哪裡有出錯嗎??10/25 16:22
29F推: 了解，非常感謝~10/25 17:09
30F推: andrew大~不好意思，想再請教你一個基本問題，把結果10/26 17:28
31F→: 填到一個矩陣裡，matrix(mean(precip[i,j,time=(1:31),10/26 17:28
32F→: 1440,400))這樣寫對嗎??10/26 17:28
33F推: 我懂了，首先先創一個數字由1~24，2列3行共4個(矩陣)10/27 17:21
34F→: 的陣列a，之後再創一個矩陣m，列和行的數目和陣列a的10/27 17:21
35F→: 第一和第二個維度一樣(2列3行)，最後再把陣列a每個列和10/27 17:21
36F→: 行各別的值(共4個)相加取平均後放到矩陣m裡(如:(1+7+1310/27 17:21
37F→: +19)/4=10)，我這樣的理解應該沒有錯吧 :)?!10/27 17:21
38F推: 而我的程式應該改為:precip1 <- matrix(precip,1440,4010/27 17:30
39F→: 0)10/27 17:30
40F→: for(i in (1:1440)){10/27 17:30
41F→: for(j in 1:400)){10/27 17:30
42F→: precip1[i,j] <- mean(precip[i,j,time=(1:31)])}}，非10/27 17:30
43F→: 常謝謝andrew大用引導式的方式教我，其實我也比較喜歡10/27 17:30
44F→: 用這種方式來學習。10/27 17:30
49F推: 抱歉，andrew大，剛剛才看到你最後的回覆，所以我要改11/02 16:31
50F→: 成precip1 <- matrix(precip1,1440,400)這樣才對嗎?!11/02 16:31
56F推: 了解XD，所以我的precip1要改成NA_real_，先把這個建11/02 18:16
57F→: 立出來空的矩陣(precip1)先(隨便)填個值(NA或0)，之後11/02 18:16
58F→: 再把算出來的平均降雨值填到這個矩陣裡。11/02 18:16
60F推: andrew大~我想請教你一個問題，如果畫圖的部分我想畫特11/07 17:19
61F→: 定區域的話(如:東亞)，改image裡的經.緯度範圍和月平11/07 17:19
62F→: 均降雨裡的經.緯度範圍可以畫的出來，但是之後再疊加世11/07 17:19
63F→: 界地圖(有設成同樣範圍)的時候畫不出來(它有顯示有多11/07 17:20
64F→: 少個警告訊息)，但我打warnings( )，它也沒列出警告訊11/07 17:20
65F→: 息，想請教你這個問題該怎麼解決呢??11/07 17:20
67F推: 好的。11/08 12:01