Подмножество кадра данных на основе нескольких условий в R - PullRequest
Новая
       9

Подмножество кадра данных на основе нескольких условий в R

0 голосов
/ 05 января 2019

Я пытаюсь разместить большой массив данных (следы циклонов в северном полушарии) при нескольких условиях: данные ниже 

centro <- read.table("https://forms.naturwissenschaften.ch/imilast/_ERAinterim_1.5_1979_MTEX/ERAinterim_1.5_NH_M02_19790101_20121231_MTEX.txt?_ga=2.18919096.1825595846.1546710263-1112023567.1546710263", sep="", fill = T, nrows = 500,
                 header = F, skip = 2) # read here only the first 500 rows

centro <- na.omit(centro)

colnames(centro) <- c("Code","CycloneNo","StepNo","DateI10","Year","Month","Day","Time","LongE","LatN","Intensity1","Intensity2","Intensity3")

Я хочу установить подмножество только циклонов (основанных на уникальном столбце CycloneNo), которые сформированы в пространственном поле (например, -4 и 40 E долготы и 32-45 N широты), когда столбец StepNo == 1. Обычно, это было бы легко сделать: 

centro_subs <- centro[centro$LongE>=-4 & centro$LongE <= 40 & centro$LatN>= 32 & centro$LatN <= 45,]

Однако я хочу сохранить только циклоны, сформированные в этом блоке (когда StepNo == 1), но остальные их треки тоже вне этого блока.

Я попытался получить это, выполнив: 

df_s <- centro[1,]
df_s[1,] <- NA # create an empty dataframe to be filled in the iteration


for (i in 1:length(unique(centro$CycloneNo))){
print(i)
a <- centro[centro$LongE[centro$StepNo==1]>= -4 & 
centro$LongE[centro$StepNo==1] <= 40 & 
centro$LatN[centro$StepNo==1]>= 32 & centro$LatN 
<=45[centro$StepNo==1],]
df_s <- rbind(a, df_s)
}

Однако это заканчивается пустым кадром данных, заполненным NA. Я знаю, что это довольно сложно описать здесь. Я чувствую, что я как-то близок, но я также очень устал сейчас, чтобы попытаться найти новые пути.

Ответы [ 2 ]

0 голосов
/ 07 января 2019

Иосиф дает хороший ответ. В качестве альтернативы можно использовать это в data.table, что может обеспечить немного большую читабельность за счет некоторой скорости. 

centro <- data.table(centro)
centro[CycloneNo %in% CycloneNo[StepNo == 1 & 
                                  LongE %between% c(-4,40) & 
                                  LatN %between% c(32,45)]]
0 голосов
/ 06 января 2019

Я не думаю, что вы хотите петлю. Я уверен, что это не самый элегантный подход, но я думаю, что он работает. 

step1s <- subset(centro, StepNo == 1) # only take step 1 of all cyclones
keeps <- step1s$CycloneNo[step1s$LongE>=-4 & step1s$LongE <= 40 & step1s$LatN>= 32 & step1s$LatN <= 45] # find cyclone numbers for cyclones meeting the condition
centro_sub <- centro[centro$CycloneNo %in% keeps, ] # keep all steps of cyclones meeting the conditions
...