[R] Creating a loop with code from the mblm package

[Bailey Hewitt]

Hello,


I am a novice in coding in R and have come across an error I am having a hard time fixing. I am trying to use the mblm package to run a Theil-Sen linear model. The code for this function is:

mblm(Y ~ X, dataframe, repeated = FALSE)

My goal is to put this into a loop so that I can calculate the Theil-Sen slope of each column in my csv. file. The file contains one column of years (x value) and 3 columns of days of the year (y values). All columns are the same length. The code I currently have is:


read.csv("~/Documents/NH- Lake Mendota_SenSlope_Data2.csv", header = TRUE, sep = ",")

mydata= read.csv("~/Documents/NH- Lake Mendota_SenSlope_Data2.csv", header = TRUE, sep = ",")


attach(mydata)


install.packages("mblm")

library("mblm")


for (i in c(1:3)){

  x <- mblm(mydata[,i] ~ Year, mydata, repeated = FALSE)

  print(x)

}


Which gives me the following error:

Error in names(res$residuals) = as.character(1:length(res$residuals)) :

  'names' attribute [2] must be the same length as the vector [0]


Which I cannot seem to solve although as I understand it it is an error that I am causing in the mblm package. If anyone has any insight into how I could start fixing this that would be greatly appreciated!


Bailey


	[[alternative HTML version deleted]]