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