search for: 0.006

I've gone over the data and do not see my error; the dput() output of the data frame and the pdf output of cenboxplot() are attached. The command used: cenboxplot(sb.t$quant, sb.t$ceneq1, range=1.5, main='Total Recoverable Antimony', xlab='Pre-Mining Era', ylab='Concentration (log mg/L)') (on a single line in emacs). The RL on the plot is drawn at 0.01 rather

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

I can't get my simultaneous equations to work using system fit. Please help. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) str(Empdata) Empdata$gnipc<-as.numeric(Empdata$gnipc) install.packages("systemfit") library("systemfit") pdata <- plm.data(Empdata,

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip

Hi, I want to produce two parallel rotated strips with color gradient. So far, the sample strip is something produced by this: pushViewport(viewport(x = unit(0.638, "npc"), y =unit(0.386, "npc"), width=.62, height=0.006, angle=137.2)) grid.rect(y=100:1/100, just="top", gp=gpar(col=NA, fill=colorRampPalette(c("lightgray",

Hi R user I need to read in some values from a computer program output. I can't change the output format because the developer of the program doesn't allow to change the format of output. There are two formats. First one looks like this if I have 10 variables, ------------------------------------------------------------------------------------------------------ [ 1]

Folks, The contour() function wants x and y to be in increasing order. I have a situation where I have a grid in x and y, and associated z values, which looks like this: x y z [1,] 0.00 20 1.000 [2,] 0.00 30 1.000 [3,] 0.00 40 1.000 [4,] 0.00 50 1.000 [5,] 0.00 60 1.000 [6,] 0.00 70 1.000 [7,] 0.00 80 0.000 [8,] 0.00 90

I having problem with dendrogram leaf names when I read a tab delimited file into dataframe; I have a text file, tab delimited, using read.table into a data frame as follows: > test1<-read.table("c:\\R\\data\\Tremont4.txt", header=TRUE, sep="\t") When I do this the "test1" data frame is picking up my first column names as part of the data and not the case

Dear R-help reader, it would be great if someone knows what I'm doing wrong. I have (shorten) dataframe, which consists of a group identification and a number >ex UID REL 1 R1.B8.31 0.000 2 R1.B8.31 0.000 3 R1.B8.31 0.000 4 R1.B8.31 0.000 5 R1.B8.38 0.010 6 R1.B8.38 0.060 7 R1.B8.38 0.006 8 R1.B8.38 0.010 9 R1.B8.48 0.080 10 R1.B8.48 NA 11 R1.B8.48 0.006 I'm

Hi all, I'm an R newbie, I did this script to create a scatterplot using the "tree" matrix from "datasets" package: library('datasets') with(trees, { plot(Height, Volume, pch=3, xlab="Height", ylab="Volume") symbols(Height, Volume, circles=Girth/12, fg="grey", inches=FALSE, add=FALSE) } ) I'd like to use the column Named

Folks, A few days ago, I had asked a question on this mailing list about making a contour plot where a function z(x,y) is evaluated on a grid of (x,y) points, and the data structure at hand is a simple table of (x,y,z) points. As usual, R has wonderful resources (and subtle complexity) in doing this, and the gurus of the list showed me the way. Here's a complete working example. One might

Hola Colegas: Me reintegro al estudio de R y quisiera me pudieran resolver el problema que planteo en el archivo adjunto. Gracias anticipadas. *MANOLO MÁRQUEZ P.* ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130511/b0c9f15d/attachment.html> ------------ próxima parte ------------ An embedded

Recently there was a question on using lme with large data sets. As an experiment I fit a linear mixed-effects model to a data set with about 350,000 observations on 6 predictors, a numerical response, and a single grouping factor. The timings shown below were on a 1.2 GHz Athlon with 1 GB of PC133 memory and 2 GB of swap. The operating system is Debian 3.0 GNU/Linux. The kernel is 2.4.14.

I have a data frame with 1 factor, one date, and 37 numeric values: str(waterchem) 'data.frame': 3525 obs. of 39 variables: site : Factor w/ 64 levels "D-1","D-2","D-3",..: 1 1 1 1 1 ... $ sampdate : Date, format: "2007-12-12" "2008-03-15" ... $ CO3 : num 1 1 6.7 1 1 1 1 1 1 1 ... $ HCO3 : num 231 228 118 246

Hi list, I have the following code to compute the acf of a time series acfresid <- acf(residfit), where residfit is the series when I type acfresid at the prompt the follwoing is displayed Autocorrelations of series ?residfit?, by lag 0.0000 0.0833 0.1667 0.2500 0.3333 0.4167 0.5000 0.5833 0.6667 0.7500 0.8333 1.000 -0.015 0.010 0.099 0.048 -0.014 -0.039 -0.019 0.040 0.018

I've data that look like: Outcome V1_AA V1_EU V1_NA V2_AA V2_EU V2_NA 0 0.046 1.001 0.954 0.045 1.001 0.954 0 0.007 1 0.993 0.007 1 0.993 1 1.774 0.217 0.009 1.774 0.217 0.009 1 0.004 1.996 0 0.004 1.996 0 1

Consider the following little "benchmark" > require(Matrix) > tmp <- Matrix(c(rep(1,1000),rep(0,9000)),ncol=1) > ind <- sample(1:10000,10000) > system.time(tmp[ind,]) user system elapsed 0.004 0.001 0.005 > ind <- sample(1:1000,10000,replace=TRUE) > system.time(tmp[ind,]) user system elapsed 0.654 0.006 0.703 >

I have the following small dataset: > stardata NS HE EB CW RW PW 1 0 0.000 0.000 0.042 0.006 0 2 0 0.006 0.000 0.013 0.005 0 3 0 0.000 0.011 0.000 0.000 0 I have plotted the star diagrams as follows: stars(stardata, key.labels = dimnames(stardata)[[2]], labels = NULL, key.loc = NULL, draw.segments=TRUE, col.segments="gray", lty="blank") I am