similar to: Passing Bablok

Hi Eik, thank you so much - it works perfectly! Thank you and best wishes Alain   Eik Vettorazzi <E.Vettorazzi@uke.de> hat am 6. Februar 2013 um 17:01 geschrieben: > Hi Alain, > here is a one-liner for a df without the rowSum column > >

Dear R-helpers, does anyone know of some package/function that can build a network from the functions that are implemented in a package, i.e. visualize the cross-references from one function to another in the same or some dependent package? An example would be a function like 'nls' on top of the hierarchy and then a network of nodes from the functions that are called within 'nls'

Dear R-helpers and graphics gurus, I have two problems with plotting confidence bands: 1. First is relatively simple. I am using the Passing-Bablok procedure to obtain "unbiased" regression coefficients. This procedure yields the "a" & "b" coefficient values along with their confidence intervals. I then plot the raw data with the regression line, but I would

Dear R-helpers, I´m in the context of writing a general function for error propagation in R. There are somehow a few questions I would like to ask (discuss), as my statistical knowledge is somewhat restricted. Below is the function I wrote, the questions are marked. Many thanks in advance. propagate <- function(expr, varList, type = c("stat", "raw"), cov = TRUE) {

Hi, here is a practical problem we would like to solve. In a pneumatic post the acceleration of the capsule is measured and plotted over time. From the graph achieved we would like to derive some kind of statistic value that describes the stress the capsule, or what is in it, is exhibited to. The amount of stress introduced to the capsule will probably depend on two things, the maximum

Hi, I''m new to XEN and have the following Problem: If I shutdown my domU the network of the dom0 leaks. A ping from another server to my xen dom0 shows a Packetloss > 80%. Then if I restart the domU after some seconds the ping show a loss from 0%. I retryed that 3 times but it is really hard to go to the dom0 via ssh and recreate the domU if the packetloss is that high. I used google

Dear all, I encounter some discrepancies when comparing the deviance of a weighted and unweigthed model with the AIC values. A general example (from 'nls'): DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) This is the unweighted fit, in the code of 'nls' one can see that 'nls' generates a vector

Dear All, I encounter some discrepancies when comparing the deviance of a weighted and unweigthed model with the AIC values. A general example (from 'nls'): DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) Now for a weighted fit: fm2DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal),

Hi Eric, you can use an annotate-layer, eg ind<-which(sig>0,arr.ind = T) ggplot(m1.melted, aes(x = Month, y = Site, fill = Concentration), autoscale = FALSE, zmin = -1 * zmax1, zmax = zmax1) + geom_tile() + coord_equal() + scale_fill_gradient2(low = "darkred", mid = "white", high = "darkblue",

Hi Steven, It's not entirely clear what you actually want to achieve in the end. As soon as you "know" x1, and assuming that the different "xi" do not differ in length in the real application, you know the length of the target vector. Instead of the loop, you can use 'Reduce' without having to initialize a starting vector. # generate sample vectors, put them in

I've been trough the R documentation for about half an hour and it's not clear to me how to do this: I need to format to character a series of integers from 1 to 1000, and I like them to look like "0001" "0002", "0059", "0123" and so on. Padded with zeroes to have four digits. Cheers! Mario. r-help-request at r-project.org wrote: > Send

Hi Anne, I would suggest to change the linear model to lm(BloodPressure~Age), as this model makes more sense in biological means (you would assume that age influences pressure, not vice versa) and also obeys the statistical assumption of weak exogeneity, that age can be measured without error, at least compared to error-prone bp measures. Cheers Am 18.04.2018 um 16:07 schrieb Gerrit Eichner:

Hello all, I'm trying to replicate a confirmatory factor analysis done in Amos. The idea is to compare a one-factor and a two-factor model. I get the following warning message when I run either model: "Could not compute QR decomposition of Hessian. Optimization probably did not converge." I have no idea what to do here. I believe posters reported the same problem. It seems

Hello, I'm trying to fit a pretty simple confirmatory factor analysis using the sem package. There's a CFA example in the examples, which is helpful, but the output for my (failing) model is hard to understand. I'd be interested in any other ways to do a CFA in R, if this proves troublesome. The CFA is replicating a 5 uncorrelated-factor structure (for those interested, it is a

I tried to bootstrap the correlation between two variables x1 and x2. The resulting distribution has two distinct peaks, how should I interprete it? The original code is attached. Y. C. Tao ---------------- library(boot); my.correl<-function(d, i) cor(d[i,1], d[i,2]) x1<-c(-2.612,-0.7859,-0.5229,-1.246,1.647,1.647,0.1811,-0.07097,0.8711,0.4323,0.1721,2.143,

Hello! I have 2 variables - predictor "pred" and response variable "DV": pred<-c(439635.053, 222925.718, 668434.755, 194242.330, 5786.321, 115537.344, 100835.368, 7133.206, 159058.286, 4079991.629, 3380078.060, 2661279.136, 2698324.478, 1245213.965, 1901815.503, 1517019.451, 1396857.736, 1034030.988, 432249.574, 342329.325, 1831335.792, 2209578.859, 1641709.754,

Hello I have recently installed R 2.14.2 on a brand new pc running Windows 64. One of things I would like to do is to recompile my R package (with C and fortran source codes) for this new environment (it works on a mac and on a linux box). To this end, I had rtools 2.14 downloaded and installed (I have also added 2.14.2/bin/x64 to PATH so that Rcmd works on a terminal). On a cmd terminal I

reading ?plotmath you might notice that "_" isn't the propper syntax for subscripts. This will work: ggtern(data=x11,aes(A,B,C,xend = c(0.7,.00,0.7),yend = c(.30,.50,.0),zend =c(.0,.50,0.3)))+ geom_point()+ theme_showarrows()+geom_segment(size=.5)+ geom_text_viewport(x=c(.45,.27,.37),y=c(.32,.29,.22),label=c("P[a]","P[b]","P[c]"), parse=TRUE)

D'oh! Thanks for pointing this out. I blame caffeine depletion at the time... -pd > On 4 Dec 2017, at 15:48 , Eik Vettorazzi <E.Vettorazzi at uke.de> wrote: > > reading ?plotmath you might notice that "_" isn't the propper syntax for > subscripts. This will work: > > ggtern(data=x11,aes(A,B,C,xend = c(0.7,.00,0.7),yend = > c(.30,.50,.0),zend

Hi, I have trouble with Sweave (I think) cutting of my newlines. As stated in the help of Sweave, I generate tex code straight from R for dynamically computed reports. If I do this in R: for (i in 0:4) {cat("\n",i,"\n")};cat("\n 3") 0 1 2 3 4 3 The output looks correct. However, Sweave for some reason seems to trim everything outside