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

Bailey Hewitt bailster at hotmail.com
Sun Sep 4 03:43:45 CEST 2016


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 = ",")




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

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



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!


	[[alternative HTML version deleted]]

More information about the R-help mailing list