Dennis Alexis Valin Dittrich [EMAIL PROTECTED]writes:
ps2pdf14 -dPDFSETTINGS=/prepress figure.eps tmp.pdf
pdftops -eps tmp.pdf figure_with_embeded_fonts.eps
The first script is delivered with ghostscript. The second program
belongs to the xpdf package.
Thanks a lot.
That was a miracle. The
On Wed, Oct 11, 2006 at 03:24:13PM +0900, Yury Yuryev wrote:
Dennis Alexis Valin Dittrich [EMAIL PROTECTED]writes:
ps2pdf14 -dPDFSETTINGS=/prepress figure.eps tmp.pdf
pdftops -eps tmp.pdf figure_with_embeded_fonts.eps
The first script is delivered with ghostscript. The second program
victor wrote:
Uwe Ligges wrote:
Dominique Katshunga wrote:
Can someone help interprete the error message below? i was trying to
load the package copula from the R command prompt.
Error in loadNamespace(package, c(which.lib.loc, lib.loc), keep.source =
keep.source) :
in 'copula' methods
Dear all,
I want to see if the treatment of an animal with a specific compound has an
effect on the expression of certain genes. Though my question is based in
biology, it really is all about how to deal with the standard deviation in
normalised data.
I have three groups of animals; untreated,
The arch test requires a varest object and I am trying to
write one that
will use R-metrics arima, arch, Garch objects, or at least a vector.
Also the arch function has the following line of code that I can not
find the function anywhere
archs.resids - apply(resids, 2, function(x)
Hi
You definitely shall consult an introduction manuals you can find in
doc directory of your R installation. Good starting point is also
http://www.r-project.org/ and its documentatation section.
boxcox - function(x,min,max,step) {
lambda - seq(min,max,step) # you set lambda1
s -
I'd suggest looking at the boxcox function in the MASS package, which might
show you how to change your code.
It would be easier to diagnose the problem if you included a small example
showing the errors. Some possibilities are that if x is a vector, nrow(x)
will return NULL, that you are using
Hi!
I'd like to plot things with axes going from the
highest to the lowest value, so that e.g. high values
on the xaxis are plotted to the left and low values to
the right.
Could anyone tell me how this is done? I couldn't find
anything in the documentation.
Thanks a lot!
Silli
Silvia Lipski [EMAIL PROTECTED] writes:
Hi!
I'd like to plot things with axes going from the
highest to the lowest value, so that e.g. high values
on the xaxis are plotted to the left and low values to
the right.
Could anyone tell me how this is done? I couldn't find
anything in the
On Wed, 2006-10-11 at 02:27 -0700, Silvia Lipski wrote:
Hi!
I'd like to plot things with axes going from the
highest to the lowest value, so that e.g. high values
on the xaxis are plotted to the left and low values to
the right.
Could anyone tell me how this is done? I couldn't find
I'd like to plot things with axes going from the
highest to the lowest value
You can use the xlim and ylim parameters of the plot function. E.g:
x = 1:100
plot(x, x^2) # normal
plot(x, x^2, xlim=c(100,0)) # reversed x axis
cu
Philipp
--
Dr. Philipp Pagel
Hello,
I'd like to overlay two histograms for comparison. I know that this
can be done, in fact I have done it myself several times. But right
now I can't neither remember nor find any documentation.
Is there anywhere a comprehensible guide for making graphics / plots
with R? I know there are
Hi
something like that?
y-1:10
x-11:20
plot(x,y)
plot(x[10:1],y) # reverse x
HTH
Petr
On 11 Oct 2006 at 2:27, Silvia Lipski wrote:
Date sent: Wed, 11 Oct 2006 02:27:34 -0700 (PDT)
From: Silvia Lipski [EMAIL PROTECTED]
To:
Joe Byers schrieb:
JGR, when I manually compiled it for my system by specifying my location
of java 1.5.0_06, worked fairly nicely. The autoinstall does not find
my installation of this java version since Redhat EL4 does not support
it. Redhat EL5 will in the near future. Maybe JGR will
On Wed, 11 Oct 2006, Stefan Grosse wrote:
Joe Byers schrieb:
JGR, when I manually compiled it for my system by specifying my location
of java 1.5.0_06, worked fairly nicely. The autoinstall does not find
my installation of this java version since Redhat EL4 does not support
it. Redhat EL5
Dear all,
I would like to estimate the parameters of a mixture of two beta
binomial distributions. I saw the function fitdistr from the MASS
library could be a solution but I have some difficulties. Does anyone
know how to specify to fitdistr another function as these which are
predefined? Or is
Hello all,
I am trying to explore random forest in R. What I want to do is get the node
number in which the case falls in the tree of random forest. For that I am
calling the predict method as:
learn.pred - predict (learn.rf, newdata=learn.data.x,norm.votes=
TRUE,predict.all = TRUE, nodes=
On Wed, 11 Oct 2006, Caroline TRUNTZER wrote:
Dear all,
I would like to estimate the parameters of a mixture of two beta
binomial distributions. I saw the function fitdistr from the MASS
library could be a solution but I have some difficulties. Does anyone
know how to specify to fitdistr
Do provide a reproducible example, as the Posting Guide suggests.
Try:
library(randomForest)
example(predict.randomForest)
iris.pred - predict(iris.rf, iris[ind == 2,], nodes=TRUE)
str(iris.pred)
attr(iris.pred, nodes)
Andy
From: Rupendra
Hello all,
I am trying to explore random forest
Dear All,
I am willing to test homogenity of variances using Levene's test. I know
how to do it with simple model, but I am lost when dealing with a
little bit more complex model.
I have response variable (y) and four factors (a, b, c, d). One of
these four factors (d) is nested within another
January == January Weiner [EMAIL PROTECTED] writes:
Dear all, I am having troubles importing values written as
scientific notation using read.table(). I'm sure this is a
frequent problem, as many people in my lab have this
problem as well, so I'm sure that I just have
I am a teaching asistant at Ankara University. I try to analize ABD stock
market returns.I want to compare coefficients of variation. But some of means
are positive, some of means are negative. In this instance , how can compare
coefficients of variation?
Could you help me,please?
Thanks...
Hello all,
I am trying to explore random forest in R. What I want to do is get the node
number in which the case falls in the tree of random forest. For that I am
calling the predict method as:
learn.pred - predict (learn.rf, newdata=learn.data.x,norm.votes=
TRUE,predict.all = TRUE, nodes=
Prof Brian Ripley schrieb:
Which arch and how did you install R? For me, it works with R
compiled from the sources (but not installed from RPM) on FC5 i686,
but not x86_64.
The issue is a run-time one:
I have i686, no x86_64 ...
OK, so I have to compile R from source ...
Thanks for the
Hi,
is there anyone who successfully used JGR (Java GUI for R) on Linux
machines, more specific on 64bit AMD? Tried to get JGR running but after
installing it as recommended (all installs ok), got a segfault with a
very nice implement me message:
[EMAIL PROTECTED]: R
Hi Petr,
Thanks for your response. I have data that looks like the following:
sample 1 sample 2 sample 3
red candy400 300 2500
green candy1000 200
black candy 300
On 11 October 2006 at 16:05, Stefan Grosse wrote:
| Prof Brian Ripley schrieb:
|
| Which arch and how did you install R? For me, it works with R
| compiled from the sources (but not installed from RPM) on FC5 i686,
| but not x86_64.
| The issue is a run-time one:
|
| I have i686, no x86_64
Note: this is advocacy for education in clear quantitative
language and is a border-line off topic rant...
The other day I read a paper from a student who used notation
like 2e-4 in the text - blech! I sent it back for revisions.
You have sent it back for revisions just because the student
Hi, my post about failing JGR on amd64 solved and here is how I got it
running, hopefully it helps others.
Suppose your JDK (and you need JDK, not only JRE) is in
/usr/lib/jvm/java-1.5.0-sun (thus JRE in
/usr/lib/jvm/java-1.5.0-sun/jre). Set JAVA_HOME to point there to JDK.
Run R CMD
Hi, my post about failing JGR on amd64 solved and here is how I got it
running, hopefully it helps as well.
Suppose your JDK (and you need JDK, not only JRE) is in
/usr/lib/jvm/java-1.5.0-sun (thus JRE in
/usr/lib/jvm/java-1.5.0-sun/jre). Set JAVA_HOME to point there to JDK.
Run R CMD
Hi
a litle bit different story. But
x1 - sample(c(rep(red,400),rep(green, 100),
rep(black,300)),100)
is maybe close. With data frame (if it is not big)
DF
color sample1 sample2 sample3
1 red 400 3002500
2 green 100 0 200
3 black 3001000 500
x -
Hi
I don't know if this quaifies as an R question, but it relates to using R.
I use R under Linux and use Kate as my script editor. I would like to be
able to pipe the selected block to an R instance which is in an xterm
session (not the console in Kate!).
Is it possible to configure Kate
Here's a way using apply(), and the prob= argument of sample():
df - data.frame(sample1=c(red=400,green=100,black=300),
sample2=c(300,0,1000), sample3=c(2500,200,500))
df
sample1 sample2 sample3
red 400 3002500
green 100 0 200
black 3001000 500
Hi,
A = rep(rep(c(0,1),4),2)
B = rep(rep(c(0,1),each=4),2)
C = rep(rep(rep(c(0,1),2),each=2),2)
X = data.frame(A,B,C,rnorm(16))
X = data.frame(A,B,C,D=rnorm(16))
X
A B C D
1 0 0 0 0.829987
2 1 0 0 1.319203
3 0 0 1 0.223752
4 1 0 1 0.017353
5 0 1 0 1.288284
6 1 1 0
Hello,
I'm reposting from last week as I had no response. I'm considering to
post a bug report instead, but out of respect for the developers I'd
apreciate if someone could confirm the problem before I commit a bug
report.
Thanks.
Forwarded Message
Subject: [R] RODBC: longest
Hi,
First of all, thanks for the great work on R in general, and MASS in
particular. It's been a life saver for me many times.
However, I think I've discovered a bug. It seems that, when I use
weights during an initial least-squares regression fit, and later try to
add terms using
I have had a problem in finding the minimum of a function, the function in
cause is:
curv - function(a,b){
date - read.table(bessa.csv,header=T,sep=;,dec=,)
calP - (22000)/(1+exp(-(a*date$v+b)))
err - (calP-date$P)^2
return(sum(err))
}
It would like to know which function that I must use to
Hello,
I would like to know whether is possible to install R in a Itanium 64
bits computer. We cannot neither compile source code because we have
Intel compilers.
Specifically, we have an Intel Itanium II, Linux Red Hat 64 bits and
Intel Fortran and C compilers: ifort, icc.
Thanks
--
jm~
Greetings:
I've searched the R archives with no luck.
I want to print this to the screen as part of on-screen instructions as an
example:
default.FACTOR.labels - c(Probe1, Probe2, Probe3)
I can't seem to trick gsub()
gsub(', \, default.FACTOR.labels - c('Probe1', 'Probe2', 'Probe3')))
[1]
On Wed, 2006-10-11 at 13:30 -0400, Charles Annis, P.E. wrote:
Greetings:
I've searched the R archives with no luck.
I want to print this to the screen as part of on-screen instructions as an
example:
default.FACTOR.labels - c(Probe1, Probe2, Probe3)
I can't seem to trick gsub()
D'oh!
I've been using cat() but somehow never got the bigger picture.
Thanks!!!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: Marc Schwartz [mailto:[EMAIL PROTECTED]
Sent: Wednesday,
Hello all,
I'll bet this is not complicated, but I really did hunt for information
on it, and tried some ideas without success.
Running Windows XP and R 2.3.1
I am using the Sound package to read in .wav files, and this works fine.
The R object produced is of class sample. I want to use some
On 10/11/06, Marc Schwartz [EMAIL PROTECTED] wrote:
On Wed, 2006-10-11 at 13:30 -0400, Charles Annis, P.E. wrote:
Greetings:
I've searched the R archives with no luck.
I want to print this to the screen as part of on-screen instructions as an
example:
default.FACTOR.labels -
Ricardo Bessa wrote:
I have had a problem in finding the minimum of a function, the
function in cause is:
The first, obvious, optimization step is that you
are reading the whole table every time you call the function.
You should read the table just once. Instead of:
curv - function(a,b){
I'm not familiar with seewave package, but I imagine what you have to do is
extract the sound component of the sample object and convert that to a data
frame.
waveLeft - 2*((seq(0,80,length=88200)%%1^2)-.5)
s - as.Sample(waveLeft,44100,16)
ds - data.frame(s$sound)
is.data.frame(ds)
[1] TRUE
I tried all of the approaches below.
the problem with:
x - data.frame(matrix(NA,100,3))
for (i in 2:ncol(DF)) x[,i-1] - sample(rep(DF[,1], DF[,i]),100)
if you want result in data frame
or
x-vector(list, 3)
for (i in 2:ncol(DF)) x[[,i-1]] - sample(rep(DF[,1], DF[,i]),100)
is that this code
Thanks, Marc,
I tried it and it didn't really work.
x.num-as.numeric(x)
is.numeric(x.num)
[1] TRUE
x.num-readClipboard()
is.numeric(x.num)
[1] FALSE
is.character(x.num)
[1] TRUE
I use readClipboard to take in data, and it seems after the data is taken
in, x.num changed to character.
I just
Yulei,
It would appear that the default mechanism for the function (which
appears to be Windows specific) is to read in the data from the current
system clipboard as a _character vector_. Thus, perhaps something like:
x.num - as.numeric(readClipboard())
would be helpful.
With respect to
Two things to note:
(1) rep() can be vectorized:
rep(1:3, 2:4)
[1] 1 1 2 2 2 3 3 3 3
(2) you will likely get much better performance if you work with
integers and convert to strings after sampling (or use factors), e.g.:
c(red,green,blue)[sample(rep(1:3,c(400,100,300)), 5)]
[1] red blue
Marc, Thanks! x.num - as.numeric(readClipboard()) works perfect.
Meanwhile, I have tried to use odbcConnectExcel from RODBC pack. It seems
to work as well. I am able to read in the data, but I need to construct a
sql table from the data, and use sqlFetch to retrieve and assign them to
Hi,
I have questions about object attributes, and how they are handled when
subsetted. My examples will use:
tm - (1:10)/10
ds - (1:10)^2
attr(tm,units) - sec
attr(ds,units) - cm
dat - data.frame(tm=tm,ds=ds)
attr(dat,id) - test1
When a primitive class object (numeric, character, etc.) is
Hello,
I get this following error when try the configure command. Any suggestions
would be helpful
configure: error: --with-readline=yes (default) and headers/libs are not
available
Not sure whats happening here. I am trying to install on a Redhat Linux ES v4.0
on a AMD 64 opteron machine.
[EMAIL PROTECTED] writes:
Hello,
I get this following error when try the configure command. Any suggestions
would be helpful
configure: error: --with-readline=yes (default) and headers/libs are not
available
Not sure whats happening here.
Well, readline headers and libraries are
Yes, we have an SGI Altix that is similar and R installs just fine. You
should be able to compile
the source code with either the gnu or Intel compilers.
See the R Installation and Administration manual at
http://www.cran.r-project.org/doc/manuals/R-admin.pdf,
and pay particular attention to
You can define your own class then define [ to act any way you would like:
[.myobj - function(x, ...) {
y - unclass(x)[...]
attributes(y) - attributes(x)
y
}
tm - structure(1:10, units = sec, class = myobj)
tm
tm[3:4] # still has attributes
On 10/11/06, Michael
On Thursday 05 October 2006 02:06, Roger Bivand wrote:
On Wed, 4 Oct 2006, Dylan Beaudette wrote:
Greetings:
As I am not a windows user, I cannot try this: is it possible to install
rgdal on windows without having to compile it from source ?
Andy Jaworski already replied that rgdal
In qqmath, how would one go about having 'distribution' change with
panel.number? I've tried
set.seed(1)
mydata - data.frame(ind = factor(rep(2:4, each = 100)))
mydata$val - rt(300, df=rep(2:4, each = 100))
plot-qqmath(~ val | ind,
layout=c(3,1),
data = mydata,
On 10/11/06, Benjamin Tyner [EMAIL PROTECTED] wrote:
In qqmath, how would one go about having 'distribution' change with
panel.number? I've tried
set.seed(1)
mydata - data.frame(ind = factor(rep(2:4, each = 100)))
mydata$val - rt(300, df=rep(2:4, each = 100))
plot-qqmath(~ val | ind,
On Wed, 2006-10-11 at 14:25 -0400, Brian Frappier wrote:
I tried all of the approaches below.
the problem with:
x - data.frame(matrix(NA,100,3))
for (i in 2:ncol(DF)) x[,i-1] - sample(rep(DF[,1], DF[,i]),100)
if you want result in data frame
or
x-vector(list, 3)
for (i in
Duncan Murdoch's definition is _the_ only one that I know. X is Uniform on A
means E phi(X) = \int_A phi(x) dx / \int_A dx, so that the probability
density is equal to 1/ \int_A dx everwhere on the set A.
By the way, another way to simulate X ~ Dirichlet(A1, A2, ..., Ad)
is to generate d
I don't have the previous messages, but it seems to me that
the solutions didn't quite get the requirements of the problem.
For example, it's obvious that for n = 2, the
random variables X1 and X2 should be X1 - runif(1)
and X2 - 1 - X1; while the solution X - runif(2);
X1 - X[1] / sum(X); X2 -
Thanks Manuel,
My only problems with the approach you suggest is that it does not seem to
result in a random sample without replacement as it generates a sample based
on the a priori probabilities, not physical selection and deletion from
subsequent sampling. I beleive the sample function would
On 10/11/2006 5:04 PM, Alberto Monteiro wrote:
I don't have the previous messages, but it seems to me that
the solutions didn't quite get the requirements of the problem.
For example, it's obvious that for n = 2, the
random variables X1 and X2 should be X1 - runif(1)
and X2 - 1 - X1; while
Thanks for the info. As it turned out i had to update a lot of packages. Finally
the configure worked and created the Makefile. But my make flopped at the end
with this error:
Any help is appreciated.
Thank you
-
/usr/bin/ld:
Dear Jerome,
Sorry, I'm not running your configuration. I cannot reproduce this
problem with Sybase Adaptive Server IQ 12.6 and R 2.3.1 + RODBC 1.1.7 on
WinXP SP2 or Linux 2.4.21-32.0.1.ELsmp #1 SMP Tue May 17 17:52:23 EDT
2005 i686 i686 i386 GNU/Linux.
-- ASIQ uses STRING() to mean CONCAT()
--
I'm not sure if this is the place to post this question, but, I am
having trouble compiling the source code. I do have a suitable C
compiler and f2c but I get this error when I run ./configure
configure: error: --with-readline=yes (default) and headers/libs are
not available
Any ideas? Thanks.
Michael Toews wrote:
Hi,
I have questions about object attributes, and how they are handled when
subsetted. My examples will use:
tm - (1:10)/10
ds - (1:10)^2
attr(tm,units) - sec
attr(ds,units) - cm
dat - data.frame(tm=tm,ds=ds)
attr(dat,id) - test1
When a primitive class object
On 10/11/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:
You can define your own class then define [ to act any way you would like:
[.myobj - function(x, ...) {
y - unclass(x)[...]
Careful, this is not the same as
y - NextMethod([, x, ...)
/Henrik
attributes(y)
(redirected to r-help)
On 11 October 2006 at 17:52, James Bullard wrote:
| Gavin Simpson wrote:
| On Wed, 2006-10-11 at 17:58 -0400, T C wrote:
| I'm not sure if this is the place to post this question, but, I am
| having trouble compiling the source code. I do have a suitable C
| compiler
On Wed, 11 Oct 2006 18:00:12 -0400
Yeh, Richard C [EMAIL PROTECTED]
wrote:
Dear Jerome,
Sorry, I'm not running your configuration. I cannot
reproduce this
problem with Sybase Adaptive Server IQ 12.6 and R 2.3.1
+ RODBC 1.1.7 on
WinXP SP2 or Linux 2.4.21-32.0.1.ELsmp #1 SMP Tue May 17
I changed the path to include C:\Program Files\R\R-2.3.1\bin, but Rcmd check
[package] just returns an error:
'Rcmd' is not recognized as an internal or external command, operable
program or batch file.
Any ideas? I want to be able to compile this package on Windows.
[[alternative HTML
Don't use the pathname with embedded blanks
C:\Program Files\R\R-2.3.1\bin
Use the 8.3 equivalent
C:\Progra~1\R\R-2.3.1\bin
You can find the 8.3 name with
dir /x
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
x - c(0,7,8,14,15,120,242)
y - c(122,128,130,158,110,110,92)
lowess(x,y)
$x
[1] 0 7 8 14 15 120 242
$y
[1] 122. 128. 132.2857 158. 110. -4930.
110.
R version 2.2.1, 2005-12-20, i486-pc-linux-gnu
attached base packages:
[1] methods stats
That's a windows message which says it can't find the command
you typed anywhere on its path.
If you can't figure it out get Rcmd.bat from batchfiles:
http://cran.r-project.org/contrib/extra/batchfiles/
and place that file anywhere on your path. It will
find R in the registry and run
Hello,
I have lots of data in zoo format and would like to do some time
series analysis. (using library(zoo), library(ts) )
My data is usually from one year, and I try for example stl() to find
some seasonalities or trends.
I have now accepted, that I might have to convert my series
75 matches
Mail list logo