similar to: subtraction based on two groups in a dataframe

Hi all, I used bwplot in lattice to create a 6-panel boxplot grouped by a conditioning variable (param) that displays concentration (conc) in response to treatment (trtmnt). Here is the functional part of my code followed by my three questions: library(lattice); ww<-read.csv(file="c:/Rdata/lattice_boxplot_prep.csv",header=TRUE,sep=","); attach(ww);

Full_Name: Somkiat Apipattanavis Version: 2.1.1 OS: Windows Submission from: (NULL) (128.138.44.123) Bug found in predict.locfit for density estimation # Example of bug found in prdict.locfit (Locfit) library('locfit') # generate data y =c(4281,2497,4346,5588,5593,3474,4291,2542,5195,4056, 3114,2864,4904,7625,3377,4001,4999,7191,8062,5668) x1=c( 0.258729, 1.460156, 0.192323,

Apologies for the coming to this late... 1. By now I hope Somkiat has realized that R-bugs is not the place to report problems in contributed packages. Please direct such reports to the package maintainer. 2. This is really user error. predict() expect the newdata to be a data frame containing variables with the same names as those used in the fitting process. E.g., you fitted the model with

https://bugs.freedesktop.org/show_bug.cgi?id=103033 Bug ID: 103033 Summary: nouveau can't start my discrete nvidia card (BIOS / ACPI error) Product: xorg Version: unspecified Hardware: Other OS: All Status: NEW Severity: normal Priority: medium Component:

I currently work on a draft of an aquatic bioassessment. The conditions tested are the following: ER river water T dechlorinated water control 0.5 + 0.5mg / L of malate T + 1 dechlorinated water control + 1g / L of malate T ED dechlorinated water control SED + ER + river water sediment SED ED + sediment + water dechlorinated. It is the result of AChE in muscle (fillet of fish). The production of

Hi all I ran a workload that reads & writes within 10 files each file is 256M, ie, (10 * 256M = 2.5GB total Dataset Size). I have set the ARC max size to 1 GB on etc/system file In the worse case, let us assume that the whole dataset is hot, meaning my workingset size= 2.5GB My SSD flash size = 8GB and being used for L2ARC No slog is used in the pool My File system record size = 8K ,

What do I need to do to fix this? [jonsmirl@smirl klibc-0.184]$ make make[1]: Entering directory `/home/dri/klibc-0.184/klibc' gcc -Wp,-MT,syscalls.nrs,-MD,./.syscalls.nrs.d -mregparm=3 -DREGPARM=3 -march=i386 -Os -g -falign-functions=0 -falign-jumps=0 -falign-loops=0 -nostdinc -iwithprefix include -D__KLIBC__ -DBITSIZE=32 -I../include/arch/i386 -I../include/bits32 -I../include

Dear all,I am a new user of WinBUGS and need your help. After running the following code, I got parameters of beta0 through beta4 (stats, density), but I don't know how to get the prediction of the last value of h, the variable I set to NA and want to model it using the following code.Does anyone can given me a hint? Any advice would be greatly appreciated.Best

R-help, I'm trying to evaluate the misclasification rate of an estimated discrimated function (hyDA) hyDA <- with(hyB, lda(Maturity ~ Length + Weight)) , so here is what I do (code from S-plus 6.1 Guide to Statistics Vol. 2 p 101) test <- predict(hyDA) tbl <- table(hyB$Maturity, test$class) tbl <- tbl[-2,-2] tbl1 <- cbind(tbl, error = (apply(tbl, 1, sum) - diag(tbl)) /

Hi, I'm trying to build several SCL perl516 modules for using assp spamfilter. Some modules have dependencies which spend a lot of errors or can not be build on a x86_64 system. perl-Email-Valid-0.184 for example needs Win32::IPHelper and Win32::TieRegistry. If I build this rpm for standard perl510, the Win32 modules are not required. Why are theses dependencies only for perl516? Can I

If I have a table (we'll call it, "test") containing two columns (as below): i x1 x2 x3 x4 x5 y 0 1.125 0.232 7.160 0.0859 8.905 1.5563 7 0.920 0.268 8.804 0.0865 7.388 0.8976 15 0.835 0.271 8.108 0.0852 5.348 0.7482 22 1.000 0.237 6.370 0.0838 8.056 0.7160 29 1.150 0.192 6.441 0.0821 6.960 0.3130 37 0.990 0.202 5.154 0.0792 5.690 0.3617 44 0.840 0.184 5.896 0.0812 6.932 0.1139 58

Hi there, Everytime when I paste My R output in WORD, the columns couldn't line up nicely like they appear in R console. Is there a way to fix this problem? Thanks for any help! time n.risk n.event survival std.err lower 95% CI upper 95% CI 6 21 3 0.857 0.0764 0.707 1.000 7 17 1 0.807 0.0869 0.636 0.977 10 15 1

Hi, all There are 62 samples in my data and I tested 3 times for each one, then I want to use ICC(intraclass correlation) from irr package to test the consistency among the tests. *combatexpdata_p[1:62] is the first text results and combatexpdata_p[63:124] * is the second one and *combatexpdata_p[125:186]* is the third. Here is the result:

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060711/1da40382/attachment.pl

Hello, I am trying to recursively append some data from multiple files into a common object For this, I am using in a loop NewObject <- rbind(NewObject,tempObject) For the first loop, obviously there is no NewObject ... so I wanted to do NewObject <- tempObject[0,] Now when it loops again I want to put the statement do "NewObject <- tempObject[0,]" inside a if statement

Dear all, I am attempting to explain patterns of arthropod family richness (count data) using a regression model. It seems to be able to do a pretty good job as an explanatory model (i.e. demonstrating relationships between dependent and independent variables), but it has systematic problems as a predictive model: It is biased high at low observed values of family richness and biased low at

Some help with how to re-label the vertical axis in a histogram would be appreciated. qplot(off.sc,weight=rel.freq,binwidth=.29,main="test Figure"+ylab("New from inside"))+ylab("New from outside")+ xlab("off.sc\nAggregated frequency plots for 17 equal intervals.") The code

Hello, I hesitate to call this a bug, because I could have forgotten something important, but the MarkowitzPortfolio example in fPortfolio does not work for me. Here's my code: > library(fPortfolio) > >xmpPortfolio("\nStart: Load monthly data set of returns > ") > data(berndtInvest) > # Exclude Date, Market and Interest Rate columns from data

days=Sys.Date()-1:70 price=abs(rnorm(70)) regular=rep(c(0,0,0,0,1,0,1,0,0,1),c(7,7,7,7,7,7,7,7,7,7)) y=data.frame(cbind(days,price,regular)) y is like days price regular 1 14990 0.16149463 0 2 14989 1.69519358 0 3 14988 1.57821998 0 4 14987 0.47614311 0 5 14986 0.87016180 0 6 14985 2.55679229 0 7 14984 0.89753533 0 the output I want:

Dear R helpers, I have been using the loadings function from the multiv library and I get the typical output (see below). When I try to export these results to a file using a write.table() I get the following error message "Error in as.data.frame.default(x[[i]], optional = TRUE) : can't coerce loadings into a data.frame" Any idea why write.table is doing that and any