klebyn a écrit :
my objective:
1) to save PNG files;
- i don't know the best way to make this;
?png
(also bmp(), jpg() available)
hih
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Claire Lee wrote:
I'm using R in Windows XP. I created a package myself.
I've used R CMD check to check it. Everything seems OK
except the latex. I get the error message:
* checking bbHist-manual.tex ... ERROR
LaTeX errors when creating DVI version.
This typically indicates Rd problems.
Hi Everybody,
i again have the problem in using tune.knn(), its giving an error saying
missing values are not allowed again here is the script for
BreastCancer Data,
library(e1071)
library(mda)
trdata-data.frame(train,row.names=NULL)
attach(trdata)
xtr - subset(trdata, select = -Class)
Paul Baer [EMAIL PROTECTED] writes:
Is it possible to write functions in such a way that, rather than
having to write a=function(a), one can just write function(a) and
have the variable passed as the argument be modified?
My real interest here is being able to invoke the editor by writing
Dear list members,
I am working with a Cox ph model for the duration of unemployment. The event of
interest
in my analysis is getting employed. I have various background variables
explaining this
event: age, sex, education etc. I have multiple unemployment spells per person.
I use a model with
Marc Schwartz wrote:
On Wed, 2005-10-05 at 22:19 +1000, sosman wrote:
I have some data in a CSV file:
time,pos,t,tl
15:23:44:350,M1_01,4511,1127
15:23:44:350,M1_02,4514,1128
15:23:44:350,M1_03,4503,1125
...
15:23:44:491,M2_01,4500,1125
15:23:44:491,M2_02,4496,1124
Hi,
As part of a batch process I run R with option --gui none.
This works nicely on one maching running
Version 1.9.0 Under development (unstable) (2004-03-06)
but on another machine using the newer Version 2.1.0 (2005-04-18) I
get the following fatal error message:
ERROR: unknown GUI none
Version 2.x no longer has a -gui=none option. Running BATCH jobs
supposedly disables a gui though.
Luc Vereecken
At 10:43 AM 10/6/2005, Morten Lindow wrote:
Hi,
As part of a batch process I run R with option --gui none.
This works nicely on one maching running
Version 1.9.0 Under development
On Thu, 6 Oct 2005, Luc Vereecken wrote:
Version 2.x no longer has a -gui=none option. Running BATCH jobs
supposedly disables a gui though.
There is no --gui=none option in 2.1.x (but there is in 2.0.x).
The gui is not disabled in 2.1.x either according to the help page nor in
practice
# Not sure why this did not work
#sapply(tw, setlen, len)
It probably did, but you discarded the result. Try
tw - sapply(tw, setlen, len)
On Thu, 6 Oct 2005, sosman wrote:
Marc Schwartz wrote:
On Wed, 2005-10-05 at 22:19 +1000, sosman wrote:
I have some data in a CSV file:
Hi all,
I am trying to compare two distance matrices with R. I would like to
create a XY plot of these matrices and do some linear regression on
it. But, I am a bit new to R, so i have a few questions (I searched in
the documentation with no success).
The first problem is loading a distance
you could first make the frames to put together in R, then outside of R
glue the frames into a gif.
eg:
frames-10
for(i in 1:frames) {
jpeg(paste(ani_, i, .jpg, sep = ))
plot(1:10,1:10, col = i)
dev.off()
}
then use an image editing program to glue the jpgs together -- eg
Can anybody help me write a code on the following data example, which
fills out all NA values by using a linear interpolation with the two
closest values?
Doy is day of year (%j).
Code example:
yr-c(rep(2000,14))
doy-c(16:29)
dat-c(3.2,NA,NA,NA,NA,NA,NA,5.1,NA,NA,NA,NA,NA,4.6)
you might want to see
Douglas Bates. Least squares calculations in R. R News,
4(1):17-20, June 2004. http://CRAN.R-project.org/doc/Rnews/
he gives some rules of thumb, eg
use solve(A,b) not solve(A) %*% b
use crossprod(X) not t(X) %*% X
use crossprod(X,y) not t(X) y
Katharine Mullen
Very great job
On 10/4/05, Tuszynski, Jaroslaw W. [EMAIL PROTECTED] wrote:
Hi,
I was playing with Mandelbrot sets and come up with the following code, I
thought I would share:
library(fields) # for tim.colors
library(caTools) # for write.gif
m = 400 # grid size
C = complex(
I've rolled up R-2.2.0.tar.gz a short while ago. This version contains
several changes and additions, mostly incremental. See the full list
of changes below.
You can get it from
http://cran.r-project.org/src/base/R-2/R-2.2.0.tar.gz
or wait for it to be mirrored at a CRAN site nearer to you.
Anyone know why I would get an Error: couldn't find function write.gif
despite loading library(caTools) with no errors in R 2.1.1 under XP?
Thanks,
Roger
library(fields) # for tim.colors
fields is loaded use help(fields) for an overview of this library
library(caTools) # for write.gif
Loading
The binary version of caTools on CRAN (1.0) does not have the write.gif
function but the source version (1.4) does ...
hth, ingmar
From: roger bos [EMAIL PROTECTED]
Reply-To: roger bos [EMAIL PROTECTED]
Date: Thu, 6 Oct 2005 08:14:42 -0400
To: Tuszynski, Jaroslaw W. [EMAIL PROTECTED]
Cc:
See write.gif function in caTools.
Jarek
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of klebyn
Sent: Wednesday, October 05, 2005 8:05 PM
To: r-help@stat.math.ethz.ch
Subject: [R] playing with R: make a animated GIF file...
Hello all
I am playing
Is this what you want?
yr-c(rep(2000,14))
doy-c(16:29)
dat-c(3.2,NA,NA,NA,NA,NA,NA,5.1,NA,NA,NA,NA,NA,4.6)
ta-cbind(yr,doy,dat)
ta
yr doy dat
[1,] 2000 16 3.2
[2,] 2000 17 NA
[3,] 2000 18 NA
[4,] 2000 19 NA
[5,] 2000 20 NA
[6,] 2000 21 NA
[7,] 2000 22 NA
[8,] 2000 23 5.1
[9,] 2000 24 NA
http://cran.r-project.org/src/contrib/Descriptions/caTools.html has current
versions (1.4) of both source and binary.
Jarek
\
Jarek Tuszynski, PhD. o / \
Science Applications International Corporation \__,|
Prof,
My apologies if you believe I was contradicting you. I though I was
concurring with you that 1) R does not have anything that is as powerful as
NuOpt and 2) there is no 'one' optimization package in R. I did, however,
want to point out that R does have separate packages that offer a fairy
Well, since knn() can't handle incomplete data as it says, you can
choose to either omit incomplete observations (e.g., using na.omit()),
or to impute the data if the conditions are met (missingness at random,
...); see, e.g., packages cat, mix, norm, and e1071 for that.
HTH,
David
Hi, hi all,
I am trying to compare two distance matrices with R. I would like to
create a XY plot of these matrices and do some linear regression on
it. But, I am a bit new to R, so i have a few questions (I searched in
the documentation with no success).
The first problem is loading a
I use the following code to keep my packages up to date. Is there some
reason why I am not getting full function availablility with this updating
code:
x - packageStatus(repositories=http://cran.r-project.org/src/contrib;)
st - x$avai[Status]
install.packages(rownames(st)[which(st$Status==not
Dear Peter,
See the effects package, described in
http://www.jstatsoft.org/counter.php?id=75url=v08/i15/effect-displays-revi
sed.pdf.
I hope this helps,
John
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
check out the effects package on CRAN.
Kjetil
Peter Dunn wrote:
Hi all
I'm doing some things with a colleague comparing different
sorts of models. My colleague has fitted a number of glms in
Genstat (which I have never used), while the glm I have
been using is only available for R.
Dear All,
I have written the following programs to find a non-singular (10*10)
covariance matrix.
Here is the program:
nitems - 10
x - array(rnorm(5*nitems,3,3), c(5,nitems))
sigma - t(x)%*%x
inverse - try(solve(sigma), TRUE)
while(inherits(inverse, try-error))
{
x -
Dear all,
I try for long to understand exactly what is the factor type and especially
how it works, but it seems too difficult for me
I read paragraphs about it, and I understand quite well what it is (I think)
but I still can't figure how to deal with.
Especially these 2 mysteries (for me) :
On 10/6/2005 9:14 AM, Florence Combes wrote:
Dear all,
I try for long to understand exactly what is the factor type and especially
how it works, but it seems too difficult for me
I read paragraphs about it, and I understand quite well what it is (I think)
but I still can't figure how to
Don't you have the dimension of x backward? Try:
set.seed(1)
x - matrix(rnorm(50, 3, 3), 10, 5)
vinv - solve(crossprod(x))
vinv
[,1] [,2] [,3] [,4] [,5]
[1,] 0.019918251 -0.006247646 0.006600209 0.003687249 -0.018670806
[2,] -0.006247646
'x' is a rank 5 matrix from which you want to create a rank 10
crossproduct! Try the following instead:
x - array(rnorm(100 * nitems, 3, 3), c(100, nitems))
sigma - crossprod(x)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public
I will be out of the office starting 10/05/2005 and will not return until
10/07/2005.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Well, cool,
why do you need to tell this to more than 3000 readers of R-help ??
Kouros I will be out of the office starting 10/05/2005 and will not
return until
Kouros 10/07/2005.
Kouros __
Kouros R-help@stat.math.ethz.ch mailing list
On Thu, 6 Oct 2005 16:10:15 +0200 Anette Nørgaard wrote:
This is exactly what I requested, thank you!! However I do actually
have several columns in my data sheet where I need to do the same
thing, then how do I come about that?
Look at na.approx() in package zoo.
Best,
Z
e.g.
Is doy intended to represent the number of days since the beginning
of the year? In that case convert the first two columns to class Date
and interpolate using approx. See ?approx for variations:
tt - as.Date(paste(yr, 1, 1, sep = -)) + doy - 1
ta[,dat] - approx(tt, dat, tt)$y
Even better
2d I can't manage to deal with factors, so when I have some, I
transform
them in vectors (with levels()), but I think I miss the power and
utility
of
the factor type ?
levels() is not the conversion you want.
in fact I use
'as.numeric(levels(f))[f]'
(from the ?factor
On 10/6/2005 10:20 AM, Florence Combes wrote:
2d I can't manage to deal with factors, so when I have some, I
transform
them in vectors (with levels()), but I think I miss the power and
utility
of
the factor type ?
levels() is not the conversion you want.
in fact I use
a last question, and thanks a million for your patience and your
explanations ...
I tried with a df called merged and a column named Pcc_0h_A (which is
numeric values):
length(as.vector(merged$Pcc_0h_A))
[1] 12202
as.numeric(as.vector(merged$Pcc_0h_A)[1:10])
[1] 12.276 11.958 14.098 13.843
na.approx(zoo(ta[,-seq(2)], tt))
where tt is as before.
On 10/6/05, Anette Nørgaard [EMAIL PROTECTED] wrote:
This is exactly what I requested, thank you!! However I do actually have
several columns in my data sheet where I need to do the same thing, then
how do I come about that?
e.g.
On 10/6/2005 10:50 AM, Florence Combes wrote:
a last question, and thanks a million for your patience and your
explanations ...
I tried with a df called merged and a column named Pcc_0h_A (which is
numeric values):
length(as.vector(merged$Pcc_0h_A))
[1] 12202
head(merged)
ID Name Pcc_0h_A Pcc_0h_swapped_A
3302 301495 Q0010_01 |Q0010||Hypothetical ORF 12.276 11.716
6943 309175 Q0010_01 |Q0010||Hypothetical ORF 11.958 11.271
14065 298935 Q0017_01 |Q0017||Hypothetical ORF 14.098 13.122
6420 306615 Q0017_01 |Q0017||Hypothetical ORF 13.843 13.061
5066
Marc == Marc Bernard [EMAIL PROTECTED]
on Thu, 6 Oct 2005 15:11:09 +0200 (CEST) writes:
...
Marc Here is the program:
Marc nitems - 10
Marc x - array(rnorm(5*nitems,3,3), c(5,nitems))
Marc sigma - t(x)%*%x
Marc inverse - try(solve(sigma), TRUE)
Hi all user R,
My simple question is...I have a vector of names of predictors,
text-c(datem,cola,eslom)...I try to plot the model with this predictor
in sequence loop,
for(i in 1:3){
png(paste(fig_,i,sep=))
plot(preplot.gam(mod9)[[i]],se=T,rug=F,main=,xaxt=n,ylab=,xlab=)
From: Graham Williams
Received Tue 04 Oct 2005 7:26pm +1000 from Karin Lagesen:
I have read and reread the boxplot and the boxplot stats page, and I
still cannot understand how and what boxplot shows. I realize that
this might be due to me not knowing enough statistics, but anyway...
See ?get.
Andy
From: Fernando Espíndola
Hi all user R,
My simple question is...I have a vector of names of predictors,
text-c(datem,cola,eslom)...I try to plot the model
with this predictor in sequence loop,
for(i in 1:3){
png(paste(fig_,i,sep=))
I do have full access to that directory. I have the
bbHist package in c:/PROGRA~1/R/rw2011/library/bbHist
directory. Then under the library directory I did
check and build. Here's what I got:
$ R CMD check bbHist
* checking for working latex ... OK
* using log directory
Claire Lee wrote:
I do have full access to that directory. I have the
bbHist package in c:/PROGRA~1/R/rw2011/library/bbHist
directory. Then under the library directory I did
check and build. Here's what I got:
Now I understand:
You *must not* have the sources in that library!
You want to
bady == bady [EMAIL PROTECTED]
on Thu, 06 Oct 2005 14:39:27 +0200 writes:
bady Hi, hi all,
I am trying to compare two distance matrices with R. I would like to
create a XY plot of these matrices and do some linear regression on
it. But, I am a bit new to R, so i have a
Can someone, please, explain the difference is results below (notice
the isdst value)
unlist(as.POSIXlt('2005-7-1'))
sec min hour mday mon year wday yday isdst
0 0 0 1 6 105 5 181 1
unlist(as.POSIXlt(as.Date('2005-7-1')))
sec min hour mday mon
On Thu, 6 Oct 2005, Pyy-Martikainen Marjo wrote:
I specify a model with strata by covariate interactions. I would like to
conduct a Wald test
for the null hypothesis no differences between any covariate effects in the 3
groups.
The survey package has a function regTermTest (which isn't
See the Avoiding Errors section of the Help Desk article
in R News 4/1.
On 10/6/05, Omar Lakkis [EMAIL PROTECTED] wrote:
Can someone, please, explain the difference is results below (notice
the isdst value)
unlist(as.POSIXlt('2005-7-1'))
sec min hour mday mon year wday yday isdst
Dear R-help,
When R starts to execute some code, all the graphics windows
associated with that R session go blank, and remain blank until the R
prompt returns.
I am using the ion window manager
(http://modeemi.cs.tut.fi/~tuomov/ion/) under Debian linux.
I notice that under the gnome desktop
On Thu, 6 Oct 2005, Omar Lakkis wrote:
Can someone, please, explain the difference is results below (notice
the isdst value)
Timezones. (That is carefully documented on the help page for
as.POSIXlt.)
unlist(as.POSIXlt('2005-7-1'))
sec min hour mday mon year wday yday isdst
0
Your R is an obselete version. There are some changes in 2.2.0 which
affect the hinting given to X11 WMs, so please try the current version of
R. (I know it had only recently been released, but the posting guide
refers you to the development versions.)
I don't know if this will solve it (I
I am running R 2.1.1 and R Cocoa GUI 1.12(1622) on a Mac Dual G5 with
OS 10.4.2
Anytime I try to use the syntax color preference for the built-in
editor, R crashes within a few minutes of working on a file. What
information can I provide that might facilitate investigation of this?
Thank
A related comment - don't rely (too much) on boxplots. They show only
a few things, which may be limiting in many cases and completely
misleading in others. Here are a couple of suggestions for plots which
you may find more useful than the standard box plots:
- figure 3.27 from
Hi all,
I'm new to the list here, and I have what I think is a simple
question. Using the circular package, is there a way to plot the
mean and variance on top of a rose diagram or other plot of the data?
Thanks in advance...
- Jason
Jason Horn
Boston University Department of Biology
5
You are right,
in the help page of mat2dist, dist2mat we found:
mat2dist and dist2mat are local copies of as.dist and as.matrix.dist of mva.
now these functions are in stats. We have to check it.
Thanks.
Quoting Martin Maechler [EMAIL PROTECTED]:
bady == bady [EMAIL PROTECTED]
on
I got the 'hgu133a' for R 2.0.0, would it be useful ?
On 10/6/05, Francesca Buffa [EMAIL PROTECTED] wrote:
Dear all,
does anybody have, or know where could I find, a version of the hgu133b
library compatible with R 2.0.1?
(I need it sooner than I can update to the new R version)
thank you
Hi,
I have been using the lqs function with method='lms'. However the
results I get are a little different from the results noted by Rousseeuw
Leroy (Robust Regression and Outlier Detection) and I was wondering
how to use these results for outlier detection.
I'm using the stackloss dataset,
I keep getting this error when I try to use the sem package. I and
another person who has successfully used the sem package for similar
analysis (fMRI effective connectivity) cannot figure out what is
wrong with my code. I would appreciate any suggestions.
The error message:
Error in
Suzanne Witt wrote:
I keep getting this error when I try to use the sem package. I and
another person who has successfully used the sem package for similar
analysis (fMRI effective connectivity) cannot figure out what is
wrong with my code. I would appreciate any suggestions.
It is
I am trying to use sem to measure effective connectivity among four
brain regions. I have pasted the code that I am trying to run since
that seems easier than trying to come up with another example.
The input data is time series data taken from SPM; they are each
1x121 columns of numbers.
Dear Suzanne,
Take a look at your model specification:
p.ram
[,1][,2][,3]
[1,] LM1 - LSMA LM1 - LSMA NA
[2,] LSMA - RSMA LSMA - RSMA NA
[3,] RSMA - RM1 RSMA - RM1 NA
[4,] LSMA - LSMA LSMA - LSMA NA
[5,] RSMA - RSMA RSMA - RSMA NA
[6,] RM1 - RM1
Jason Horn wrote:
Hi all,
I'm new to the list here, and I have what I think is a simple
question. Using the circular package, is there a way to plot the
mean and variance on top of a rose diagram or other plot of the data?
The package seems to be very well designed:
I did not know
Dear all,
I am facing a problem that seems to me more tricky now that it did at first
sight.
I have a collection of data which consist in frequency distributions: 6
patches had been proposed to female insects (for oviposition), 3 of them
corresponding to one treatment (A), the other 3 to
Selon bogdan romocea [EMAIL PROTECTED]:
A related comment - don't rely (too much) on boxplots. They show only
a few things, which may be limiting in many cases and completely
misleading in others. Here are a couple of suggestions for plots which
you may find more useful than the standard box
Dear all
I have yet to find code for making and showing error bars. Is there some
idea that the simple presentation of confidence intervals is not good
statistics?
Any leads appreciated.
Chris Buddenhagen, Botany Department, Charles Darwin Research Station, Santa
Cruz,Galapagos. Mail: Charles
Chris,
computation of confidence intervals is left to the user. Showing them
can be achieved using the arrows() function.
I hope that this helps,
Andrew
On Thu, Oct 06, 2005 at 08:19:27AM -0600, Chris Buddenhagen wrote:
Dear all
I have yet to find code for making and showing error bars.
I haven't seen a reply to this, so I will offer a comment:
It looks to me like the correlation you question is the correlation
between the estimated intercept and slope of the gls regression line.
This is different from the correlation between ts.mar and ts.anr. To
Dear Dr. Hartley:
RSiteSearch(Inverse autocorrelation function) produced 11 hits,
none of which seemed to relate to your question. If you have a
reasonable algorithm for computing the IACF, it might not be difficult
to program in R. If you have compiled code, e.g., C++ or Fortran,
I'm using R 2.0.1 under Windows XP. I get the following message when I
run the code listed below. I don't seem to have any problems using the
function slice outside integrate.
Error in integrate(slice, 0, Inf, , , , , , , a, b) * delta :
non-numeric argument to binary operator
[ By
Murray Jorgensen [EMAIL PROTECTED] wrote:
I'm using R 2.0.1 under Windows XP. I get the following message when I
run the code listed below. I don't seem to have any problems using the
function slice outside integrate.
Error in integrate(slice, 0, Inf, , , , , , , a, b) * delta :
Andrew Robinson wrote:
Chris,
computation of confidence intervals is left to the user. Showing them
can be achieved using the arrows() function.
I hope that this helps,
Andrew
Also do ?xYplot after installing the Hmisc package. -Frank
On Thu, Oct 06, 2005 at 08:19:27AM -0600,
75 matches
Mail list logo