R/writeFit.R

"writeFit" <- function (multimodel, multitheta, plotoptions) 
{
    m <- multimodel@modellist
    t <- multitheta
    res <- multimodel@fit@resultlist
    for (i in 1:length(m)) {
        x2 <- m[[i]]@x2
        x <- m[[i]]@x
        fitted <- matrix(nrow = m[[i]]@nt, ncol = m[[i]]@nl)
        irfmu <- vector()
        for (j in 1:length(x2)) {
            fitted[, j] <- res[[i]]@fitted[[j]]
            if (m[[i]]@weight) 
                fitted[, j] <- fitted[, j]/m[[i]]@weightM[, j]
        }
	if(!plotoptions@writefitivo) {	
	  write.table(fitted, file=paste(plotoptions@makeps, 
	  "_fit_dataset_", i, ".txt", sep=""),  quote=FALSE,
	  row.names = x, col.names = x2)
        }
	else {
	   cat("description line 1 \n", file = paste(plotoptions@makeps, 
	  "_fit_dataset_", i, ".txt", sep="")) 
	    cat("description line 2 \n", file = paste(plotoptions@makeps, 
	  "_fit_dataset_", i, ".txt", sep=""), append=TRUE)
	   cat("Wavelength explicit \n", file = paste(plotoptions@makeps, 
	  "_fit_dataset_", i, ".txt", sep=""), append=TRUE) 
	  cat(paste("Intervalnr", m[[i]]@nl, "\n"), 
	  file = paste(plotoptions@makeps, "_fit_dataset_", i, ".txt", 
	  sep=""), append=TRUE) 
	  write.table(fitted, file=paste(plotoptions@makeps, 
	  "_fit_dataset_", i, ".txt", sep=""),  quote=FALSE,
	  row.names = x, col.names = x2, append = TRUE)


	}

     }
}

Try the TIMP package in your browser

Any scripts or data that you put into this service are public.

TIMP documentation built on May 2, 2019, 5:55 p.m.