similar to: contour figure from rsm

Dear all I try to persuade if statement in my function to work as I want (but I am not very successful:( My question is why my function with if statement is evaluated correctly in case of atomic variables and incorrectly in case of vector variables. Here is an example: #function with if statement fff <- function(otac,sklon) { test <- (otac *

Sorry everybody, temporary dyslexia. Sent from my BlackBerry wireless smartphone

Hola compañeros de la lista. Los molesto con la siguiente duda. En un diseño central compuesto (CCD) con dos factores (V1 y V2) y una variable de respuesta (R), utilizando valores codificados (-1.4142, -1, 0, 1, 1.4182), al aplicar la orden: rsm.segundo.orden <- rsm(R ~ Bloque + SO(V1, V2), data = DATOS.Codificados) Obtengo el siguiente modelo: R = 103.92 -2.16

Hello R experts I am looking for: - advice - sample commands - pointers in MASS 4E - packages That help to produce 1D and 2D response plots from quadratic RSM linear regression models that have been produced in R. Let me explain a little what I mean. For 1D plots I mean the variable range (either scaled or actual) on the X axis and the Y axis showing the model prediction (preferably with

Hi all (and especially Frank), I'm trying to use x=T, y=T in order to run a validated stepwise cox regression in rsm, having multiply imputed using mice. I'm coding model.max<-fit.mult.impute(baseform,cph,miced2,dated.sexrisk2,x=T,y=T) baseform is baseform<-Surv(si.age,si=="Yes")~ peer.press + copy.press + excited + worried + intimate.friend + am.pill.times +

Hello, I have a data.frame called spec containing data about samples. But I don't have these data for all my samples. So if I have data (that is code of the sample is in spec$Code), I would like to assign data1 to the variable m. If I don't have this data, I would like to assign 1 to m. I tried this : m<-spec$data1[spec$Code==code]*(code %in%specmodules$Code) + 1*(!code %in%

Hola compañeros de la lista. Estoy aprendiendo a usar el paquete "rsm" para superficies de respuesta. Siguiendo este ejemplo todo va bien, hasta que trato de obtener el valor máximo predicho para la variable de respuesta en los valores de x1 y x2 estimados. ------------------------------------------------------------- library("rsm") ChemReact CR <- coded.data(ChemReact,

Hello all, I apologize for a terribly simple question. I'm used to using Stata and am trying to `switch' over to R. I would like to recode a vector in a data frame when the value of it meets the following condition: if (dataframe$factor=='3'){dataframe$thevector<-(an arithmetic mean). What I would like to result is the creation of a new variable within the data frame for

I'm trying to make a loopo with some file.... > > library(dplR) > files.rwl <- system("ls *.rwl", intern=TRUE) > files.rwl [1] "cimfasy.rwl" "rocquce.rwl" > for (i in files.rwl) assign(gsub("\\.rwl$", "_rwl", i), read.rwl(i, header=0)) There are 70 series There are 21 series > > cimfasy_rwl 1990 1.55 0.00

Hello I have a data set that looks like this; > b2 dato chr status PRRSvac PRRSsanVac PRRSsanDk PRRSdk 33 2007-12-03 090432 R?d SPF 34 2007-02-09 090432 R?d SPF+sanDK 35 2002-12-17 090432 R?d SPF+DK 36 2002-11-27 090432 R?d SPF+sanDK 37 2002-07-23

Attached is a snapshot of the new version of the FAQ. What is still missing is something on eval and .Options versus options(). As always, feedback is greatly appreciated. Best, -k *********************************************************************** R FAQ Kurt Hornik v0.2-0, 1997/09/01 This document contains answers to some of the most frequently asked questions about R. Feedback

An updated version of the R FAQ to accompany the new 0.60 release is now available at the usual site, http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html A plain text version of the FAQ is appended below. -kh ****** snip snip snip ************************************************** R FAQ Kurt Hornik v0.60-6, 1997/12/08 This document contains answers to some of the most frequently asked

An updated version of the R FAQ to accompany the new 0.60 release is now available at the usual site, http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html A plain text version of the FAQ is appended below. -kh ****** snip snip snip ************************************************** R FAQ Kurt Hornik v0.60-6, 1997/12/08 This document contains answers to some of the most frequently asked

An updated version of the R FAQ to accompany the new 0.60 release is now available at the usual site, http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html A plain text version of the FAQ is appended below. -kh ****** snip snip snip ************************************************** R FAQ Kurt Hornik v0.60-6, 1997/12/08 This document contains answers to some of the most frequently asked

Estimado Gilsanz Algo que puede ser, ¿tiene instalado tcltk en Windows? Yo ante un problema parecido lo solucioné realizando un contador, en mi caso conozco la cantidad de consultas, entonces simplemente ?imprimo? en pantalla 1 de 100, 2 de 100, etc. No es gráfico, pero R me informa dos cosas, que está trabajando y cuánto falta, feo a la vista pero muy simple. Javier Rubén Marcuzzi De:

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems