Question

我尝试为硅获得Cauchy's dispersion equation的前三个系数。使用包含某些波长折射率的csv（您可以找到here），我尝试拟合以下模型：

library(readr)
library(tidyverse)
library(magrittr)
library(modelr)
library(broom)
library(splines)

# CSV parsing
RefractiveIndexINFO <- read_csv("./silicon-index.csv")

# Cleaning the output of the csv-parsing
indlong = tibble(RefractiveIndexINFO$`Wavelength. µm`,RefractiveIndexINFO$n)
names(indlong) = c('w','n')

# Remove some wavelengths that might not fit
indlong_non_uv = indlong %>% filter(indlong$w >= 0.4)

# Renaming variables
w = indlong_non_uv$w
n = indlong_non_uv$n

# Creating the non linear model
model = nls(n ~ a + b*ns(w,-2) + c*ns(w,-4), data = indlong_non_uv)

# Gathering informations on the fitted model
cor(indlong_non_uv$n,predict(model))
tidy(model)

出现以下错误：

Error in c * ns(w, -4) : non-numeric argument to binary operator

如何绕过这种情况并连续得到三个系数（a，b，c）？

显然，使用model = nls(n ~ a + b*ns(w,-2), data = indlong_non_uv)不会出错。

Answer 1

试试这个：

library(readr)
library(tidyverse)
library(magrittr)
library(modelr)
library(broom)
library(splines)

# CSV parsing
RefractiveIndexINFO <- read_csv("aspnes.csv")
RefractiveIndexINFO <- RefractiveIndexINFO[1:46,]
RefractiveIndexINFO <- as.data.frame(apply(RefractiveIndexINFO,2,as.numeric))
names(RefractiveIndexINFO) <- c('w','n')

indlong_non_uv = RefractiveIndexINFO %>% filter(RefractiveIndexINFO$w >= 0.4)

# Creating the nonlinear model
model <- nls(n ~ a + b*w^(-2) + c*w^(-4), data = indlong_non_uv,
   start=list(a=1, b=1, c=1))

# Gathering informations on the fitted model
cor(indlong_non_uv$n,predict(model))
# [1] 0.9991006

tidy(model)
#   term    estimate   std.error statistic      p.value
# 1    a  3.65925186 0.039368851 92.947896 9.686805e-20
# 2    b -0.04981151 0.024099580 -2.066904 5.926046e-02
# 3    c  0.05282668 0.003306895 15.974707 6.334197e-10

或者，您可以使用线性回归：

model2 <- lm(n ~ I(w^(-2)) + I(w^(-4)), data = indlong_non_uv)
summary(model2)

# Coefficients:
#              Estimate Std. Error t value Pr(>|t|)    
# (Intercept)  3.659252   0.039369  92.948  < 2e-16 ***
# I(w^(-2))   -0.049812   0.024100  -2.067   0.0593 .  
# I(w^(-4))    0.052827   0.003307  15.975 6.33e-10 ***

与R非线性拟合

1 个答案: