similar to: 2^2 problem revisited

Hi, I am very new to R and I've been trying to work through the R book to gain a better idea of the code (which is also completely new to me). Initially I imputed my data from a text file and that seemed to work ok, but I'm trying to examine linear relationships between gdist and gair, gdist and gsub, m6dist and m6air, etc. This didn't work and I think it might have something to do

Greetings, I have a data set look like this: > print(MyResult) 2 5 9 10 11 13 mean 6.108394 3.272211e+01 3.951703e+02 2.086217e+03 79.838867 2.592272e+04 median 1.288889 2.743511e+00 4.061491e+00 4.943710e+01 -2.821667 -3.040625e+00 var 2885.350595 2.343144e+05 1.022680e+06 1.200589e+07 365179.508831

Hi all, I have been trying to reproduce an analysis from Douglas Montgomery?s book on design and analysis of experiments. Table 6.10 of example 6.2 on page 246, gives a table as follows: > NPK <- expand.grid(A=mp,B=mp,C=mp,D=mp) > Rate <- c(45,71,48,65,68,60,80,65,43,100,45,104,75,86,70,96) > filtration <- cbind(NPK,Rate) > filtration A B C D Rate 1 - - - - 45 2

Hello, I am getting output from anova() and summary(aov()) that depends on the order of the factors in the fitted model object, and this has me baffled. I see this dependency with the data.frame below but not with an example (table 6.4) from Montgomery's DOE book. This is with R 1.3.0 on Debian GNU-Linux. Where have I gone wrong? > centerpts run sample CH50mg 1 day1 dev126 0.56 2

Hello, I have a file with very big values. I want to display the real values with classic ways (without exposant) > summary(a[,1]); Min. 1st Qu. Median Mean 3rd Qu. Max. 1.198e+09 1.199e+09 1.200e+09 1.200e+09 1.201e+09 1.202e+09 Can you help please? Regards Adel -- PhD candidate in Computer Science Address 3 avenue lamine, cité ezzahra, Sousse 4000 Tunisia tel:

I am scraping data from a web page using XML (excellent package BTW - that's scraping data the easy way!). So far, I've got the code: tables <- readHTMLTable(theurl) rhf <- tables$tabResHistFull div1 <- rhf[which(rhf$V1=="Div ps"),] div1 which is giving me the result:        V1 V2    V3    V4    V5    V6    V7          V8    V9   V10   V11   V12   V13   V14  V15 15

I have some doubts about the validity of my procedure to estimeate linear contrasts ina a factorial design. For sake of semplicity, let's imagine a one way ANOVA with three levels. I am interested to test the significance of the difference between the first and third level (called here contrast C1) and between the first and the seconda level (called here contrast C2). I used the following

sink() isn't behaving as i expect, when used inside a function, eg: x<-data.frame(F=c("O","O")) f<-"foo.txt" sink(f); format(x); sink(); # foo.txt looks great! foo<-function(x,f) { sink(f); format(x); sink(); } foo(x,f=f) # foo.txt is empty! why is this, and how can i successfully sink() within a function? my real function does some

sink() isn't behaving as i expect, when used inside a function, eg: x<-data.frame(F=c("O","O")) f<-"foo.txt" sink(f); format(x); sink(); # foo.txt looks great! foo<-function(x,f) { sink(f); format(x); sink(); } foo(x,f=f) # foo.txt is empty! why is this, and how can i successfully sink() within a function? my real function does some

I want to fit discrete data that was measured on a wavegenerator. In this minimal example i generate some artificial data: testsin <- 2+ 5 * sin(1:100) #generate sin data testsin <- testsin+ rnorm(length(testsin), sd = 0.01) #add noise mydata <- list(X=1:100, Y=testsin) # generate mydata object nlmod <- nls(X ~ A+B*sin(C* Y), data=mydata, start=list(A=2, B=4, C=1), trace=TRUE) #

Colleagues, I am fitting an Emax model using nls. The code is: START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL <- list(maxiter=1000, warnOnly=T) #FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula FORMULA <- as.formula(YVAR ~ EMAX - EFFECT / (1 + (C50/XVAR)^GAMMA)) FIT <-

I have been having the same problem as poster Hodgess, below. It appears that her question was never answered, so I would like to share a solution with the community. The problem is the (apparent?) inability to produce moving range process behavior (a.k.a. "control") charts with individuals data in the package "qcc" (v. 2.0). I have also struggled with the same limitation in

I have a set of modeled climate data recorded at irregular intervals. The format of the data is such that there are monthly measurements for the years 2000, 2020, 2050, 2080, etc. Therefore I have 12 regular records, a skip of some number of years, then 12 more monthly records, another skip, and so on.... I created a dataframe from .txt with the read.table() command. For starters I need

Hi R-users,   I have this code below and use some code from sn package for multivariate skew t (rmst) to generate a set of random numbers.   I try to do a looping to produce ‘rand’ for different values of degrees of freedom (degf).  Actually, I’m not sure how to specify the attribute for the rand.  As we know, rand will give us a set of results + all the parameters used, so just wonder whether we

hello, I'd like know how to sort a data frame in R for example how I should do to sort by Catholic with swiss data frame like below thanks Fertility Agriculture Examination Education Catholic Infant.Mortality Courtelary 80.2 17.0 15 12 9.96 22.2 Delemont 83.1 45.1 6 9 84.84 22.2

##I want to remove the rows where the row sums are zero and this is as far as I have gotten ffg <- (structure(list(CD = c(0, 0, 0, 0, 3.125, 0, 0, 0, 0, 1.6, 3.125, 0, 0, 6.25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.125, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.125, 0, 0, 0, 0, 0, 0, 0, 0,

Estimada comunidad tengo un problema del que no encuentro datos que me ayuden mucho en la web. Estoy haciendo calculos por grupos con data,table. Tengo un archivo (zp.res) con tres columnas que clasifican los datos (sol, con, dia) y una columna de datos numericos (media), de la siguiente forma: sol con dia media 1: con 0 1 -22.6 2: con 0 1 -36.6 3: con 0 1 -35.6 y

All, I would like to apply a set of rules to each row of the sample data set below. The rule sets are the guidelines for determining an individual's date for retirement eligibility. The rules are found in this document, http://www.opm.gov/feddata/RetirementPaperFinal_v4.pdf. I am only interested in the top two categories for retirement eligibility, the CSRS and FERS plans. The data set has

Dear R users, I seek for a more elegant manner to manipulate my data that that I produced so far. Data is in a data frame THC515.DATA -sample at the end of this mail- of 5 variables, a string and 4 continuous numeric variables. I need to get the mean of two variables "Q" and "T" at the levels of a third one "dP", individually for each level of "n". I

Hello, I'm trying to get an energenie EG-UPS-033 (http://energenie.com/item.aspx?id=8087) to work with nut. It standard comes with linux software (UPSmart), which shows the input/output voltage, load level, input frequency. and battery and case temperature. The values change about every second. The selected monitor mode in this tool is Mega(USB). The bad thing of this tool that it