หมายเหตุ: หากไม่มีคำตอบที่ถูกต้องหลังจากผ่านไปหนึ่งเดือนฉันได้โพสต์ข้อความใหม่เป็นSO
พื้นหลัง
ฉันมีแบบจำลองที่Y = f ( X )
คือเมทริกซ์ n × mของตัวอย่างจากพารามิเตอร์ mและ Yคือ n × 1เวกเตอร์ของเอาต์พุตโมเดล
มีการคำนวณอย่างเข้มข้นดังนั้นฉันจึงต้องการประมาณ fโดยใช้ลูกบาศก์หลายตัวแปรผ่านจุด ( X , Y )เพื่อให้ฉันสามารถประเมิน Yที่จุดจำนวนมาก
คำถาม
มีฟังก์ชั่น R ที่จะคำนวณความสัมพันธ์โดยพลการระหว่าง X และ Y หรือไม่?
โดยเฉพาะฉันกำลังมองหารุ่นหลายตัวแปรของsplinefun
ฟังก์ชั่นซึ่งสร้างฟังก์ชั่น spline สำหรับกรณี univariate
เช่นนี้เป็นวิธีการsplinefun
ทำงานสำหรับกรณี univariate
x <- 1:10
y <- runif(10)
foo <- splinefun(x,y)
foo(1:10) #returns y, as example
all(y == foo(1:10))
## TRUE
สิ่งที่ฉันได้ลอง
ฉันได้ตรวจสอบแพ็กเกจmdaและดูเหมือนว่าสิ่งต่อไปนี้จะใช้งานได้:
library(mda)
x <- data.frame(a = 1:10, b = 1:10/2, c = 1:10*2)
y <- runif(10)
foo <- mars(x,y)
predict(foo, x) #all the same value
all(y == predict(foo,x))
## FALSE
แต่ฉันไม่สามารถหาวิธีที่จะใช้ลูกบาศก์อิสระ mars
อัปเดตตั้งแต่เสนอเงินรางวัลฉันเปลี่ยนชื่อ - ถ้าไม่มีฟังก์ชั่น R ฉันจะยอมรับตามลำดับที่ต้องการ: ฟังก์ชั่น R ที่ส่งออกฟังก์ชั่นกระบวนการเกาส์เซียนหรือฟังก์ชั่นการแก้ไขหลายตัวแปรอื่นที่ผ่านจุดการออกแบบ ใน R, Matlab อื่น