Modification des valeurs dans le fichier txt à l'aide de R

Question

Je souhaite écrire un code logiciel qui modifiera les valeurs des paramètres requis de mon fichier txt d'entrée de modèle hydrologique. certains paramètres sont fixes et certains je vais changer aux valeurs observées. Par exemple HYDRUS_Version = 4 est fixé et WaterFlow = 1 et SoluteTransport = 0 etc sont les paramètres que je veux changer leurs valeurs. Je veux assigner les valeurs de waterFlow = 5 ou WaterFlow = 3.1 et SoluteTransport = 2 ou = 2.2

J'ai essayé ce code pour l'écoulement d'eau mais les valeurs dans mes dossiers de txt ne sont pas changées. Je suis tout nouvel apprenant de R.

lines <- readLines("G:/Rlearning/HYDRUS1D.txt") 

     library(gsubfn) 
     i1 <- grepl("Vertical Conductivity", lines) 
     lines[i1] <- gsubfn("[0-9.]+", ~format(as.numeric(x)*2, 
              scientific = FALSE), lines[i1]) 

Voici le dossier complet txt.

; 
[Main] 
HYDRUS_Version=4 
WaterFlow=1 
SoluteTransport=0 
Unsatchem=0 
Unsatchem=0 
HP1=0 
HeatTransport=0 
EquilibriumAdsorption=1 
MobileImmobile=0 
RootWaterUptake=1 
RootGrowth=0 
MaterialNumbers=1 
SubregionNumbers=1 
SpaceUnit=cm 
TimeUnit=days 
PrintTimes=160 
NumberOfSolutes=0 
InitialCondition=1 
; 
[Profile] 
NumberOfNodes=101 
ProfileDepth=1.2E+02 
ObservationNodes=5 
GridVisible=1 
SnapToGrid=1 
ProfileWidth=80 
LeftMargin=40 
GridOrgX=0 
GridOrgY=0 
GridDX=5.E+00 
GridDY=5.E+00 

Answer 1

Créer une fonction personnalisée à chaîne correspondant à rechercher et remplacer la valeur:

replaceFn <- function(phrase, value, file){ 
    line <- file[grep(phrase, file)] 
    value.pos <- as.numeric(gregexpr("=", line)) 
    file[grep(phrase, file)] <- paste0(substring(line, 1, value.pos), value) 
    return(file) 
} 


dat <- replaceFn("WaterFlow", 3.1, dat) 
dat <- replaceFn("SoluteTransport", 2.2, dat) 

[1] ";"      "[Main]"     "HYDRUS_Version=4"  "WaterFlow=3.1"   "SoluteTransport=2.2"  "Unsatchem=0"    
[7] "Unsatchem=0"    "HP1=0"     "HeatTransport=0"   "EquilibriumAdsorption=1" "MobileImmobile=0"  "RootWaterUptake=1"  
[13] "RootGrowth=0"   "MaterialNumbers=1"  "SubregionNumbers=1"  "SpaceUnit=cm"   "TimeUnit=days"   "PrintTimes=160"   
[19] "NumberOfSolutes=0"  "InitialCondition=1"  ";"      "[Profile]"    "NumberOfNodes=101"  "ProfileDepth=1.2E+02" 
[25] "ObservationNodes=5"  "GridVisible=1"   "SnapToGrid=1"   "ProfileWidth=80"   "LeftMargin=40"   "GridOrgX=0"    
[31] "GridOrgY=0"    "GridDX=5.E+00"   "GridDY=5.E+00" 

Données d'entrée:

dat <- c(";", "[Main]", "HYDRUS_Version=4", "WaterFlow=1", "SoluteTransport=0", 
"Unsatchem=0", "Unsatchem=0", "HP1=0", "HeatTransport=0", "EquilibriumAdsorption=1", 
"MobileImmobile=0", "RootWaterUptake=1", "RootGrowth=0", "MaterialNumbers=1", 
"SubregionNumbers=1", "SpaceUnit=cm", "TimeUnit=days", "PrintTimes=160", 
"NumberOfSolutes=0", "InitialCondition=1", ";", "[Profile]", 
"NumberOfNodes=101", "ProfileDepth=1.2E+02", "ObservationNodes=5", 
"GridVisible=1", "SnapToGrid=1", "ProfileWidth=80", "LeftMargin=40", 
"GridOrgX=0", "GridOrgY=0", "GridDX=5.E+00", "GridDY=5.E+00")