Newsgroup members,
I appreciate the help on this topic.
David Duffy provided a solution (below) that was quite helpful, and came
close to what I needed. It did a great job creating two vectors of
dichotomous variables with a known correlation (what I referred to as a
phi-coefficient).
My
You will have to show us the error! It will be shown in
d-p-q-r-tests.Rout.fail (unless this was a segfault or similar).
It is not OK to skip the test, but note that this test is random and does
fail about 1 in 50 times, so you could just try rerunning it.
On Mon, 26 Sep 2005, Jeff Ross
Dear Giacomo,
what you mean precisely saying to create a Trend
variable in a regression? If it concerns about time
series analysis you could give a look to my
contribute:
http://cran.r-project.org/doc/contrib/Ricci-ts-italian.pdf
Best Regards
Vito
Hi,
my name is Giacomo.
I would like to
nikonia == nikonia [EMAIL PROTECTED]
on Mon, 26 Sep 2005 12:15:55 +0200 writes:
nikonia I want to compute the distances in a mixed variable
nikonia matrix using the Gower coefficient. I understand it
nikonia is possible to calculate distances in a matrix with
nikonia mixed
Knut Krueger wrote:
.
I am looking for a histogram or box plot with the adding normal
distribution curve
I think that must be possible, but I am not able to find out how to do.
Regards Knut
__
R-help@stat.math.ethz.ch mailing list
concerning question 1: -- labeling of outliers --
you can get the outliers by boxplot(...)$out
try:
set.seed(17)
x-rexp(99)
names(x)-paste(x,1:99)
out-boxplot(x)$out
text(rep(1.1,length(out)), out, names(out))
Peter Wolf
Stephen D. Weigand wrote:
Dear Yulei,
On Sep 26, 2005, at 6:56 PM,
Le 27.09.2005 10:32, Peter Wolf a écrit :
Knut Krueger wrote:
I am looking for a histogram or box plot with the adding normal
distribution curve
I think that must be possible, but I am not able to find out how to do.
Regards Knut
*There are a lot of answers to add a histogram.
I've got a comlaint about the signature in the last post.
I answered form another computer and there is the automatic signature of
our second business.
On this place is also a part of the equine resarch program in progress.
I did not realize that there was the siganture below.
Sorry for that.
I retract the siggestion I proposed last night -- it was based
on a bad hunch! Sorry for wasting time.
Best wishes,
Ted.
E-Mail: (Ted Harding) [EMAIL PROTECTED]
Fax-to-email: +44 (0)870 094 0861
Date: 27-Sep-05
Dear R users:
I am using the package coda (the last verison in CRAN) to analyse the output
from a MCMC Bayesian analysis. And I get unconsitented results. I have export
the chain using the read.table function and after I have transformed this data
frame to an mcmc object using the mcmc
Hi
***PLEASE do read the posting guide!***
see
?par
?points
?lines
?split.screen
?layout
?matplot
?lattice
HTH
Petr
On 26 Sep 2005 at 18:08, C Tate wrote:
__
R-help@stat.math.ethz.ch mailing list
Dear Juan Pablo,
It is best to send package-specific queries the package maintainer (me
in this case) which you can find out by typeing library(help=coda)
The inconsistency comes from two different ways of estimating the
spectral density at frequency 0: spectrum0() adapted from Heidelberger
This is a good example of where you did not give us enough information to
answer your question. as.Date has a lot of formatting options for working
with dates. If you are still having problems, please provide a short example
(i.e. reproducable code).
Thanks,
Roger
On 9/26/05, Chris Buddenhagen
I'm trying to produce empirical bayes estimates based on the lognormal
model in disease mapping
Is there a way this can be done in R?
thanks
Oarabile
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hello,
I am using the 'regsubsets' function
(from leaps package)
to get the best linear models
to explain 1 variable
from 1 to 5 explanatory variables
(exhaustive search).
Is there anyone who can tell me
on which criterion is based
the 'regsubsets' function ?
Thank you.
samuel
Samuel
Is there anyone who can tell me
on which criterion is based
the 'regsubsets' function ?
the leaps help page says:
cp or adjr2 or r2 is the value of the chosen model selection statistic for
each model
so each of these can be chosen as the selection statistic using the method
argument
hth,
2005/9/26, Dimitris Rizopoulos [EMAIL PROTECTED]:
AFAIK the deriv() is for symbolic derivatives, so I don't know if it
will work in your case.
numericDeriv surely computes numerical gradient. The problem here is
its interface, which requires the vector of symbols involved in the
expression to
I am trying to create some graphs with R and it seems to be able to do
what I need. However, I have so far not been able to find any sort of
explanation of how the graphics system works. I am for instance trying
to create a multiple figure, and I seem to have to call plot.new()
before every new
Hi,
Karin Lagesen wrote:
I am trying to create some graphs with R and it seems to be able to do
what I need. However, I have so far not been able to find any sort of
explanation of how the graphics system works. I am for instance trying
to create a multiple figure, and I seem to have to call
Karin Lagesen [EMAIL PROTECTED] 9/27/2005 9:27:29 AM
I am trying to create some graphs with R and it seems to be able to do
what I need. However, I have so far not been able to find any sort of
explanation of how the graphics system works. I am for instance trying
to create a multiple figure,
Hi.
I have two vectors of gene expression for each of
several days. I want to plot both vectors on the same
plot for a visual representation of up versus down
regulation. I've tried using add=T but that doesn't
work.
eg
plot(Day, gene1)
plot(Day, gene2, add=T)
Any help would be appreciated.
On 1:03:39 am 09/27/05 Prof Brian Ripley [EMAIL PROTECTED] wrote:
You will have to show us the error! It will be shown in
d-p-q-r-tests.Rout.fail (unless this was a segfault or similar).
It is not OK to skip the test, but note that this test is random and
does fail about 1 in 50 times, so
Hi
Just one of those niggles...
I've just been trying to plot a filled circle.
I thought that this would do it
plot(1,1,type=n)
points(1,1,pch=1,bg=blue,cex=5)
#bg: background (fill) color for open plot symbols
But I need to do this instead
points(1.2,1,pch=19,col=blue,cex=5)
Am I
Richard Nixon wrote:
Hi
Just one of those niggles...
I've just been trying to plot a filled circle.
I thought that this would do it
plot(1,1,type=n)
points(1,1,pch=1,bg=blue,cex=5)
#bg: background (fill) color for open plot symbols
But I need to do this instead
Richard Nixon [EMAIL PROTECTED] writes:
Hi
Just one of those niggles...
I've just been trying to plot a filled circle.
I thought that this would do it
plot(1,1,type=n)
points(1,1,pch=1,bg=blue,cex=5)
#bg: background (fill) color for open plot symbols
But I need to do this
Dear R-help,
I was recently wanting to use GEE with the negative binomial family. It
seems that this is lacking in the otherwise excellent implementations of
the GEE methodology ( packages: gee, yags, geepack).
I would have thought it a simple step to allow the creation of a family,
i.e
I am working on automatic optimization of ARIMA parameters.
That takes a lot of computing power, which I would like to reduce by aggregating
and smoothing.
Any thoughts on the subject?
Suggested algorithms?
What is the best order? aggregate then smooth or smooth then aggregate?
Many thanks in
Dear R-users,
i'm a very newbie in linux, but decided to build R from source.
Following the R Installation and Administration manual, i did this:
./configure --enable-R-shlib # this option is here because i intend to
build the GNOME console after...
make
make check
no problems in make check,
This is what happens if you don't have a usable X11 display. Did you
perhaps use a root account on a console owned by a normal user?
On Tue, 27 Sep 2005, Fernando Mayer wrote:
Dear R-users,
i'm a very newbie in linux, but decided to build R from source.
Following the R Installation and
Fernando Mayer [EMAIL PROTECTED] writes:
Dear R-users,
..
Erro: .First.lib failed for 'tcltk'
...
The problem is in the tcltk package. I have searched in the archives,
but didin't find anything similar to this problem (maybe there is, but i
was not able to identify it). I really don't know
On Tue, 27 Sep 2005, Richard Nixon wrote:
Hi
Just one of those niggles...
I've just been trying to plot a filled circle.
I thought that this would do it
plot(1,1,type=n)
points(1,1,pch=1,bg=blue,cex=5)
#bg: background (fill) color for open plot symbols
But I need to do this instead
fem-4.0.1 [and all the other modules, except post] compiles OK, AFAIK,
but when I do a
make check
only 1 test is reported, as follows:
making checksrc
making checkview3d
making checkviewaxis
making checktests
make check-TESTS
$ELMER_HOME undefined, setting it to ../src
CC:
Prof Brian Ripley [EMAIL PROTECTED] writes:
This is what happens if you don't have a usable X11 display. Did you
perhaps use a root account on a console owned by a normal user?
More likely on SUSE is to get done in by a lack of -devel packages.
I've already mentioned tcl-devel and tk-devel
Dear Prof. Ripley, that was exactly what i did. Should i login on the
system as root, instead of login as user and became root via console?
What should i do to have X11 usable?
Thanks,
Fernando Mayer.
Prof Brian Ripley escreveu:
This is what happens if you don't have a usable X11 display.
On Tue, 27 Sep 2005, Fernando Mayer wrote:
Dear Prof. Ripley, that was exactly what i did. Should i login on the system
as root, instead of login as user and became root via console? What should i
do to have X11 usable?
Well, I don't worry about this, but it you want to do a complete test
Hi everyone,
I would like some package of R or any help to solve the next problem
with a weighting fatcors:
Giving a data matrix with dichotomous (2 or more) variables in columns
and individuals in rows, and also a theorical distribution of the
dichotomous variables I would like a row weight
On Tue, 27 Sep 2005, Simon.Bond wrote:
Does anyone have any suggestions for a way round this problem for the
future (I had to resort to using Stata), or maybe more realistically, how
much work it would take to build an extendible version of the gee
algorithm?
I don't think it would take
On 9/27/2005 12:13 PM, Chris Buddenhagen wrote:
Dear Duncan Murdoch
Thanks for responding. I will provide more complete info with future
questions. The data for a small pilot study are below.
The values for each month are the percentage cover of an invasive plant here
in the Galapagos
In data OME in MASS I would like to extract the first 5 observations per
subject (=ID). So I do
library(MASS)
OMEsub - split(OME, OME$ID)
OMEsub - lapply(OMEsub,function(x)x[1:5,])
unsplit(OMEsub, OME$ID)
- which results in
[[1]]
[1] 1 1 1 1 1
[[2]]
[1] 30 30 30 30 30
[[3]]
[1] low low low
Try
plot(Day, gene1)
lines(Day, gene2)
see
?lines
for more details.
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of IAIN
Søren Højsgaard [EMAIL PROTECTED] writes:
In data OME in MASS I would like to extract the first 5 observations per
subject (=ID). So I do
library(MASS)
OMEsub - split(OME, OME$ID)
OMEsub - lapply(OMEsub,function(x)x[1:5,])
unsplit(OMEsub, OME$ID)
- which results in
[[1]]
[1] 1
On Tue, 2005-09-27 at 19:12 +0200, Søren Højsgaard wrote:
In data OME in MASS I would like to extract the first 5 observations
per subject (=ID). So I do
library(MASS)
OMEsub - split(OME, OME$ID)
OMEsub - lapply(OMEsub,function(x)x[1:5,])
unsplit(OMEsub, OME$ID)
- which results in
On Tue, 27 Sep 2005, Ferran Carrascosa wrote:
Hi everyone,
I would like some package of R or any help to solve the next problem
with a weighting fatcors:
Giving a data matrix with dichotomous (2 or more) variables in columns
and individuals in rows, and also a theorical distribution of the
Hi,
I have the following codes to replace missing using median, assuming missing
only occurs on continuous variables:
trn1-read.table('trn1.fv', header=F, na.string='.', sep='|')
# median
m.trn1-sapply(1:ncol(trn1), function(i) median(trn1[,i], na.rm=T))
#replace
trn2-trn1
for (each in
Can someone please show me what I need to get something like this to
work
for(a in c(1:5)){
data$a-c(a:10)
}
so that i end up with a structure
data$1-[1,2,3,4,5,6,7,8,9,10]
data$2-[2,3,4,5,67,8,9,10]
data$3-[3,4,5,67,8,9,10]
data$4-[4,5,67,8,9,10]
data$5-[5,67,8,9,10]
thanks loads
Tom
I include my workaround for object.size failing on RODBC objects.
Code adapted from Petr Pikal's code
(http://tolstoy.newcastle.edu.au/R/help/04/06/1228.html).
ls.objects-function (pos = 1, order.by,...)
{
napply - function(names, fn) sapply(names, function(x) fn(get(x,pos=pos)))
names -
Le 27 Septembre 2005 10:39, tom wright a écrit :
Can someone please show me what I need to get something like this to
work
for(a in c(1:5)){
data$a-c(a:10)
}
so that i end up with a structure
data$1-[1,2,3,4,5,6,7,8,9,10]
data$2-[2,3,4,5,67,8,9,10]
data$3-[3,4,5,67,8,9,10]
Use 'which(...arr.ind=T)'
x.1
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 6 10 3 4 10 7 9 8 4 10
[2,] 8 7 4 7 4 8 3 NA 3 4
[3,] 7 7 10 10 3 5 3 2 2 2
[4,] 3 4 5 10 10 2 6 9 4 5
[5,] 3 5 9 5 6 NA 3 NA 6 7
[6,] 9 6 10 5 10 4 2 10 NA 5
[7,] 5 2 5 10 3 7 6 4 6 8
[8,] 2 6 1 8 9 2 7 8 3 8
When sapply is used on a character vector it will use those
as the names so:
data - sapply(as.character(1:5), function(x) seq(as.numeric(x),10))
will give data with the required names.
On 9/27/05, tom wright [EMAIL PROTECTED] wrote:
Can someone please show me what I need to get something like
Dear all,
I'm trying to combine in R the two following datasets.
I'm would be very grateful if someone could help me with this problem.
Dataset 1:
X Y Z
1 27 48
2 25 50
3 27 40
4 28 56
Where X is a unique number and Y and Z some variables
Hi Jerome,
look at the merge() function
Nolwenn
On Tue, 27 Sep 2005, Jérôme Lemaître wrote:
Dear all,
I'm trying to combine in R the two following datasets.
I'm would be very grateful if someone could help me with this problem.
Dataset 1:
X Y Z
1 27 48
2 25
Hmm my bad,
Thanks for your replies but I think my example was a little to simple
the actual code I'm using is:
f_haardisolve-function(v_dataset){
#pad data to make length a power of 2
v_dataset-f_paddata(v_dataset)
l_cooef-list() #holder for cooefficents
i_count-1
On Mon, Sep 26, 2005 at 07:05:33PM +0200, Göran Broström wrote:
On Mon, Sep 26, 2005 at 09:27:56AM +0200, Martin Maechler wrote:
I'm crossposting to the ESS-help mailing list which is slightly
more appropriate here. [This may be a rare case where
crossposting seems to make much sense.]
Got it thanks
the trick was to create a variable with the new name
str_idx-paste('c',i,sep='')
then use double square brackets
l_cooef[str_idx]]-whatever()
f_haardisolve-function(v_dataset){
#pad data to make length a power of 2
v_dataset-f_paddata(v_dataset)
l_cooef-list() #holder
Hello ,
How about cbind or rbind?
You could create a new data frame and push yours into it.
Greetz
Jérôme Lemaître schrieb:
Dear all,
I'm trying to combine in R the two following datasets.
I'm would be very grateful if someone could help me with this problem.
Dataset 1:
X Y Z
1
Your question seems almost too broad, but I will attempt a brief
response.
Have you reviewed the material on matrices in An Introduction to R,
avaliable via, e.g., help.start()? Also, have you reviewed LAPACK
(http://www.netlib.org/lapack/)?
If the answer to
Hi,
I'm trying to understand environment object in R.
I used the example:
f - function(x) {
y - 10
g - function(x) x + y
return(g)
}
h - f()
h(3)
then i saw that f return an environment
h
function(x) x + y
environment: 01B28570
but I coudn't access to x
Try 'get(x, env = environment(h))'
-roger
Ron Ophir wrote:
Hi,
I'm trying to understand environment object in R.
I used the example:
f - function(x) {
y - 10
g - function(x) x + y
return(g)
}
h - f()
h(3)
then i saw that f return an
Ron Ophir [EMAIL PROTECTED] writes:
Hi,
I'm trying to understand environment object in R.
I used the example:
f - function(x) {
y - 10
g - function(x) x + y
return(g)
}
h - f()
h(3)
then i saw that f return an environment
h
function(x) x +
Hej,
but in your function you add x and y to 10 and 3, so your values are
merged to one value available in g variable.
And now you want to see what was your y and your x?
I guess I do not get the idea of your question.
Well, then you could return y and x as well as g.
Greetz n god luck.
Roger
Uau!
This package is very interesting, very thank you Thomas.
Thanks to all!
2005/9/27, Thomas Lumley [EMAIL PROTECTED]:
On Tue, 27 Sep 2005, Ferran Carrascosa wrote:
Hi everyone,
I would like some package of R or any help to solve the next problem
with a weighting fatcors:
Giving
I am afraid this does look like a real problem, if a minor one. We have
for the first problem
x - 10^(ex - c(1,2,5*(1:5),50,100,200,300,Inf))
ex - -c(rev(1/x), ex)
qcauchy(ex, log=TRUE)
The first entry of the result should be Inf and the last -Inf. From the
output you have shown us I would
On Tue, 27 Sep 2005, Bliese, Paul D LTC USAMH wrote:
My situation is a bit more complicated and I'm not sure it is easily
solved. The problem is that I must assume one of the vectors is
constant and generate one or more vectors that covary with the constant
vector.
One way is to sample from
On Sun, Sep 25, 2005 at 12:07:02PM +0100, Prof Brian Ripley wrote:
On Sun, 25 Sep 2005, Andrew Robinson wrote:
Hi R-helpers,
I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R
Version 2.1.1.
I have constructed a package using package.skeleton(), when I try
$ R
Hi
I have a largish dataset (26 columns 35000 rows) which I have been
subjecting to logistic regression and support vector machine analysis.
I have noticed that R easily copes with using the data in either
technique. Now I have to try and see what the best modeling technique
to use is.
I
Shengzhe Wu [EMAIL PROTECTED] wrote:
In the Unix environment, I open a window by x11(). May I specify the
position of this window by specifying the position of the top left of
the window as in Windows environment?
I use xwit (version 3.4), a system command which manipulates existing
X
On Tue, 27 Sep 2005, Prof Brian Ripley wrote:
I am afraid this does look like a real problem, if a minor one. We have for
the first problem
x - 10^(ex - c(1,2,5*(1:5),50,100,200,300,Inf))
ex - -c(rev(1/x), ex)
qcauchy(ex, log=TRUE)
The first entry of the result should be Inf and the
Hello,
I would like to plot my data in a fashion similar to a boxplot, but
plot the true data points without a box, just overlay lines of the
summary generated by the boxplot. I have less than 10 observations
within each group, so I think showing the actual data would be more
effective than
Try this:
boxplot(Sepal.Length ~ Species, iris)
with(iris, stripchart(Sepal.Length ~ Species, vertical = TRUE, add = TRUE))
On 9/27/05, Aric Gregson [EMAIL PROTECTED] wrote:
Hello,
I would like to plot my data in a fashion similar to a boxplot, but
plot the true data points without a box,
Hi R Users,
How to use the function polygon()
together with the package scatterplot3d?
I am trying to color below of the curves
defined for the function density().
I tried to use the site: R GRAPH GALLERY
as tutorial.
I tried to adapt the example of this page:
[figure]:
70 matches
Mail list logo