R 不可见数据的多项式拟合

最初我有表格中的数据

m n
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16

我使用以下代码将其转换为

input<-file('stdin', 'r')

mn <- read.table(input, nrows = 1, as.is = TRUE)

DF <- read.table(input, skip = 0)



m <- mn[[1]]

n <- mn[[2]]

x1<- DF[[1]]

y1<-DF[[2]]

x2<-DF[[3]]

y2<-DF[[4]]

fit1<-lm(x1 ~ poly(y1, 3, raw=TRUE))

fit2<-lm(x2 ~ poly(y2, 3, raw=TRUE))

我想预测给定值之后n值的x1 y1 x2 y2的所有值

我尝试拟合lm，但我不确定如何继续处理未来要预测的所有数据点的值。如果丢失数据点，仅根据其他数据点获取系数是不够的，因为所有数据点都需要预测

x1=  1  5  9 13

x2= 2 6 10 14