Nguyen Dinh Nguyen wrote:
Dear all,
A friend of mine, who just installed Rcmdr package. When calling the
package, the error as following comes up. And actually, I've checked in CRAN
packages, the tcltk package doesn't exist anymore, but tcltk2.
Any help is appreciated
Regards
Nguyen D Nguyen
A more generic solution is
allComb - expand.grid(rep(list(letters[1:5]), 7))
whichComb - sapply(seq_len(ncol(allComb) - 1), x = allComb, function(i,
x){
whichCombination - sapply(seq(i + 1, ncol(x)), y = x, function(j,
y){
as.numeric(y[, i]) = as.numeric(y[, j])
})
Try to add options(stringsAsFactors = FALSE) in your Rprofile.site (in
the etc directory). Using as.is = TRUE seems safer than stringsAsFactors
= FALSE in the read.fwf function. Because as.is is set to FALSE by
default and stringsAsFactors is not set.
HTH,
Thierry
I am trying to debug a program, and I think tryCatch will help. The
functions involved
process through so many times before I encounter the error, things
are a bit slow to
use debug and browser().
I've read the help file and postings on conditions, and am still
having trouble.
Hi Yuan, Lucien, List.
try this:
f -
function (...)
{
args - list(...)
if(length(args)==0){
return(NULL)
}
if (length(args) == 1) {
return(args[[1]])
}
if (length(args) 2) {
jj - do.call(Recall, c(args[-1]))
return(do.call(Recall,
Hello everybody,
I don't understand the warning I get when executing the following code:
x - rnorm(100)
h - hist(x, plot= FALSE, freq = FALSE)
If I don't use the plot-argument, I don't get a warning...
Any explanations for me?
Antje
__
On 01-Sep-08 08:20:25, ONKELINX, Thierry wrote:
Try to add options(stringsAsFactors = FALSE) in your Rprofile.site
(in the etc directory). Using as.is = TRUE seems safer than
stringsAsFactors = FALSE in the read.fwf function. Because as.is
is set to FALSE by default and stringsAsFactors is
Dear R-colleagues,
this is a question from a R-newbie medical doctor:
I am evaluating data on inter-observer-reliability in endoscopy. 20
medical doctors judged 42 videos filling out a multiple choice survey
for each video. The overall-data is organized in a classical way:
observations (items
Dear Ted,
I noticed that as.is was set by default in read.fwf. So if the user sets
stringsAsFactor it is passed through ... to read.table. But I'm not sure
how as.is is passed to read.table when onlye stringsAsFactors is set. If
it's the default (FALSE) then it might be conflicting with
Hi Thierry,
On 01-Sep-08 09:45:27, ONKELINX, Thierry wrote:
Dear Ted,
I noticed that as.is was set by default in read.fwf. So if the
user sets stringsAsFactor it is passed through ... to read.table.
But I'm not sure how as.is is passed to read.table when onlye
stringsAsFactors is set. If
try something like this:
dat - read.table(textConnection(video 1 2 3 4 5 6 7 8 9 10 11 12 13
14 15 16 17 18 19 20 21
1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
3 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Dear Felix,
I have no idea about the calculation of your accuracy, sensitivity, etc, but
the sums:
dat - read.table(file=clipboard) #read in your data as dataframe dat
dat$comp - apply(dat,1,function(x) sum(x[-c(1,22)]==as.numeric(x[22])))
Good luck
Bart
drflxms wrote:
Dear R-colleagues,
It's a problem in read.fwf. It should not set a default for as.is, and in
R-devel will not.
On Mon, 1 Sep 2008, [EMAIL PROTECTED] wrote:
On 01-Sep-08 08:20:25, ONKELINX, Thierry wrote:
Try to add options(stringsAsFactors = FALSE) in your Rprofile.site
(in the etc directory). Using as.is =
I don't understand the warning I get when executing the following code:
x - rnorm(100)
h - hist(x, plot= FALSE, freq = FALSE)
The warning is
Warning message:
In hist.default(x, plot = FALSE, freq = FALSE) :
argument ‘freq’ is not made use of
The explanation is pretty simple.
freq
Hi All,
I would like to be able to use play3d() to translate the view along a
surface, and zoom in and out. All the examples in play3d() function
description are of rotations.
I would also like to know if there is a method of creating a legend of for
the 3d surface using rgl package.
Thanks in
On 01/09/2008 7:11 AM, Chibisi Chima-Okereke wrote:
Hi All,
I would like to be able to use play3d() to translate the view along a
surface, and zoom in and out. All the examples in play3d() function
description are of rotations.
I would also like to know if there is a method of creating a
Try this:
pairs - data.frame(pred = unlist(DF[2:21]), lab = DF[,22])
library(caret)
pred - factor(pairs$pred)
lab - factor(pairs$lab)
table(pred, lab)
sensitivity(pred, lab)
specificity(pred, lab)
Lines - video 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
1 1 0 0 0 0 0 0 0 0 0
Some junk got in at the beginning. It should be:
Lines - video 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
3 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
I am a new user of R package and I needed a help
regarding plotting confidence limits for FFT ( chisquare based) and spectrum
using the monte carlo simulations of background red noise . How do I do it
in R can any one help in this issue?
sudheer
[EMAIL PROTECTED] schrieb:
I don't understand the warning I get when executing the following code:
x - rnorm(100)
h - hist(x, plot= FALSE, freq = FALSE)
The warning is
Warning message:
In hist.default(x, plot = FALSE, freq = FALSE) :
argument ‘freq’ is not made use of
The explanation is
Date: Sun, 31 Aug 2008 21:29:38 -0400
From: ivo welch
Subject: Re: [R] SQL Primer for R
stumped again by SQL... If I have a table named main in an SQLite
data base, how do I get the names of all its columns? (I have a mysql
book that claims the SHOW command does this sort of thing, but it
Hi there,
I hope this question is not as stupid as the one before ...
I tried to shorten my histogram (because the distribution is quite skewed and I
simply don't want to see the long tail but still use the histogram plot). How
can I do something like this? (The example does not work but I
Antje,
Have a look at ?hist. It seems like you want this:
data - rnorm(100)
hist(data, xlim = c(0,6))
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and
You can do this:
hh - lapply(h, head, mh)
hh$breaks - head(h$breaks, mh + 1)
class(hh) - histogram
plot(hh)
On Mon, Sep 1, 2008 at 10:44 AM, Antje [EMAIL PROTECTED] wrote:
Hi there,
I hope this question is not as stupid as the one before ...
I tried to shorten my histogram (because the
I tried to shorten my histogram (because the distribution is quite skewed
and I simply don't want to see the long tail but still use the histogram
plot). How can I do something like this? (The example does not work but I
don't know why...)
data - rnorm(100) # as example, of course this is
I am a new user of R package and I needed a help
regarding plotting confidence limits for FFT ( chisquare based) and
spectrum
using the monte carlo simulations of background red noise . How do I do
it
in R can any one help in this issue?
sudheer
--
Henrique Dallazuanna schrieb:
You can do this:
hh - lapply(h, head, mh)
hh$breaks - head(h$breaks, mh + 1)
class(hh) - histogram
plot(hh)
Thanks a lot, that's exactly what I need :-) (I did not consider the class...)
Putting an xlim could also be an solution but I think this is much more
wow! the answer seems to be pragma table_info(main); thanks, Gabor.
__
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 provide commented,
Troels Ring wrote:
dear friends - I have an ordered factor, ID, in a data.frame, labs,
and make a table,
aa - with(labs,table(val 150,ID)) - and now the order in the table
is according to the alphabetic order in ID, and not as I want it. I
have tried adding
[order(unique(labs$ID)] but it
Hi
I am using mathematical anotation to make the labels in a barplot in
italics. as an example (adapted from
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/129808.html):
tt - c(italic(First~~set), italic(Second~~set), italic(third~~set)
barplot(1:3, names.arg=parse(text=tt))
But I would like
If I understand:
tt - c(expression(italic('First\nset')), expression(italic('Second\nset')),
expression(italic('third\nset')))
barplot(1:3, names.arg=tt)
On Mon, Sep 1, 2008 at 12:17 PM, Rainer M Krug [EMAIL PROTECTED] wrote:
Hi
I am using mathematical anotation to make the labels in a
On Mon, Sep 1, 2008 at 5:22 PM, Henrique Dallazuanna [EMAIL PROTECTED] wrote:
If I understand:
tt - c(expression(italic('First\nset')), expression(italic('Second\nset')),
expression(italic('third\nset')))
barplot(1:3, names.arg=tt)
Thanks a lot - that is exactly what I was looking for
On Mon, Sep 1, 2008 at 5:26 PM, Rainer M Krug [EMAIL PROTECTED] wrote:
On Mon, Sep 1, 2008 at 5:22 PM, Henrique Dallazuanna [EMAIL PROTECTED]
wrote:
If I understand:
tt - c(expression(italic('First\nset')), expression(italic('Second\nset')),
expression(italic('third\nset')))
barplot(1:3,
Dear R List,
I noticed something odd when I started saving R graphs in pdf format
and including them in LaTeX documents using (e.g.)
\includegraphics{mygraphic.pdf}. A very thin dotted border is
included in the LaTeX document around the graphic, and I cannot seem
to make it go away. I
dear friends - I have an ordered factor, ID, in a data.frame, labs, and
make a table,
aa - with(labs,table(val 150,ID)) - and now the order in the table is
according to the alphabetic order in ID, and not as I want it. I have
tried adding
[order(unique(labs$ID)] but it doesn't help.
I have R
How can I read a space-delimited file, where the data values for each case
are folded before column 80, and so appear on two lines for each case?
The first few cases look like this
loc type bio H2S sal Eh7 pH buf P K Ca Mg Na Mn Zn Cu NH4
OI DVEG 676 -610 33 -290 5.00 2.34 20.238 1441.67
I seem to be missing something here:
given a set X:{a,b,c,whatever...}
the mathematical definition of 'permutation' is the set of all possible
sequences of the elements of X.
The definition of 'combination' is all elements of 'permutation' which
cannot be re-ordered to become a different
Hello,
I am having trouble adding symbols to grid graphics. I am able to
create a lattice scatterplot using xyplot, which has a range from
-15:15 in both the X and Y directions. However, when I try to add
circles and text to this graph using grid.circle() and grid.text(),
they are offset using
Try this:
lines - read.table(
http://www.math.yorku.ca/SCS/viscollin/data/linthall.dat;, fill = TRUE,
header = TRUE)
DF - do.call(rbind,
lapply(split(lines, rep(1:(nrow(lines)/2), each = 2)),
function(x)cbind(x[1, 1:13], x[2, 1:4])))
names(DF) - names(lines)
On
Hello,
I would like to place some text that should appear in the same position
of the graphic device, preferably outside the plotting area, regardless
the x and y axes limits. How can I do this?
Best,
Dani
--
Daniel Valverde Saubí
Grup de Biologia Molecular de Llevats
Facultat de
Dani Valverde wrote:
Hello,
I would like to place some text that should appear in the same
position of the graphic device, preferably outside the plotting area,
regardless the x and y axes limits. How can I do this?
I think you can use title() with an adj= argument.
Otherwise, try mtext(),
I've thought about it, but the position should be always the same,
locator() should return me the same value.
Daniel Valverde Saubí
Grup de Biologia Molecular de Llevats
Facultat de Veterinària de la Universitat Autònoma de Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vallès- SPAIN
On 9/1/08, Jason Pare [EMAIL PROTECTED] wrote:
Hello,
I am having trouble adding symbols to grid graphics. I am able to
create a lattice scatterplot using xyplot, which has a range from
-15:15 in both the X and Y directions. However, when I try to add
circles and text to this graph using
Hi all,
the following problem is still beyond my R-knowledge:
I have one data vector containing the signal from 4 channels that are measured
subsequently and in repeating cycles (with one factor vector for cycle and
one for channel identification).
To extract the mean of each channel during
Dear Michael:
There are doubtless better ways to do this, but the I just got the
following ugly hack to work for me:
rawDat - readLines('linthall.dat')
(N - length(rawDat))
N2 - ((N-1)/2)
outDat - character(N2+1)
outDat[1] - rawDat[1]
for(i in 1:N2)
outDat[i+1] -
dear R wizards---I am not sure at what point I owe pennance for asking
so many questions. I am now wrangling debugging. I want to write a
function
assert = function( condition, ... ) {
if (!condition) {
cat(...); cat(\n);
browser();
}
stopifnot(condition);
}
assert( nrow(ds)==12,
please ignore part 1. of course, the cat works. my mistake. I just
need to learn how to step up frames, please.
regards,
/iaw
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
I posted before recently about fitting 3D data x, y, z where all have
error attached.
I want to predict z from x and y; something like
z = b0 + b1*x + b2*y
But multiple regression is not suitable because all of x, y, and z have errors.
I have plotted a 3D scatterplot of some data using rgl. I see
try:
tapply(data,factor-vectors,avg.f, part)
On Mon, Sep 1, 2008 at 1:41 PM, Wolfgang Waser [EMAIL PROTECTED] wrote:
Hi all,
the following problem is still beyond my R-knowledge:
I have one data vector containing the signal from 4 channels that are measured
subsequently and in repeating
Hi there,
Am I correct to believe that tetrachoric correlation is a special case
of polychoric correlation when there are only two levels to the ordered
factor? Thus it should be okay to use hetcor from the polycor package
to build a matrix of correlations for binary variables?
If this is
Hi,
Is there a way to launch a file from an R command?
To be more precise: I am using windows, and I have an html file in the
working directory of R, and I would like my code in R to automatically
launch the html file in the browser (that is, open the file). Basically, i
would like to simulate
On Mon, 01-Sep-2008 at 07:32PM +0200, Dani Valverde wrote:
I've thought about it, but the position should be always the same,
locator() should return me the same value.
Here's a function I wrote for this type of job many years ago. The xp
and yp paramaters are very similar to ncp units in
ivo welch wrote:
stumped again by SQL... If I have a table named main in an SQLite
data base, how do I get the names of all its columns? (I have a mysql
book that claims the SHOW command does this sort of thing, but it does
not seem to work on SQLite.)
If you are using RSQLite, then a lot
Hi
I am a new user of R, and doing clustering analysis (dendrogram). I am quite
interested in pvclust package because of facility of AU probability of each
cluster.
Now, i have a problem when handling large data set. I have data set of 183 rows
and 5546 columns, and interested to make
See ?browseURL
On Mon, Sep 1, 2008 at 5:06 PM, Applejus [EMAIL PROTECTED] wrote:
Hi,
Is there a way to launch a file from an R command?
To be more precise: I am using windows, and I have an html file in the
working directory of R, and I would like my code in R to automatically
launch the
Somewhere in between the R-Versions 2.6 and 2.7 the behaviour of the
function type.convert and therefore also read.table, read.csv, etc.
has changed (see below):
In 2.6 and before:
type.convert(c(+, -, +))
[1] + - +
Levels: + -
In 2.7 and later:
type.convert(c(+, -, +))
[1] 0 0 0
Michael Friendly wrote:
How can I read a space-delimited file, where the data values for each case
are folded before column 80, and so appear on two lines for each case?
The first few cases look like this
loc type bio H2S sal Eh7 pH buf P K Ca Mg Na Mn Zn Cu NH4
OI DVEG 676 -610 33 -290 5.00
ivo welch wrote:
please ignore part 1. of course, the cat works. my mistake. I just
need to learn how to step up frames, please.
parent.frame() will give you the frame of the caller. For example, in
the browser,
caller - parent.frame()
ls(caller)
caller$foo
etc. parent.frame() has
On Mon, 1 Sep 2008, Henrique Dallazuanna wrote:
See ?browseURL
Or more generally on Windows, ?shell.exec (which is what browseURL does by
default on Windows).
On Mon, Sep 1, 2008 at 5:06 PM, Applejus [EMAIL PROTECTED] wrote:
Hi,
Is there a way to launch a file from an R command?
To
On Mon, 1 Sep 2008, Tim Beissbarth wrote:
Somewhere in between the R-Versions 2.6 and 2.7 the behaviour of the function
type.convert and therefore also read.table, read.csv, etc. has changed (see
below):
In 2.6 and before:
type.convert(c(+, -, +))
[1] + - +
Levels: + -
In 2.7 and later:
Hi Duncan,
Yes, unfortunately my filter must have eaten your first response to my
parseHTTPHeader question. Thanks for responding to it again here on the
list and for the more efficient paste tip.
Below is a getURI call that returns a not found message. It is a
publicly available web site
Dear all,
I'm trying to interpolate a dataset to give it twice as many values (I'm giving
the dataset a finer resolution by interpolating from 1 degree to 0.5 degrees)
to match that of a corresponding dataset.
I have the data in both a data frame format (longitude column header values
along
Thank you prof,
Answered by the author of the function in question! :) Who could ask for
more?
Imad.
Prof Brian Ripley wrote:
On Mon, 1 Sep 2008, Henrique Dallazuanna wrote:
See ?browseURL
Or more generally on Windows, ?shell.exec (which is what browseURL does by
default on
On 01/09/2008 6:17 PM, Steve Murray wrote:
Dear all,
I'm trying to interpolate a dataset to give it twice as many values (I'm giving
the dataset a finer resolution by interpolating from 1 degree to 0.5 degrees)
to match that of a corresponding dataset.
I have the data in both a data frame
I've made an LDA model on some data from one source. I have some new data
that I want to see if I can place to the sources in the LDA model.
I used the predict function as follows:
predict(wak.insitu.ld, wak.alr.alluvial)
where wak.insitu.ld is an LDA model generated from some data and
I don't really understand this. You seem willing only to reveal the bit of the
process where you are convinced the problem is. This may not be helpful if the
problem lies elsewhere.
I take it you are using the lda() function from MASS, though you don't say.
Let's look at a little example:
Dear support, we are trying to export a frequency distributions obtanined by
stament freq(x,variable.labels=NULL,display.na=TRUE,levels=NULL)
of prettyR.
How can I do it?
Thanks in advance
__
R-help@r-project.org mailing list
I am estimating a multinomial model with two quantitative predictors, X1
and X2, and 3 responses. The responses are called neutral, positive and
negative with neutral being the baseline. There are actually many models
being estimated because I estimate the model over time and also for
various
Hi,
I am looking for R editor. does anyone know good editor? which tells you
syntax error and it has function to beautify format (insert TAB etc.).
Yu
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Hi Gents,
I need to get the control limits from qcc function.
As follows:
qcc(MDI, type = xbar.one)
Call:
qcc(data = MDI, type = xbar.one)
xbar.one chart for MDI
Summary of group statistics:
Min. 1st Qu. MedianMean 3rd Qu.Max.
0.3266 0.4249 0.4371
Check out:
http://www.sciviews.org/_rgui/projects/Editors.html
On Mon, Sep 1, 2008 at 8:52 PM, Yuan Jian [EMAIL PROTECTED] wrote:
Hi,
I am looking for R editor. does anyone know good editor? which tells you
syntax error and it has function to beautify format (insert TAB etc.).
Yu
Renata Mara Schneckenberg Wowk wrote:
Hi Gents,
I need to get the control limits from qcc function.
As follows:
qcc(MDI, type = xbar.one)
Call:
qcc(data = MDI, type = xbar.one)
xbar.one chart for MDI
Summary of group statistics:
Min. 1st Qu. MedianMean 3rd
The statement Globals - list() in the body of a function in a package was
intended to write an empty list to the R workspace to collect results during
the computations of the function. A package name space has not been specified.
The package appears to function correctly, but
during the R CMD
More and more I am getting warnings from packages that I install that the
package was built with 2.7.2 (I am running 2.7.1). I would like to upgrade but
don't want to loose all of the packages that I have installed and the settings.
Is there a way to just upgrade without uninstalling and
Hello, Kevin
You can get some hints by browsing in this mailist with the subject
of Upgrading R means I lose my packages, which were posted several
days ago.
HTH
Leon
[EMAIL PROTECTED] wrote:
More and more I am getting warnings from packages that I install that the package was
built
Hi R users,
I have one simple question but I really don't know why I can't get it work.
The data was adopted from Efron's An introduction to the bootstrap book.
nlaw
LSAT GPA
[1,] 576 3.39
[2,] 635 3.30
[3,] 558 2.81
[4,] 578 3.03
[5,] 666 3.44
[6,] 580 3.07
[7,] 555 3.00
if your two points are valid and you would be willing to run a least
spuares regression on then you should be fairly comfortable with
using na.approx in the zoo package. If interpolating values is going
to give you an erroneous result then maybe this is not the correct
path. But have a look.
Dear all,
I have a string (a file path in windows system) contains several \s.
When I copy and paste it into R for assignment to an variable, the \s
are always interpreted as escaping characters. How can I keep it intact?
Thanks.
Regards,
Leon
Please look at the help for boot, in particular what it says the form of
'statistic' should be. Possibly what you want is
corr.t - function(data, ind) cor(data[ind,])[1,2]
nlaw appears to be a single-column matrix: that will not work either.
On Mon, 1 Sep 2008, [EMAIL PROTECTED] wrote:
Hi
On Tue, 2 Sep 2008, Leon Yee wrote:
Dear all,
I have a string (a file path in windows system) contains several
\s. When I copy and paste it into R for assignment to an variable, the \s
are always interpreted as escaping characters. How can I keep it intact?
By reading from a file or
I do not see why you can not use regression even in this case.
To make things more simple suppose that the exact model is:
y = a + b*x, i.e.
y1 = a + b*x1
...
yn = a + b*xn
But you can not observe y and x. Instead you observe
ui = xi + ei (i=1,...,n) and
vi = yi + di (i=1,...,n)
Now you have
rkevinburton at charter.net writes:
I just realized after some tips and a little digging that what I was trying to
do manually has already been
done. I was trying to fit my data using 'lm' then taking the residual data
and trying to do a spectral
estimate (for seasonality) usiing fft and
Dear Renata,
In your particular case, something like this should do the job:
res=qcc(MDI, type = xbar.one)
names(res)
res$limits
A general approach (using the example in ?qcc) could be:
require(qcc)
data(pistonrings)
attach(pistonrings)
diameter - qcc.groups(diameter, sample)
ex =
Hi Professor Ripley,
(Thank you very much)^infinite for your response. I still have few questions
1. R shows that I have two column in nlaw data.
ncol(nlaw)
[1] 2
corr(nlaw)
[1] 0.7763745
2. Why can't I use boot::corr? I check corr in boot package and it
seems to be used for computing
84 matches
Mail list logo