在R中寻找区域年最大值

我正在处理全球最大年降水量网格数据。然而，我想通过使用一个遮罩（基于所有土地面积的145个最大值），在我的145年中，每年“仅”隔离土地面积的最大值。也就是说，当我应用遮罩时，我只接收到NA值，我无法理解为什么（当不应用遮罩时，以下过程工作正常）。以下是我迄今为止所做的工作：

Model66 <- brick("MaxPrecNOAA-GFDLGFDL-ESM2Ghistorical.nc", var="onedaymax")



#Applying the mask to isolate land areas only:

data("wrld_simpl")
b <- wrld_simpl
land <- mask(Model66,b)

#To derive highest maximum value for each layer/year for land only (145 years = 145 maximum values)

Gmax <- sapply(unstack(land), function(r){max(values(r))}) 

Gmax
[1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 
NA NA NA NA NA NA NA
[40] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 
NA NA NA NA NA NA NA
[79] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 
NA NA NA NA NA NA NA
[118] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA

Model66尝试：

Gmax <- sapply(unstack(land), function(r){max(values(r), na.rm=T)}) 

Gmax，如果您添加了。你可以让别人更容易找到并测试你问题的答案。这样你可以帮助别人来帮助你！这真的很管用！！！！谢谢你，阿曼多！这真的是一个小细节导致了这一切！非常感谢你！！！