lm
Я пытаюсь добавить набор экстраполированных "наблюдений" к матрице в R. Я знаю, как это сделать, используя обычные методы программирования (читай; кучу вложенных циклов и функций), но я чувствую, что это должно быть возможно гораздо более чистым способом, используя build in R-functionality. Приведенный ниже ...
Мой вопрос имеет отношение к оценке темпов роста численности населения в модель роста Мальтуса. В качестве примера игрушки рассмотрим набор данных игрушек df: structure(list(x= c(0L, 24L, 48L, 72L, 96L, 120L, 144L, 168L ), y = c(10000, 18744.0760659189, 35134.0387564953, 65855.509495469, 123440.067934292, ...
Возможный дубликат:нахождение места пересечения двух линейных фитингов в R Учитывая некоторые точки На графике (обычно только около 6 или 7 точек), мне нужно найти наиболее подходящее решение, где решение состоит из следующего: две линейные линии линии должны пересекаться точка пересечения (точк...
Ниже приведены 4 набора данных (я просто создал их случайным образом для обеспечения воспроизводимого кода). Я создал список из них, чтобы я мог применить "lm" к этим нескольким наборам данных сразу: H<-data.frame(replicate(10,sample(0:20,10,rep=TRUE))) C<-data.frame(replicate(5,sample(0:100,10,rep...
Я работаю с набором данных о сохранении и его влиянии на биомассу, в котором пятьдесят участков земли, каждый по одному гектару, были отобраны случайным образом с площади в десять тысяч гектаров в Северной Англии. Для каждого земельного участка были зарегистрированы следующие переменные: * Биомасса: оцен...
У меня есть вектор Y, содержащий будущие доходы, и вектор X, содержащий текущие доходы. Последнее изменение элемента на, как последнего возврата тоже в самом конце рядов. X = { 0.1, 0.3, 0.2, 0.5 } Y = { 0.3, 0.2, 0.5, NA } Other = { 5500, 222, 523, 3677 } lm(Y ~ X + Other) Я хочу убедиться, что последний...
Мой вопрос почти ответил в dplyr 0.3.0.9000 как использовать do() правильно, но не совсем. У меня есть некоторые данные, которые выглядят так: > head(myData) Sequence Index xSamples ySamples 6 0 5 0.3316187 3.244171 7 0 6 1.5131778 2.719893 8 0 7 1.9088933 3.122...
Предположим, у меня есть переменная ответа и данные, содержащие три ковариаты (как пример игрушки): y = c(1,4,6) d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2)) Я хочу подогнать линейную регрессию к данным: fit = lm(y ~ d$x1 + d$x2 + d$y2) есть ли способ, чтобы написать формулу, так что я...
Я ищу только значение коэффициента линейной модели B1(newx), а не имя. Я просто хочу значение 0,5. Я не хочу имя "newx". newx <- c(0.5,1.5.2.5) newy <- c(2,3,4) out <- lm(newy ~ newx) out выглядит так: Call: lm(formula = newy ~ newx) Coefficients: (Intercept) newx 1.5 ...
у меня есть регрессионная модель для некоторых данных временных рядов, исследующих использование наркотиков. Цель состоит в том, чтобы подогнать сплайн к временному ряду и разработать 95% CI и т. д. Модель выглядит следующим образом: id <- ts(1:length(drug$Date)) a1 <- ts(drug$Rate) a2 <- lag(a1-1) ...