小编Hav*_*v0k的帖子

科学名称通常由3个信息组成:属,种上皮和作者.一个简单的例子如下:

Acanthus ilicifolius L.

• 属:Acanthus
• 物种上皮:ilicifolious
• 作者:L.

简单.然而,当我们必须处理杂交种,亚种/品种/形态,几个作者和其他不一致时,问题变得更加复杂.在这些情况下,物种名称可能如下所示:

比照 穿心莲(Burm.f.)墙.前Nees

• cf.:物种未确定100%确定
• 属:穿心莲
• 物种上皮:paniculata
• 作者:( Burm.f.)Wall.前Nees

或这个:

Ipomoea pes-caprae(L.)DC.亚种.brasiliensis(L.)Ooststr.f

• 属:Ipomea
• 物种上皮:pes-caprae
• 物种作者:(L.)DC.
• 亚种上皮:brasiliensis
• 亚种作者:(L.)Ooststr.f

我正试图找到一种可靠的解构这些名字的方法.如果if/else语句我可以使用吨写一些hackish代码,但我正在寻找更优雅(和健壮)的东西.我想的是某种解析器的名称类似于解析数学表达式的计算器.不幸的是,我不是最复杂的程序员,我之前也没有写过真正的解析器,也不知道在这种情况下它是否有意义,因为科学名称有很多变化.你认为解决这个问题的最佳方法是什么？首选语言是R,如果它更适合任务,也许也是Julia.

我的模型包括一个响应变量,五个预测变量和一个预测变量_1和预测变量_2的交互项.我想为每个预测变量绘制部分残差图,我通常会使用crPlots包中的函数来实现car.不幸的是,该函数抱怨它不适用于包含交互项的模型.

有没有其他方式做我想要的？

编辑:我创建了一个说明问题的小例子

require(car)
R <-  c(0.53,0.60,0.64,0.52,0.75,0.66,0.71,0.49,0.52,0.59)
P1 <- c(3.1,1.8,1.8,1.8,1.8,3.2,3.2,2.8,3.1,3.3)
P2 <- c(2.1,0.8,0.3,0.5,0.4,1.3,0.5,1.2,1.6,2.1)

lm.fit1 <- lm(R ~ P1 + P2)
summary(lm.fit1)
crPlots(lm.fit1) # works fine

lm.fit2 <- lm(R ~ P1*P2)
summary(lm.fit2)
crPlots(lm.fit2) # not available

我想根据 HTML 标签查询我的 Shiny 应用程序中的元素。例如，对于以下应用程序，如何访问所有h1元素，包括其标签和 id？

library(shiny)
ui <- fluidPage(
  h1("Get", id = "1"),
  h1("this", id = "2"),
  h2("But", id = "3"),
  h3("not", id = "4"),
  h4("that", id = "5")
)

server <- function(input, output) {}

shinyApp(ui = ui, server = server)