search for: 0.0048

Dear Guys, I do converting codes from Fortran into R and got stuck in solving LOOPING procedure with R. In FORTRAN, it is (DO and END DO) for looping in the net. In R, it is (FOR with { }). I believe there is something wrong with my coding in R, do hope that you can help me solving following problems. It seems easy, but results are not the same.

Dear R-helpers, I am trying to plot a Bland-Altman-Plot using the BA.plot function from the package MethComp. While there is a function to transform the values for analysis as shown in the snippet below, I would like to have logarithmic axes for display as well. The usual log = 'xy' does not work because of the properties of the y-asxis (positive and negative values). I am sure that

Dear All,   wondering if someine can access the link to the randsamp code referenced in the R-help archive here: http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg75645.html ? I have tried but for whatever reason I can not get trough. My problem seems to be similar to what the author originally described there, so having access to it would be great. Else, if you have any thougths on sampling

Hello: Will someone PLEASE help me with this problem. This is the third time I've posted it. When I appply anova() to two equations estimated using glmmPQL, I get a complaint, > anova(fm1, fm2) Error in anova.lme(fm1, fm2) : Objects must inherit from classes "gls", "gnls" "lm","lmList", "lme","nlme","nlsList", or

I am trying to fit a generalised least squares model using gls in the nlme package. The model seems to fit very well when I plot the fitted values against the original values, and the model parameters have quite narrow confidence intervals (all are significant at p<5%). The problem is that the log likelihood is always given as -Inf. This doesn't seem to make sense because the model

Hello: Let me rephrase my question to attract interest in the problem I'm having. When I appply anova() to two equations estimated using glmmPQL, I get a complaint, > anova(fm1, fm2) Error in anova.lme(fm1, fm2) : Objects must inherit from classes "gls", "gnls" "lm","lmList", "lme","nlme","nlsList", or "nls"

Hi, I recently changed from uw imap to dovecot on the sound recommendation of a friend and have mostly succeeded in getting all of my clients up and running, but am really stuck with the iPhone which is failing to make connections. I run certificates on all of my clients and thunderbird happily connects both locally and remotely. I installed the certificate on the iPhone after great pain (pk12

On Fri, Sep 14, 2007 at 11:42:18PM -0700, Tanya Lattner wrote: > The 2.1 pre-release (version 1) is available for testing: > http://llvm.org/prereleases/2.1/version1/ > > [...] > > 2) Download llvm-2.1, llvm-test-2.1, and the llvm-gcc4.0 source. > Compile everything. Run "make check" and the full llvm-test suite > (make TEST=nightly report). > > Send

I am plotting three Pearson Type IV distributions. It looks like I have to plot the distribution with the highest value of y and then use lines() to add the two distributions that are shorter / have lower max values of y. The following code figures out which distribution has the max y value, plots it first and then uses lines for the other two distributions with a series of three if statements.

hello all, I wonder if anyone could give me a hint on which statistical technique I should use and how to carry it out in R in my case. Thanks in advance. My data is composed of two columns, the same numerical variable (continuous) from actual measurement and model prediction. My objective is to compare the data agreement (if there is significant difference) and make conclusions about the model

Hello: I am getting an error message when appplying anova() to two equations estimated using glmmPQL. I did look through the archives but didn't finding anything relevant to my problem. The R-code and results follow. Hope someone can help. ANDREW ____________________________ > fm1 <- glmmPQL(choice ~ day + stereotypy, + random = ~ 1 | bear, data = learning, family =

Hello: I am really stuck on this problem. Why do I get an error message with anova() when I compare these two equations? Hope someone can help. ANDREW ____________________________ > fm1 <- glmmPQL(choice ~ day + stereotypy, + random = ~ 1 | bear, data = learning, family = binomial) > fm2 <- glmmPQL(choice ~ day + envir + stereotypy, + random = ~ 1 |

I am trying to find the percentage of the parameters explaining the bacterial community composition. I have one data matrix with relative abundance of OTUs and one with environmental parameters. I used varpart in vegan package but the values in the venn diagram is bigger than 100% in total.How is it possible? What might be the reason? Thank you library(vegan) gotud <-

Hello list, I'm trying to generate classifiers for a certain task using several methods, one of them being decision trees. The doubts come when I want to estimate the cross-validation error of the generated tree: tree <- rpart(y~., data=data.frame(xsel, y), cp=0.00001) ptree <- prune(tree, cp=tree$cptable[which.min(tree$cptable[,"xerror"]),"CP"]) ptree$cptable

#is there a way to get NA in the table of descriptive statistics instead of the function stopping Thank you in advance #data x.f <- structure(list(Site = structure(c(9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L), .Label = c("BC", "HC", "RM119", "RM148", "RM179", "RM185",

Dear all, I'm using R 2.8.1 under Windows Vista on a dual core 2,4 GhZ with 4 GB of RAM. I'm trying to reproduce a result out of "Analysis of Financial Time Series" by Ruey Tsay. In R I'm using the fGarch library. After fitting a ar(3)-garch(1,1)-model > model<-garchFit(~arma(3,0)+garch(1,1), analyse) I'm saving the results via > result<-model

On Mon, Dec 01, 2014 at 06:17:03PM +0800, Jason Wang wrote: > Hello: > > We used to orphan packets before transmission for virtio-net. This breaks > socket accounting and can lead serveral functions won't work, e.g: > > - Byte Queue Limit depends on tx completion nofication to work. > - Packet Generator depends on tx completion nofication for the last > transmitted

On Mon, Dec 01, 2014 at 06:17:03PM +0800, Jason Wang wrote: > Hello: > > We used to orphan packets before transmission for virtio-net. This breaks > socket accounting and can lead serveral functions won't work, e.g: > > - Byte Queue Limit depends on tx completion nofication to work. > - Packet Generator depends on tx completion nofication for the last > transmitted

On Mon, Dec 1, 2014 at 6:42 PM, Michael S. Tsirkin <mst at redhat.com> wrote: > On Mon, Dec 01, 2014 at 06:17:03PM +0800, Jason Wang wrote: >> Hello: >> >> We used to orphan packets before transmission for virtio-net. This >> breaks >> socket accounting and can lead serveral functions won't work, e.g: >> >> - Byte Queue Limit depends

On Mon, Dec 1, 2014 at 6:42 PM, Michael S. Tsirkin <mst at redhat.com> wrote: > On Mon, Dec 01, 2014 at 06:17:03PM +0800, Jason Wang wrote: >> Hello: >> >> We used to orphan packets before transmission for virtio-net. This >> breaks >> socket accounting and can lead serveral functions won't work, e.g: >> >> - Byte Queue Limit depends