similar to: Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours

Displaying 20 results from an estimated 900 matches similar to: "Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours"

2017 Jul 31
2
Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours
> Hi Rosa > something like > plot(1,1, sub=expression(lambda^"2")) > So with your example, do you want something like > plot(c(1:5), CRP7raw[1,], type = "n", xlim=c(1,5), ylim=c(-10,5) , > xlab="Day in ICU", > ylab="CRP (mg/dL)", > sub = mtext(expression(lambda^2))) OOps! Either plot( ..., sub = *) or
2017 Jul 31
0
Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours
Hi Rosa something like plot(1,1, sub=expression(lambda^"2")) So with your example, do you want something like plot(c(1:5), CRP7raw[1,], type = "n", xlim=c(1,5), ylim=c(-10,5) , xlab="Day in ICU", ylab="CRP (mg/dL)", sub = mtext(expression(lambda^2))) CRP7graph <- apply(CRP7, 1, lines, col="gray") Cheers Petr >
2017 Jul 31
0
Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours
Hi Martin see in line > -----Original Message----- > From: Martin Maechler [mailto:maechler at stat.math.ethz.ch] > Sent: Monday, July 31, 2017 10:52 AM > To: PIKAL Petr <petr.pikal at precheza.cz> > Cc: Rosa Oliveira <rosita21 at gmail.com>; r-help mailing list <r-help at r- > project.org> > Subject: Re: [R] Superscript and subscrib R for legend x-axis and
2017 Jul 31
4
Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours
>>>>> PIKAL Petr <petr.pikal at precheza.cz> >>>>> on Mon, 31 Jul 2017 09:11:18 +0000 writes: > Hi Martin see in line >> -----Original Message----- From: Martin Maechler >> [mailto:maechler at stat.math.ethz.ch] Sent: Monday, July >> 31, 2017 10:52 AM To: PIKAL Petr <petr.pikal at precheza.cz> >> Cc:
2017 Jul 31
0
Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours
Hi, everyone, Before everything, thanks. Lots of thanks ;)!!!! I don?t think you understood everything I need to do. I want to write t_i instead of "Day in ICU? [i subscript for t] and y_ij instead of "CRP (mg/dL)? [ij superscript for y]. The label of the axis? :( Can you help me on that task? Thanks!!!!! Best, Rosa Oliveira > On 31 Jul 2017, at 10:28, Martin Maechler
2017 Aug 01
0
Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours
Hi Keep your messages coppied to R helplist, others could give you answers too. See in line From: Rosa Oliveira [mailto:rosita21 at gmail.com] Sent: Tuesday, August 1, 2017 4:38 PM To: PIKAL Petr <petr.pikal at precheza.cz> Subject: Re: [R] Superscript and subscrib R for legend x-axis and y-axis and colour different subjects in longitudinal data with different colours Hi Pikal, I looked
2007 Jul 25
2
using contrasts on matrix regressions (using gmodels, perhaps)
Hi, I want to test for a contrast from a regression where I am regressing the columns of a matrix. In short, the following. X <- matrix(rnorm(50),10,5) Y <- matrix(rnorm(50),10,5) lm(Y~X) Call: lm(formula = Y ~ X) Coefficients: [,1] [,2] [,3] [,4] [,5] (Intercept) 0.3350 -0.1989 -0.1932 0.7528 0.0727 X1 0.2007 -0.8505 0.0520
2010 Jul 05
1
Linux-Windows problem
Dear All, I faced the following problem. With the same data.frame the results are different under Linux and Windows. Could you help on this topic? Thanks in advance, Ildiko Linux: > d = read.csv("CRP.csv") > d$drugCode = as.numeric(d$drug) > cor(d, use="pairwise.complete.obs") PATIENT BL.CRP X24HR.CRP X48HR.CRP drug drugCode PATIENT NA
2018 May 21
2
Bootstrap and average median squared error
Dear R-experts, I am trying to bootstrap (and average) the median squared error evaluation metric for a robust regression. I can't get it. What is going wrong ? Here is the reproducible example. ############################# install.packages( "quantreg" ) library(quantreg) crp <-c(12,14,13,24,25,34,45,56,25,34,47,44,35,24,53,44,55,46,36,67) bmi
2018 May 22
2
Bootstrap and average median squared error
I forgot, you should also set.seed() before calling boot() to make the results reproducible. Rui Barradas On 5/22/2018 10:00 AM, Rui Barradas wrote: > Hello, > > If you want to bootstrap a statistic, I suggest you use base package boot. > You would need the data in a data.frame, see how you could do it. > > > library(boot) > > bootMedianSE <- function(data,
2018 May 22
0
Bootstrap and average median squared error
Hello, If you want to bootstrap a statistic, I suggest you use base package boot. You would need the data in a data.frame, see how you could do it. library(boot) bootMedianSE <- function(data, indices){ d <- data[indices, ] fit <- rq(crp ~ bmi + glucose, tau = 0.5, data = d) ypred <- predict(fit) y <- d$crp median(y - ypred)^2 } dat <-
2005 Dec 01
2
Minimizing a Function with three Parameters
Hi, I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and \beta_1, this can be achieved by solving the following three equations: n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) - \sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} 1/(psihat(i)) - (\alpha+1) \sum\limits_{i=1}^{n} ( 1 / (x_i + \psihat(i)) ) = 0 \alpha \sum\limits_{i=1}^{n} (
2017 Jul 19
2
spaghetti plot - urgent
Hi everyone, I?m trying to do a spaghetti plot and I know I?m doing all wrong, It must be. What I need: 15 subjects, each with measurements over 5 different times (t1, ..., t5), and the variable that I need to represent in the spaguetti plot is given by: PCR = b0 + b1 * ti + epsilon B0, - baseline of each subject B1 - trajectory of each subject over time (so multiply by t) Epsilon - error
2018 May 22
1
Bootstrap and average median squared error
Hello, Right! I copied from the OP's question without thinking about it. Corrected would be bootMedianSE <- function(data, indices){ d <- data[indices, ] fit <- rq(crp ~ bmi + glucose, tau = 0.5, data = d) ypred <- predict(fit) y <- d$crp median((y - ypred)^2) } Sorry, rui Barradas On 5/22/2018 11:32 AM, Daniel Nordlund wrote: > On 5/22/2018
2018 May 22
0
Bootstrap and average median squared error
On 5/22/2018 2:32 AM, Rui Barradas wrote: > bootMedianSE <- function(data, indices){ > ???? d <- data[indices, ] > ???? fit <- rq(crp ~ bmi + glucose, tau = 0.5, data = d) > ???? ypred <- predict(fit) > ???? y <- d$crp > ???? median(y - ypred)^2 > } since the OP is looking for the "median squared error", shouldn't the final line of the
2008 Jun 12
1
[7-STABLE] ping -s 4000 with ipsec panic
[FreeBSD 7-STABLE/i386] Hello, I've got a 100 % reproductible panic with ipsec when using a 'ping -s 4000'. It works without ipsec My ipsec setup is very simple, i just use setkey: /etc/ipsec.conf flush; spdflush; add 192.168.1.21 192.168.1.200 esp 1011 -E rijndael-cbc "0123456789012345"; add 192.168.1.200 192.168.1.21 esp 1012 -E rijndael-cbc
2012 Aug 07
6
Decimal number
HI >i have a little problem please help me to solve it >this is the code in R: >> beta0 [1] 64.90614 > beta1 [1] 17.7025 > beta [1] 17 64 >her beta<- c(beta0, beta1) thank you in advance hafida -- View this message in context: http://r.789695.n4.nabble.com/Decimal-number-tp4639428.html Sent from the R help mailing list archive at Nabble.com.
2008 Aug 22
2
WinBUGS with R
Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data <- list("y","X","n","m") #My variables inits.beta <- rep(0,K) inits.beta0 <- 0 inits <-
2011 May 04
1
hurdle, simulated power
Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "power" code is just returning the coefficients, as something isn't working quite right. The
2009 Nov 05
2
faxes received on mISDN
Hi, My initial setup for receiving faxes worked as follows: fax call arrives on ISDN BRI connected to a BOSCH PBX, signal sent to ALCATEL PBX via PRI QSIG then finally sent to ASTERISK via PRI EUROISDN. The Asterisk server then forwarded the call to a iaxmodem and HylaFax received the data. All worked fine. Now I got rid of both BOSCH and ALCATEL in the "fax path" and it's as