similar to: Repeated cross-validation for a lm object

Dear R users I'd like to know if it is possible to create conditional plots in the soiltexture package. I have data from the particle size distribution of soils derived from two parent materials and I want to use the parent material as conditioner. Thanks a lot. Alessandro Samuel-Rosa Post-Graduation Program in Soil Science Federal University of Santa Maria Av. Roraima, n? 1000, CEP

Hi everybody, I'm trying, without success, to avoid some variables labels to overlap on a capscale (vegan package) object graph. I was able to move the labels with the coding below, but it moves the entire vector... CAPobj<-CAPpotiFT for(i in c(2:8))CAPobj$CCA[[i]][,1]<-CAPpotiFTI$CCA[[i]][,1] CAPpotiFTI$CCA[[i]][,1][2]<-0.85 And here the graphic script

Hi! I have to import some TXT files into R, but the separation between the columns are made with different blank spaces, but each file use the same separation. Example: 31 104 5 0 11RUA SAO SEBASTIAO 25 BAIRRO FILETO 01

Hi, I compare results of a simple two-level poisson estimated using lmer and those estimated using MLwiN and Stata (v.9). In R, I trype: ------------------------------------------------------------------------------------------- m2 <- lmer(.D ~ offset(log(.Y)) + (1|pcid2) + educy + agri, male, poisson) -------------------------------------------------------------------------------------------

?s 22:54 de 11/06/2023, javad bayat escreveu: > Dear Rui; > Many thanks for your email. I used one of your codes, > "data2$LU[which(data2$Layer == "Level 12")] <- "Park"", and it works > correctly for me. > Actually I need to expand the codes so as to consider all "Levels" in the > "Layer" column. There are more than hundred

Dear Rui; Many thanks for your email. I used one of your codes, "data2$LU[which(data2$Layer == "Level 12")] <- "Park"", and it works correctly for me. Actually I need to expand the codes so as to consider all "Levels" in the "Layer" column. There are more than hundred levels in the Layer column. If I use your provided code, I have to write it

Hello everyone, I am working with plm package and I have problem with random and within models, which are giving errors which says "empty model". However, the model is not empty. In the source code for plm.fit, where the error originates it says something like (writing from the top of my head...) X <- model.matrix(formula,data, lhs=1,...) if (ncol(X) == 0) stop("empty

Dear all, I intend to do some statistics with my raster table, like correlations or regression analysis. For that, I need that each pixel be in one line, that is, there is no count > 1. So, all columns should have the same number of lines, when I compare (combine) rasters. I'd like to know if there is a way to create lines and copy the values of the lines, when count is > 1. So, even

I have this partial: 1 - bairros = bairros_select 2 - selected_id ||= nil 3 - if bairros.any? 4 - cid = nil 5 = select_tag :imovel, :bairro_id do 6 - for b in bairros 7 - sel = b.id == selected_id ? {:selected => true} : {} 8 - if cid != b.cidade_id 9 %optgroup= b.cidade.cidade 10 - cid = b.cidade_id 11 %option{sel, :value =>

Hi, I am looking to do some quantitative genetic analyses using animal models and was wondering if someone could suggest an appropriate package in R. It would help if it was similar to the ASReml genetic analyses software. Thanks, Deepa Senapathi Deepa Senapathi Centre for Agri-Environmental Research (CAER) School of Agriculture, Policy & Developement. University of Reading

Dear R community, I have a question concerning recoding of a variable. I have a data set in which there is a variable devoted to the ISCO code describing the occupation of this certain individual (http://www.ilo.org/public/english/bureau/stat/isco/isco88/major.htm). Every type of occupation begins with a number and every number added to this number describes th occupation more detailed. Now my

Desde el terminal de Ubuntu, he conseguido descargar los datos usando esto: DESCARGAR LISTA DE ESTACIONES wget --method GET --header 'cache-control: no-cache' --no-check-certificate --output-document - https://opendata.aemet.es/opendata/api/valores/climatologicos/inventarioestaciones/todasestaciones/?api_key=*vuestra api key* DESCARGAR DATOS HISTÓRICOS Diciembre 2017, ALICANTE wget

?s 21:05 de 11/06/2023, javad bayat escreveu: > Dear R users; > I am trying to fill a column based on a specific value in another column of > a dataframe, but it seems there is a problem with the codes! > The "Layer" and the "LU" are two different columns of the dataframe. > How can I fix this? > Sincerely > > > for (i in 1:nrow(data2$Layer)){ >

Estimad en s errer en s He hecho este modelo glm m1.pile<-glm(ger~tem+pot+time+I(tem^2)+I(tem^2):pot ,family="binomial" ,data=long.PILE ) Que nos da la probabilidad de germinación de una semilla en función de tem (Temperatura), pot (Humedad del suelo) y time (Tiempo que la semilla pasa en esas condiciones). Ahora quiero, para diferentes tem, pot

Hi everyone, I'm aiming to run an external executable (say filetorun.EXE) in parallel. The external executable collect needed data from a file, say "input.txt" and, in turn,generates several output files, say "output.txt". I need to generate "input.txt", run the executable and keep "input.txt" and "output.txt". I'm using Windows 7, R

Hola eRRer en s, He creado un modelo mixto (adjunto) con lme4 He dibujado la predicción usando: df <- ggpredict(m.glob.pre.anu , terms = c("pre", "area_m2", "sum.vs.win") ) plot(df) Luego he arreglado el gráfico usando: df <- ggpredict(m.glob.pre.anu , terms = c("pre", "area_m2",

Hi, I type a value, for example, "10,00". I see in the screen "10,00" but when print that value is multiplied by 100x resulting: "1000,00". What can I do? wine 0.9.55-1 Debian/Sid Thanks. -- S?vio M Ramos Arquiteto, Rio, RJ S? uso Linux desde 2000 www.debian.org

Hi R-users I recently downloaded RdbiPgSQL 1.8.0 and Rdbi 1.8.0 from Bioconductor to be installed under R 2.4.1. When requiring RdbiPgSQL an error message is showed as follows: > require(RdbiPgSQL) Loading required package: RdbiPgSQL Error in library(package, lib.loc = lib.loc, character.only = TRUE, logical = TRUE, : ‘RdbiPgSQL’ is not a valid package – installed <

Dear fellows, I need a R script that compare known curves (e.g. logistic, exponential) with my curve. That curve was generated fitting data of forest cover (hectares) measured in 10 road distances (buffers). I´d like that comparison should be done using AICc to select the best model, that is, the best fitted curve in comparison with my data. Anyone could help me? Best regards Simone.

Un ejemplo con un modelo más simple: He especificado este modelo: >formula(m2.pile) ger ~ tem + pot + time Si hago predict me da: >predict(m2.pile,newdata=data.frame(tem=25,pot=0,time=3),type="response") 0.08243262 Extraigo los coeficientes: > coef(m2.pile) (Intercept) tem pot time -1.89521331 -0.02303313 4.74499714 0.02043222 Ahora calculo la