modelingtoolkitize fails on ODEProblem where LHS is R code about modelingtoolkit.jl HOT 4 OPEN

using OrdinaryDiffEq
using ModelingToolkit
using Symbolics
using RCall

RCall.sexp(::Type{RCall.Sxp}, x::Num) = RCall.sexp(x.val)
RCall.sexpclass(::Num) = RCall.Sxp
RCall.sexp(::Type{RCall.Sxp}, x::Term{Real, Base.ImmutableDict{DataType, Any}}) = RCall.sexp(eval(x))
RCall.sexpclass(::Term{Real, Base.ImmutableDict{DataType, Any}}) = RCall.Sxp

R"""
sir_ode_r <- function(u,p,t){
    S <- u[1]
    I <- u[2]
    R <- u[3]
    N <- S+I+R
    beta <- p[1]
    cee <- p[2]
    gamma <- p[3]
    dS <- -beta*cee*I/N*S
    dI <- beta*cee*I/N*S - gamma*I
    dR <- gamma*I
    return(c(dS,dI,dR))
}
"""

function sir_ode_jl(u,p,t)
    robj = rcall(:sir_ode_r, u, p, t)
    return convert(Array,robj)
end

δt = 0.1
tmax = 40.0
tspan = (0.0,tmax)
u0 = [990.0,10.0,0.0] # S,I,R
p = [0.05,10.0,0.25] # β,c,γ
prob = ODEProblem{false}(sir_ode_jl, u0, tspan, p)
@named sys = modelingtoolkitize(prob)