Thank you for the quick answer, but i was asking if some method could do
combinations where numbers of the same vector wont combine.
A little bit like 'expand.grid', but with the possibility to choose the
length for output combinations.
Cláudio
2010/4/23 Jorge Ivan Velez jorgeivanve...@gmail.com
Clever solution. I use paste alot and didn't think to use it for this
problem. Thanks!
--
View this message in context:
http://r.789695.n4.nabble.com/Linux-System-Function-tp2063058p2063113.html
Sent from the R help mailing list archive at Nabble.com.
I am new to R and I suspect my problem is easily solved, but I haven't
been able to figure it out without using loops. I am trying to
implement Blair Karniski's (1993) permutation test. I've included a
sample data frame below. This data frame represents the conditional
means (C1, C2) for 3
Hello!
I started using R today to produce a bar chart that I can use in LaTeX.
Here is what I have got so far:
pdf(bugfix_costs.pdf)
costs - c(1, 5, 10, 20, 50, 150)
barplot(costs,
ylab=Relative cost of a bugfix,
names.arg=c(Requirements, Design, Code,
Hi,
I'm currently trying to find/define a relationship between one dependent and
several independant variables.
The problem is that i cannot use the normal multiple regression/correlation
in Spss because the data is not normal distributed.
i calculated the spearman roh and Kendalls tau
Thanks for all of the replies. For me at least, it turns out I had to manually
set permissions in the bitmap folder,
and the three folders (jpeg-8a, libpng, and libtiff). Works fine now.
Alan-
- Original Message -
From: Duncan Murdoch murdoch.dun...@gmail.com
Date: Friday, April 23,
I am new to R and I suspect my problem is easily solved, but I haven't
been able to figure it out without using loops. I am trying to
implement Blair Karniski's (1993) permutation test. I've included a
sample data frame below. This data frame represents the conditional
means (C1, C2) for 3
Zerdna,
Please note that the CRAN version 3.12 is about
to be replaced by a new cluster of packages now on R-Forge; we consider the
new bigmemory = 4.0 to be stable and recommend you start using it
immediately. Please see http://www.bigmemory.org.
In your case, two comments:
(1) Your for()
On 04/23/2010 10:51 PM, Ariane C. Böhm wrote:
Hi guys of the R-Help-Team!br /br /First of all - you do a great job!br
/I've found a lot of your mails in the internet. So I thought it would be a
good idea to ask you a question about R.br /br /R is new to me, so sorry,
if the question is too
On 04/24/2010 02:52 AM, tamas barjak wrote:
Hi All!
I have 2 plain questions:
1.)
I know that very primitive question, but that to grant it, that the drawing
on the screen divided up onto which part draw
for example:
layout(matrix(1:4,ncol=2, byrow=T))
plot(x, y, ...)--- 1. screen
plot(y,
Thanks again David,
I've made the changes. The optimize() function works alright but it does
not give me the intersection. I suspect it must have to do with the
user-defined function which I am totally clueless and am inexperienced
with. Mind showing the way.
Thanks again and thanks for
On 04/24/2010 05:48 AM, Maxim wrote:
Hi,
I have a very simple question, but I'm obviously not able to solve the
problem on my own.
I have a data.frame like
sample(c(A,B,C),size=20,replace = T)-type
rnorm(20)-value
data.frame(ty=type,val=value)-test
There must be some built in
On Apr 24, 2010, at 8:10 AM, Muhammad Rahiz wrote:
Thanks again David,
I've made the changes. The optimize() function works alright but it
does not give me the intersection. I suspect it must have to do with
the user-defined function which I am totally clueless and am
inexperienced
Hi,
I wanted the DICE coefficient (similarity measure for binary variables)
to be calculated in R and found that the igraph package has the option
of similarity.dice to do this. But, for this command, the input object
should be an igraph object. But, I have a dataframe of columns
containing
rm(list=ls())
PLEASE, DON'T DO THAT. Or rather you can do it in your workspace but
don't post it. It's not fair to a person who may not read your code line by
line before pasting it into their workspace and having it wiped out. Do you
expect us to completely clear out our workspaces just
Dear list,
just to put it in a simple way:
i read.csv from csv file to create a gdata
then, create array gdata34
however, when making a loop for assigning gdata34[1,m]-gdata[m,4],
this is what happen
gdata[1,4] 's real value is 10354, however, the gdata34[1,4] turns to be 883
then i checked
On Sat, Apr 24, 2010 at 9:05 AM, hadley wickham h.wick...@gmail.com wrote:
rm(list=ls())
PLEASE, DON'T DO THAT. Or rather you can do it in your workspace but
don't post it. It's not fair to a person who may not read your code line by
line before pasting it into their workspace and having
On Apr 24, 2010, at 9:05 AM, hadley wickham wrote:
rm(list=ls())
PLEASE, DON'T DO THAT. Or rather you can do it in your
workspace but
don't post it. It's not fair to a person who may not read your code
line by
line before pasting it into their workspace and having it wiped
out. Do
Bootstrap? You could start here:
http://www.statmethods.net/advstats/bootstrapping.html
-Ista
On Sat, Apr 24, 2010 at 6:07 AM, pauls182 pau...@gmx.net wrote:
Hi,
I'm currently trying to find/define a relationship between one dependent and
several independant variables.
The problem is that
On Apr 24, 2010, at 5:44 AM, Sebastian Bergmann wrote:
Hello!
I started using R today to produce a bar chart that I can use in LaTeX.
Here is what I have got so far:
pdf(bugfix_costs.pdf)
costs - c(1, 5, 10, 20, 50, 150)
barplot(costs,
ylab=Relative cost of a bugfix,
Hi,
Let s be a dataframe.
s
A B C
0 0 1
1 0 1
1 0 1
0 0 1
1 0 1
0 1 1
0 1 1
0 1 1
0 0 1
tab1=table(s[,c(1,2)])
tab1
B
A 0 1
0 3 3
1 3 0
tab2=table(s[,c(1,3)])
tab2
C
A 1
0 6
1 3
The problem is I need
The reason why Hmisc is not available from CRAN can be seen on the
package check summaries. Particularly for Hmisc:
http://cran.r-project.org/web/checks/check_results_Hmisc.html
which fails everywehre now. The maintainer has been notified by
automated messages.
We do not ship binaries of
If clearing out your workspace destroys *any* work, then something is
seriously wrong with your workflow.
Yes, of course. Lets all post viruses to run on each others'
machines. That will teach those users who don't run antivirus and
backup software between each posting to r-help. Also if
Perhaps, true in some respects. I am still chiseling out work using
primitive editing tools. But it still takes several minutes to load the
objects I am working on into memory and then several minutes each to build
new models. The models still reside in memory, since I do not know any
method
On Apr 24, 2010, at 8:59 AM, Shubha Vishwanath Karanth wrote:
Hi,
Let s be a dataframe.
s
A B C
0 0 1
1 0 1
1 0 1
0 0 1
1 0 1
0 1 1
0 1 1
0 1 1
0 0 1
tab1=table(s[,c(1,2)])
tab1
B
A 0 1
0 3 3
1 3 0
Just wrote a small piece of code to do so. Thanks.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Shubha Vishwanath Karanth
Sent: Saturday, April 24, 2010 6:28 PM
To: r-help@r-project.org; r-h...@stat.math.ethz.ch
Cc:
On Sat, Apr 24, 2010 at 9:31 AM, Dong GUO cgu...@gmail.com wrote:
Dear list,
just to put it in a simple way:
i read.csv from csv file to create a gdata
then, create array gdata34
however, when making a loop for assigning gdata34[1,m]-gdata[m,4],
this is what happen
So you assigned (m ==
-- Forwarded message --
From: Dong GUO cgu...@gmail.com
Date: Sat, Apr 24, 2010 at 4:50 PM
Subject: Re: [R] assign value between different type: Double vs Integer
To: Sarah Goslee sarah.gos...@gmail.com
Thanks, Sarah.
the loop is all right. this is what I just experimented.
and
Ah, you didn't provide a complete description of
your problem. gdata is a factor. You need to
read about converting factors to numeric, or
perhaps reconsider your method of creating or
importing gdata.
as.numeric(somefactor) does not extract the
printed value, it extracts the *order of the
Jay, thanks a bunch. New package seems to work just fine and great
improvement in docs by the way:). I tried the same example, new version
deals with it smoothly. In terms of usefulness of my sample code -- sure i
am writing same stuff to disk many times with only one handle -- it was some
toy
Hi, having used both the gee pacakge and the geepack package, i am unsure of
how to interpret the results.
Here are the results from the geeglm function from the geepack package
gee2-geeglm(data$erythema~data$product, data = data, id=subject,
family=binomial, corstr=independence)
Warning
Am 24.04.2010 16:01, schrieb Marc Schwartz:
If you really want the values on top of the bars, use text() instead
of mtext(). You will also need to increase the max value of the y axis
to make room for the labels using the ylim argument in barplot():
costs - c(1, 5, 10, 20, 50, 150)
mp -
Ah, Dr. Who's got nothing on me!
Beg pardon, I'm recuperating from an operation and clearly more than
usually disoriented. I meant 2.10.1 and 2.11.0 of course.
Now that I'm here, courtesy of my time slip, what I'd really love to see
in R 10.0.0, or even 2.12.0, is the choice to have output
Dear list,
just to put it in a simple way:
i read.csv from csv file to create a gdata
then, create array gdata34
however, when making a loop for assigning gdata34[1,m]-gdata[m,4], this is
what happen
gdata[1,4] 's real value is 10354, however, the gdata34[1,4] turns to be 883
then i
If I understand it correctly
par(mfrow=c(2,2))
plot(x,y)
plot(y,z)
.
should work.
On Apr 24, 2010, at 8:11 AM, Jim Lemon wrote:
On 04/24/2010 02:52 AM, tamas barjak wrote:
Hi All!
I have 2 plain questions:
1.)
I know that very primitive question, but that to grant it, that the
On Apr 24, 2010, at 10:09 AM, hadley wickham wrote:
Perhaps, true in some respects. I am still chiseling out work using
primitive editing tools. But it still takes several minutes to load
the
objects I am working on into memory and then several minutes each
to build
new models. The models
Hello there,
I am trying to run an ANOVA model using a non-Standard F ratio. Imagine that
the treatments (treatments 1 2) are applied to the row not to individual
samples. Thus the row is the experimental unit. Therefore my error term in
my ANOVA table should be the error associated with with
On Apr 24, 2010, at 8:06 AM, chrish...@psyctc.org wrote:
Ah, Dr. Who's got nothing on me!
Beg pardon, I'm recuperating from an operation and clearly more than
usually disoriented. I meant 2.10.1 and 2.11.0 of course.
Now that I'm here, courtesy of my time slip, what I'd really love to
see
Hi:
I need help with a dataframe(see pic attached). is a mix of dates and text.
I want to create a table either using latex function from hmisc or xtable. I
already know how to do this but the problem is getting the dataframe into R.
I don't have a reproducible example but I am hoping that the pic
Well, this has seriously gotten off the original topic.
While Hadley makes some sense, it is nevertheless
sometimes the case (surely so for David, I would surmise)
that one is putting together a response to an R-help
query when a new query prompts one to temporarily abandon
the first and
David Winsemius sent the following at 24/04/2010 18:24:
Now that I'm here, courtesy of my time slip, what I'd really love to see
in R 10.0.0, or even 2.12.0, is the choice to have output either in the
current plain text or to some simple formatting primitives that would
have default tabs
Hello,
I'm trying currently using barplot to summarize precipitation data. So far
I've compiled total annual snow and rain accumulation in a table
(attached). I've been successful at plotting it using the following code:
fig=barplot(t(Annual_Precip_table), horiz=TRUE, space=0, col=c(grey70,
Hi:
I need help with a dataframe. is a mix of dates and text.
I want to create a table either using latex function from hmisc or xtable. I
already know how to do this but the problem is getting the dataframe into R.
I don't have a reproducible example . If someone is interested in
helping with
Hi,
I'm having trouble seeing the added value over functions defined by setGeneric
vis-a-vis
methods defined by inheritance and polymorphism. setGeneric offers a 'clean'
call to a generic function, ie. no need to call new(), so less typing to do for
the user. But such explicit calls can also
Define a function dejunkify which removes commas and parens and
converts cols 2 onwards to numeric. Assuming that the lines near the
end of the file all start with B use B as the comment character and
read in the file and then dejunkify it. Then read in the lines again
and extract out the lines
On Sat, Apr 24, 2010 at 12:54 PM, Peter Ehlers ehl...@ucalgary.ca wrote:
Well, this has seriously gotten off the original topic.
While Hadley makes some sense, it is nevertheless
sometimes the case (surely so for David, I would surmise)
that one is putting together a response to an R-help
Hi,
I've got a csv file with scores like this:
40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,32,38,40,20,40,26,40,30,40,5,5,13,12,40,40,3,33,29,23,2,24,9,15,4,21,16,5,26,8,8,18
I used it but for example:
robots2[74]
[1] \Teoma\ or \Ask Jeeves\ or \Jeeves\
It wasn't sorted.
Sebastian.
2010/4/14 Barry Rowlingson b.rowling...@lancaster.ac.uk:
On Wed, Apr 14, 2010 at 6:20 PM, Sebastian Kruk residuo.so...@gmail.com
wrote:
I have a problem, In a few cases
quote Well it sounds like you won't lose any work, but you will lose
some time. But why not open another session of R for R-help questions?
Then you're never in any danger. I often have four or five instances of
R running for different projects.
Hadley
/quote
Well for one thing you turkey,
Hi
I do not speak in English very well, I'm sorry.
I try good:
layout(matrix(1:4,ncol=2, byrow=T)) = par(mfrow=c(2, 2))
plot(1,1) --- anddoes plot in the upper left corner
plot(2, 2) --- upper right corner
BUT, I like plot( n, m) for example lower right corner, or lower left
corner, or again
On Apr 24, 2010, at 3:53 PM, Albert Vilella wrote:
Hi,
I've got a csv file with scores like this:
40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,32,38,40,20,40,26,40,30,40,5,5,13,12,40,40,3,33,29,23,2,24,9,15,4,21,16,5,26,8,8,18
Hi,
I also want to show user's name at the lower-left corner. So I added one
more line to Deepayan's code but it didn't work. I am not sure how to use
the page argument here. Appreciate any comment.
lattice.options(default.args = list(page = function(n) {
panel.text(lab = sprintf(%s, date()),
You can use fBasics package
sample(c(A,B,C),size=20,replace = T) - type
rnorm(20) - value
data.frame(ty=type,val=value) - test
require(fBasics)
nam - rownames(basicStats(test$val))
result - do.call(cbind, with(test, tapply(val, ty, basicStats)))
rownames(result) - nam
result
Bests.
-
On Apr 24, 2010, at 4:52 PM, Jun Shen wrote:
Hi,
I also want to show user's name at the lower-left corner. So I added
one more line to Deepayan's code but it didn't work. I am not sure
how to use the page argument here. Appreciate any comment.
lattice.options(default.args = list(page =
Hello,
I am trying to calculate predicted values derived from one dataset into a
hypothetical dataset. I tried this line of code:
graphdata$fmgpredvalues - predict(Acs250.3.4, graphdata)
and received the following error message:
ERROR: ZXend[1], drop = FALSE] %*%lmeFit$beta
I
Identical Problem on Identical Platform. Only difference is R-2.11.0. The
trailing content of stats-Ex.Rout.fail is:
-
trees - window(treering, start=0)
(fit - StructTS(trees, type = level))
Call:
StructTS(x = trees, type = level)
Variances:
levelepsilon
0.0003700 0.0719877
On Apr 23, 2010, at 3:48 PM, Maxim wrote:
I have a very simple question, but I'm obviously not able to solve the
problem on my own.
I have a data.frame like
sample(c(A,B,C),size=20,replace = T)-type
rnorm(20)-value
data.frame(ty=type,val=value)-test
There must be some built in functions,
Hello,
I am new to R and am having difficulty formatting numerical output from a
regression analysis. My code iteratively performs linear regression on a
dataset while excluding certain data ranges.
My code:
rm(list = ls(all = TRUE))
sink(outfile)
dat - read.table(testdat, sep=\t, header=TRUE)
Hi, David,
Thanks. It is working perfectly now. But I still don't understand the
argument page, I couldn't find it in text() or par(). Also in the page =
function(n){..} the argument n does not appear anywhere and it seems
OK to be replaced with any other letter. Can anyone explain it a
Hello
After updating to R version 2.11.0, the xlim argument of the lattice
barchart function does no longer work as expected. For example, the
following code does not produce the same result in lattice lattice_0.18-3
and lattice_0.18-5
On 2010-04-24 13:50, hadley wickham wrote:
On Sat, Apr 24, 2010 at 12:54 PM, Peter Ehlersehl...@ucalgary.ca wrote:
Well, this has seriously gotten off the original topic.
While Hadley makes some sense, it is nevertheless
sometimes the case (surely so for David, I would surmise)
that one is
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Peter Ehlers
Sent: Saturday, April 24, 2010 3:51 PM
To: hadley wickham
Cc: r-help@r-project.org
Subject: Re: [R] Intersection for two curves
On 2010-04-24 13:50, hadley
On Apr 24, 2010, at 6:26 PM, Jun Shen wrote:
Hi, David,
Thanks. It is working perfectly now. But I still don't understand
the argument page, I couldn't find it in text() or par().
text and par are base graphics functions. page is an argument to
lattice functions:
(Page 26 of
==
y=c(9,9,17,11,7,8,15,5)
treat=c(A,C,B,C,D,A,B,D)
block=c(1,1,2,2,2,3,3,3)
q1.mod=aov(y~as.factor(treat)+as.factor(block))
q1.mod
summary(q1.mod)
par(ask=T)
plot(TukeyHSD(q1.mod))
===
On Apr 24, 2010, at 9:00 PM, David Winsemius wrote:
On Apr 24, 2010, at 6:26 PM, Jun Shen wrote:
Hi, David,
Thanks. It is working perfectly now. But I still don't understand
the argument page, I couldn't find it in text() or par().
text and par are base graphics functions. page is an
Many people seem to be reluctant to define functions,
even thought I think it is a pretty small step from
writing scripts to writing functions.
I'm not so sure - I find most students struggle to grasp that next
level of abstraction. Generalising from a specific task to a general
function is a
On 25 April 2010 08:41, Erich Studerus erich.stude...@bli.uzh.ch wrote:
Hello
After updating to R version 2.11.0, the xlim argument of the lattice
barchart function does no longer work as expected. For example, the
following code does not produce the same result in lattice lattice_0.18-3
and
Hi:
Here are a few options. The key step, it seems to me, is to transform the
data so that
you have one 'factor' consisting of the variable names and another column
with
the values; this is what the melt() function in the reshape package does.
This idea
should work irrespectively of the number of
Hello R folks,
I am encountering a problem with the following scatterplot function from the
car package:
scatterplot(y~x|z)
where y and x are continuous (interval) random variables and z is a
categorical variable. When z is a categorical variable coded 1 or 2, I
(appropriately) get a
I'm sure there are people willing to help you, but you
will have to be a little more forthcoming with information.
I don't know how many predict functions there are in the
2000-plus packages on CRAN, but there surely are quite a
few. Perhaps you could let us know what you're using
before you get
On Sat, Apr 24, 2010 at 7:43 PM, Felix Andrews fe...@nfrac.org wrote:
On 25 April 2010 08:41, Erich Studerus erich.stude...@bli.uzh.ch wrote:
Hello
After updating to R version 2.11.0, the xlim argument of the lattice
barchart function does no longer work as expected. For example, the
On Sat, Apr 24, 2010 at 3:26 PM, Jun Shen jun.shen...@gmail.com wrote:
Hi, David,
Thanks. It is working perfectly now. But I still don't understand the
argument page, I couldn't find it in text() or par(). Also in the page =
function(n){..} the argument n does not appear anywhere and it
On Fri, Apr 23, 2010 at 12:56 PM, Jon Zadra jr...@virginia.edu wrote:
I took a look at xYplot() but it is beyond over-complex, at least for my
level, and the documentation isn't too clear.
Is there a simple way to make functions like plotCI() and errbar() use the
lattice reference rather than
On 2010-04-24 21:30, Anthony Lopez wrote:
Hello R folks,
I am encountering a problem with the following scatterplot function from the
car package:
scatterplot(y~x|z)
where y and x are continuous (interval) random variables and z is a
categorical variable. When z is a categorical variable
On 04/25/2010 05:04 AM, Kevin Turner wrote:
Hello,
I'm trying currently using barplot to summarize precipitation data. So far
I've compiled total annual snow and rain accumulation in a table
(attached). I've been successful at plotting it using the following code:
If you are looping doing the assignment:
gdata34[1,m]-gdata[m,4]
why would yo expect that gdata[1,4] would have the same value as
gdata[1,4]. Do you have your indices reversed?
You need to provide at least the data in 'gdata' so we have a reproducible
example.
On Sat, Apr 24, 2010 at 9:13
You do not say what OS you are using and the differences in graphics
operations can be substantial. Had this been on a Mac, the answer
would be to page backward through the plots with cmd-left-arrow and
save them with cmd-s or File/Save As.. I suspect that similar
facilities exist on
77 matches
Mail list logo