Yes, Doug is correct and I'm wrong. In fact, his comment jogged MY memory --
I actually used BMDP a bit in the late 70's(I think it was).
Thanks to Doug for corrected chronology.
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [ma
chipmaney schreef:
I have run a kruskal.test() using the by() function, which returns a list of
results like the following (subset of results):
Herb.df$ID: 4-2
Kruskal-Wallis chi-squared = 18.93, df = 7, p-value = 0.00841
---
Shawn Morrison-2 wrote:
>
> # The paper reports a 95% CI of 0.79 - 1.10
> # "My" reproduced result for the CIs is much larger, especially on the
> upper end. Why would this be?
> # The authors report using the 'delta' method (Caswell, 2001) to
> calculate the CI in which the
>
>
Shawn,
I p
Try this:
sapply(x, '[', 'p.value')
On Thu, Feb 18, 2010 at 5:21 PM, chipmaney wrote:
>
> I have run a kruskal.test() using the by() function, which returns a list of
> results like the following (subset of results):
>
> Herb.df$ID: 4-2
> Kruskal-Wallis chi-squared = 18.93, df = 7, p-value
Thanks, as a follow-up, how do i extract the list element name (ie, 4-2 or 44-1)
thanks,
chipper
Date: Thu, 18 Feb 2010 11:56:45 -0800
From: ml-node+1560750-540257540-69...@n4.nabble.com
To: chipma...@hotmail.com
Subject: Re: Extracting values from a list
Try this:
sapply(x, '[', 'p.value
hi, I am new to Linux and R environment. I have a existing R script. I
wonder how to open my R script on Linux platform and execute selected
written R command?
thanks
Xin
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
(Corrections/additional information welcome!)
My recollection is that the BMD programs (which, in a later version,
became BMDP) predated SAS and were specifically for BioMeDical
analysis.
How could I forget those! Yes, my old (as in 1980-1985) boss at the
University hospital even had the
I have data (each Y_i is a vector) in the form of
Y_i = X_i \beta_i + Z_i b_i + epsilon_i
Were it not for the measurement error (the epsilon_i) it's a very
simple model --- nice and balanced, compound symmetry, and I'd just
use lme(y ~ x1 + x2, random=~1|subj, ...) but the measurement
At the very simplest level, you can open the R script in a text editor.
Then start R at the command line in the terminal by typing
R
and then copy and paste the desired lines from the text file into the
terminal.
There are plenty of more elegant and complex solutions, but that
one is easy and will
Hello,
if You run R in terminal, You can type R to start the R software, then use
command source("/path/to/your/script") to load the script and then You call
functions from the script as usual. Optionally, You can move to the directory
with script with command setwd("/path/to/the/directory"). Or
Hello All,
I'm having some issues controlling graphics in R. I was wondering if anyone
may help me tackle this problem:
Given a data frame "X" with variables "Year", "Factor" (w/ n groups), and
"Freq"
How do I create a single graphic with the following plots aligned in a
vertical stack?
1. box
Hi,
I have a nonlinear constrained optimization code that can handle nonlinear
(and linear) constraints. Send me an email if you are interested.
There is also a package available on R-forge, called Rsolnp.
Hope this is helpful,
Ravi.
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Peter Dalgaard
> Sent: Thursday, February 18, 2010 12:44 PM
> To: Douglas Bates
> Cc: r-help@r-project.org; Bert Gunter
> Subject: Re: [R] Use of R in clinical trials
>
>
> >> (C
Use outer margins. Try something like:
> par(mfrow=c(3,1), mar=c(0,4,0,2)+0.1, oma=c(5,0,3,0)+0.1 )
Then do your plots without resetting margins.
Also you can use xaxt='n' rather than axes=FALSE to suppress just the x axis
and not have to do the y axis and box by hand.
Hope this helps,
--
G
Hi Folks,
I want to apply cluster analysis on a categorical data set, could you
recommend me some R package and suggestion?
Thanks!
Dong
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
Hi everyone!
Is it possible to save an image of the workspace where
1) Packages
2) Classes
are saved along with the image?
Until now I only managed to save an workspace image that contained all
variables (including functions). When loading this image back into a new
session, p
Thanks for helping!
--
View this message in context:
http://n4.nabble.com/plot-more-series-or-more-colors-tp1556751p1560863.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
Without know what your data set really looks like, I'd look to decision
trees - specifically package rpart and use method = classify.
Your problem may not be appropriate in that environment, but it is hard to
say with limited explanation of issues.
good luck
Steve Friedman Ph. D.
Spatial Statist
Hi,
I'm generating a mixture of 2 univariate normal distributions using norMix and
rnorMix and would like to put a constraint on Mean (Equal). here is my code
snippet:
library(nor1mix)
X <- norMix(mu=c(50, 60 ), sig2=c(5,10), w=c(0.5, 0.5))
mixData <- rnorMix(1000, X)
plot(mixData, type="l"
On Feb 18, 2010, at 3:50 PM, Nordlund, Dan (DSHS/RDA) wrote:
>> -Original Message-
>> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
>> Behalf Of Peter Dalgaard
>> Sent: Thursday, February 18, 2010 12:44 PM
>> To: Douglas Bates
>> Cc: r-help@r-project.org; Bert
Anyone have any recollection of Prophet software, from the National (US)
Center for Research Resources?
--Chris
Christopher W. Ryan, MD
SUNY Upstate Medical University Clinical Campus at Binghamton
425 Robinson Street, Binghamton, NY 13904
cryanatbinghamtondotedu
"If you want to build a ship,
On Feb 18, 2010, at 4:26 PM, Christopher W. Ryan wrote:
> Anyone have any recollection of Prophet software, from the National (US)
> Center for Research Resources?
>
> --Chris
A quick Google search comes up with a very dated site at Northwestern:
http://www.basic.northwestern.edu/biotools/pr
What steps can be take to append data to a list?
Suppose I have the following list and want to append a z axist to the list?
pts <- list(x=cars[,1], y=cars[,2])
z<-rnorm(max(dim(cars)))
How would I go about appending z to an existing list?
Thanks a ton...
___
Hi r-users,
I have 2 sets of data and I would like to superimpose this cumulative density
in one graph.
I know how to put the 2 graphs in one same graph but my problem is the data are
different.
> z[1:20]
[1] 2.02347388 3.19514379 0.05188875 1.41333812 3.50249892 4.34272676
6.6563
On Feb 18, 2010, at 6:15 PM, Jason Rupert wrote:
What steps can be take to append data to a list?
Suppose I have the following list and want to append a z axist to
the list?
pts <- list(x=cars[,1], y=cars[,2])
z<-rnorm(max(dim(cars)))
How would I go about appending z to an existing list?
On Thu, 18 Feb 2010, Richard Valliant wrote:
Should the svyby function be able to work with svyquantile? I get the
error below ...
It works, but you need to either specify ci=TRUE or keep.var=FALSE. The
problem is that svyquantile() by default does not produce standard errors.
svyby(~api0
?c
but you have to make sure z is a list:
c(pts,z) ## probably is not what you want
c(pts,list(z)) ## probably is, but z will be unnamed
c(pts,list(z=z) ## names z "z"
Cheers,
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org
Hi
Thank you very much for your reply.
I meant I got different RMSE with different runs to the same program without
any change for the dataset or the parameters of SVM.
This is my code:
library(e1071)
readingmydata <- as.matrix(read.delim("mydataset.txt"))
train.x <- readingmydata[,-1]
trai
## The grouped boxplot is one of the features included in the HH package.
## You will need to install HH if you do not yet have the HH package
## A similar example is posted on my website
##http://astro.ocis.temple.edu/~rmh/HH/bwplot-color.pdf
## This is fake data which I hope mimics the struc
Greg Snow wrote:
Use outer margins. Try something like:
par(mfrow=c(3,1), mar=c(0,4,0,2)+0.1, oma=c(5,0,3,0)+0.1 )
Then do your plots without resetting margins.
Thanks. Perfect! This little detail has been bothering me for quite a
while...
Also you can use xaxt='n' rather than axes=FALS
For the following:
Bldgid<-c(1000,1000,1000,1001,1002,1003,1003,1003)
Maplot<-c(2,20001,20002,3,30001,4,40001,40002)
Area<-c(40,170,160,50,100,100,90,110)
#Construct Sample dataframe
MultiLotBldgs..<-data.frame(Bldgid,Maplot,Area)
CondoLots_ <- tapply(MultiLotBldgs..$Maplot, Mul
I think you can just use
CondoLots[CondoLots > 1] <- CondoLots[CondoLots > 1] -1
-Ista
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and pr
Thanks for all the replies!
I think I posted a poor example, but now I understand that you can name entries
in a list. That is really neat.
Before I was just letting them default, which is kind of problematic.
I think I also came across another way and also used the suggestion from Bert
Ok, took me a while, but I figured it out. Because my running mean had less
years than my standard rainfall graph, when I overlaid the running mean onto
the rainfall it was trying to stretch out. So I just plotted both onto the
same graph., like so:
barplot(Ann,main=title, xlab="Year",ylab="Rainf
On Tue, 16 Feb 2010, Esmail wrote:
On 16-Feb-10 09:03, Karl Ove Hufthammer wrote:
On Tue, 16 Feb 2010 08:00:09 -0500 Esmail wrote:
And along the same lines, any type of interactive debugging
utility for R?
See this article in R News:
'Debugging Without (Too Many) Tears'
http://cran.r-proje
Thanks guys. I ended up doing as you suggested Dieter. Thanks for the idea :)
--
View this message in context:
http://n4.nabble.com/sql-query-variable-tp1558189p1561158.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.
I am old enough to have lived through this particular transition.
Prior to the advent of SAS, trials were analyzed by in-house written
programs (usually in Fortran maybe with the help of IMSL). These
programs were huge card decks. Having the card reader eat a card
half way through reading the deck
Hello there,
I have a bunch of histogram bars that I'd like the first to be a certain
colour, second to be another colour, third to be a third colour, and repeat
for all my 39 bars.
I thought this was the way to go, but I get the same cyan coloured bars for
all the bars. I did a vector of 3 colo
Dear List,
I would like to ask about package RGtk2 with which I have a problem.
I will very much appreciate if somebody could tell me what I need to do.
I need to install a package scdMicro and it depends on gWidgetsRGtk2.
I am working on MAC, version 10.5.8.
When I try to load gWidgetsRGtk2 (or R
Have you tried install GTK2 manually first?
It can be downloaded from http://r.research.att.com/gtk2-runtime.dmg
HTH
Ronggui
On 19 February 2010 10:02, Anna Oganyan wrote:
> Dear List,
> I would like to ask about package RGtk2 with which I have a problem.
> I will very much appreciate if somebo
Hi,
I would like to distinguish my plots using colors but I got error message. How
do I correct that?
plot(ecdf(z), main ="CDF for observed and simulated weighted
sum",type="l",lwd=2,col="blue",
xlab="Weighted sum (mm)", ylab="Cumulative Percent", xlim=c(0,15), xaxs ='i',
yaxs ='i',ylim=c
I take back what I said about not being able to use RGtk2 with the 64-
bit mac. I (just today, for the first time ) can now get RGtk2 to load
and then to get rgl to run properly. I had been able to run rgl in the
32 bit system but that was not satisfactory. I had earlier installed
the 2.18
Dear R helpers,
I have two queries.
(1) If the dataset contains some variables having all the entries = 0 and while
analysing I want to delete those pericular columns, how do acheive this. i.e.
dataset1
sr_no var1 var2 var3 var4 var5
1 5 0
On Feb 19, 2010, at 12:38 AM, Roslina Zakaria wrote:
Hi,
I would like to distinguish my plots using colors but I got error
message. How do I correct that?
plot(ecdf(z), main ="CDF for observed and simulated weighted
sum",type="l",lwd=2,col="blue",
xlab="Weighted sum (mm)", ylab="Cumula
Dear Anna,
19.02.2010 08:17, Anna Carter wrote:
> (1) If the dataset contains some variables having all the entries = 0
> and while analysing I want to delete those pericular columns, how do
> acheive this. i.e.
Let's suppose 'df' is your data frame, then:
subset(df, select=which(colSums(df)!=0)
On Feb 19, 2010, at 1:36 AM, K. Elo wrote:
Dear Anna,
19.02.2010 08:17, Anna Carter wrote:
(1) If the dataset contains some variables having all the entries = 0
and while analysing I want to delete those pericular columns, how do
acheive this. i.e.
Let's suppose 'df' is your data frame, the
On Feb 18, 2010, at 11:32 PM, Wincent wrote:
Have you tried install GTK2 manually first?
It can be downloaded from http://r.research.att.com/gtk2-runtime.dmg
That is for Tiger builds.
HTH
Ronggui
On 19 February 2010 10:02, Anna Oganyan
wrote:
Dear List,
I would like to ask about pac
> (2) Suppose I have variable no of datasets 'say n = 10'. I wish to write a
> loop assigning each of these datasets to diffrent csv files e.g.
>
> for (i in 1:10)
> {
> write.csv(data.frame(dataset[,,i]), 'data_set[i].csv', row.names = FALSE)
> }
>
> The result of this command is generation of
(Forgot to cc. reply to K. Elo, apologies if you get it twice)
K. Elo wrote:
Dear Anna,
19.02.2010 08:17, Anna Carter wrote:
(1) If the dataset contains some variables having all the entries = 0
and while analysing I want to delete those pericular columns, how do
acheive this. i.e.
Let's sup
101 - 149 of 149 matches
Mail list logo
