as.geoComp <- function(original, # Um data frame contendo Areia Silte Argila X Y perc=TRUE # Se esta ou nao em porcentagem ){ if(perc == TRUE){ malha <- transform(original,Areia=(Areia/100),Silte=(Silte/100),Argila=(Argila/100))} malha <- original malha <- transform(malha,y1=log(Areia/Argila),y2=log(Silte/Argila)) seq11 <- seq(1,length(malha$y1)*2,by=2) seq22 <- seq(2,length(malha$y2)*2,by=2) y <- c() y[seq11] <- malha$y1 y[seq22] <- malha$y2 saida <- list() saida[[1]] <- malha[c(6,7)] saida[[2]] <- data.frame(y) saida[[3]] <- data.frame(original[4],original[5]) names(saida[[1]]) <- c("Y1","Y2") names(saida[[2]]) <- c("Y") names(saida[[3]]) <- c("Coord.X","Coord.Y") class(saida) <- "geoComp" return(saida)}