Hi,
Try print(p) instead of plot(p)
HTH,
baptiste
2009/12/29 Bryan Hanson han...@depauw.edu:
I¹m trying to build a simple formula interface to work with a function using
ggplot2. The following scheme ³works² up until the plot(p) request, at which
point there are complaints about xlim¹s and
Hi,
I think you can also use plyr for this,
dft - read.table(textConnection(P1idVeg1Veg2AreaPoly2 P2ID
1 p p 1 1
1 p p 1.5 2
2 p p 2 3
2 p h 3.5 4),
On 12/28/2009 10:38 PM, Dean1 wrote:
How can I resolve this problem?...
As a general example,
plot (1:4)
polygon(c(0,0,5,5),c(0,5,5,0), border=lavenderblush1, col =
lavenderblush1)
###see how this overlays the axes lines
#I have tried...
for (k in 1:4) axis(k, lwd.ticks=0, label=F)
#...but
Nancy -
Please notice that ** is not an R operator.
The caret (^) is the exponentiation operator in R.
- Phil
On Tue, 29 Dec 2009, Nancy Adam wrote:
Hi everyone,
I tried to write the code of computing R2 for a regression system but I failed.
This
Dear All,
I am given 15 different data sets and I would like to generate a panel
showing all of them.
Each dataset will be presented either as a boxplot or as a histogram.
There are several possible ways to achieve this (as far as I know)
(1) using plot and mfrow()
(2) using lattice
(3) using
Hi,
Here is some artificial data followed by minimal ggplot2 and lattice examples,
makeUpData - function(){
data.frame(x=sample(letters[1:4], 100, repl=TRUE), y=rnorm(100))
}
datasets - replicate(15, makeUpData(), simplify=FALSE)
names(datasets) - paste(dataset, seq_along(datasets), sep=)
I forgot the base graphics way,
## divide the window in 4x4 cells
par(mfrow=n2mfrow(length(datasets)))
## loop over the list of datasets and plot each one
be.quiet - lapply(datasets, function(ii) boxplot(y~x, data=ii))
ggplot2 has a website with many examples,
http://had.co.nz/ggplot2/
as well
Nancy,
There is a function in the caret package called defaultSummary that
will compute both RMSE and R2 for you.
Max
On Dec 28, 2009, at 11:51 PM, Nancy Adam nancyada...@hotmail.com
wrote:
Hi everyone,
I tried to write the code of computing R2 for a regression system
but I
Hi:
To make a more legible ordering of the plots, it might be a good idea to
make m$L1 an
ordered factor rather than a character variable:
m$L1 - ordered(m$L1, levels = paste('dataset', 1:15, sep = ''))
# Lattice
bwplot(value ~ x | L1, data = m, as.table = TRUE)
# ggplot2
ggplot(m)+
Dear list,
I have a balance connected to the serial port of a windows machine (COM1) and
I read the text
output of the balance with
scan(COM1, what=character, sep=\n, n=1)
after calling the previous line I press the print key on the balance which
triggers sending one line
of text to the
Hi,
I currently combine multiple processed data (data frame) into a list and
save the list as .rda using the save command. When new data come, I load
the rda file, process the new data into a data frame, append the data
frame to the end of the list, and save the whole list to the disk. The
Hi,
I have a n*m matrix and would like to count the number of elements not equal to
NA in a ROW.
e.g.
x 1 2 3 NA 10
y 2 NA 8 9 NA
Which function can I use to obtain
4 for row x and
3 for row y?
Could you help me? I found some functions for columns but not for rows...
Thank you very much!
Try this:
rowSums(!is.na(x))
On Tue, Dec 29, 2009 at 11:49 AM, Verena Weber verenawe...@gmx.de wrote:
Hi,
I have a n*m matrix and would like to count the number of elements not equal
to NA in a ROW.
e.g.
x 1 2 3 NA 10
y 2 NA 8 9 NA
Which function can I use to obtain
4 for row x and
For a single row where mat is your matrix and r is the row
sum(!is.na(mat[r,]))
For every row in the matrix
rowSums(!is.na(mat))
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Verena Weber
Sent: Tuesday, December 29, 2009
On Sun, Dec 27, 2009 at 3:33 PM, Carlos J. Gil Bellosta
c...@datanalytics.com wrote:
I tried Amazon EC2 with R recently and wrote an entry about it to a blog I
collaborate with:
http://analisisydecision.es/probando-r-sobre-el-ec2-de-amazon/
(Unfortunately, it is in Spanish...)
Google's
after pressing up-arrow.
http://www.neilkodner.com/images/skitch/R_mystery-20091229-093157.jpg
Hopefully someone knows what to do about this.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
Jim Bouldin wrote:
When I try to run the following non-linear regression with variables
index1 and prl3:
beta = 4
nls(index1~beta*(1/prl3),start = list(beta = 4))
I get this error message:
Error in nls(index1 ~ beta * (1/prl3), start = list(beta = 4)) :
REAL() can only be applied to a
Hi Nancy,
Comments in line:
On Sun, Dec 27, 2009 at 3:34 AM, Nancy Adam nancyada...@hotmail.com wrote:
Hi everyone,
Can anyone please tell whether there is a difference between the code for
using svm in regression and code for using svm in classification?
This is my code for regression,
You can keep each dataframe as a separate file and process them that way.
You can look into storing in a relational database or using filehash. It
all depends on how you want to process the data later.
On Tue, Dec 29, 2009 at 8:40 AM, Hao Cen h...@andrew.cmu.edu wrote:
Hi,
I currently
** works as exponentiation for me:
3**2
[1] 9
On Tue, Dec 29, 2009 at 6:10 AM, Phil Spector spec...@stat.berkeley.edu wrote:
Nancy -
Please notice that ** is not an R operator. The caret (^) is the
exponentiation operator in R.
- Phil
On Tue, 29
Anny Huang annylhu...@gmail.com wrote:
Is there a way that I can import R functions into Fortran? Especially, I
want to generate random numbers from some not-so-common distributions (e.g.
inverted chi square) but did not find any routines written in Fortran that
deal with distributions other
Gabor Grothendieck wrote:
** works as exponentiation for me:
3**2
[1] 9
Yes, but see ?** that it is just documented in a note and has been
deprecated in S since 20 years...
Uwe Ligges
On Tue, Dec 29, 2009 at 6:10 AM, Phil Spector spec...@stat.berkeley.edu wrote:
Nancy -
Please
This problem may well be due to the repeated R_alloc calls, and might be
fixable by refactoring the code, but I am not sure that is the answer here.
This package is designed for the 3' biased arrays (and works well enough for
them), and isn't intended to be used for the newer random primer
Hi,
I’m new to R and object-oriented programming. I work primarily with
satellite datasets using IDL. I used IDL to ingest my data, run a
spectral mixture algorithm and output a timeseries of a snow index in
16-bit integer band interleaved by line format (for most of China). In
R, I
Hi,
I've been trying out the plm package, which seems like a great boon to
those who want to analyze panel data in R. I haven't started to use the
estimation functions themselves - for now I am just interested in having
a robust way to deal with lags in unbalanced panel data, since it is
On Tue, 29 Dec 2009, Phil Spector wrote:
Nancy -
Please notice that ** is not an R operator. The caret (^) is the
exponentiation operator in R.
- Phil
Actually, no,
2**3
[1] 8
Indeed ^ is the preferred exponentiation operator, but ** has 'always'
Dear All,
I am not sure how to remove double quotation marks in a string, e.g.,
paste(variable, 1). Can anybody please help me solve it? Thank you in
advance.
Lisa
--
View this message in context:
http://n4.nabble.com/Remove-double-quotation-marks-tp990502p990502.html
Sent from the R help
On 29/12/2009 1:16 PM, Lisa wrote:
Dear All,
I am not sure how to remove double quotation marks in a string, e.g.,
paste(variable, 1). Can anybody please help me solve it? Thank you in
advance.
I think you need to tell us what is wrong with what you get from that.
When I look at the result:
Dear R experts,
I would like to substitute the values of a vector, say A, with the average
values taken over a factor
For example, lets assume
A
[1] 1 2 3 4 5 6 7 8 9 10
B
[1] 1 2 2 4 4 4 7 7 7 7
Levels: 1 2 4 7
I need to have
1.0 2.5 2.5 5.0 5.0 5.0 8.5 8.5 8.5 8.5 8.5
Thank you
Thank you for your reply. But in the following case, “cat()” or “print()”
doesn’t work.
data.frame(cbind(variable 1, variable 2, cat(paste(variable, x), \n))),
where x is a random number generated by other R script.
Lisa
Duncan Murdoch wrote:
On 29/12/2009 1:16 PM, Lisa wrote:
Dear All,
have a look at ?ave(), e.g.,
A - 1:10
B - factor(c(1,2,2,4,4,4,7,7,7,7))
ave(A, B)
I hope it helps.
Best,
Dimitris
farida...@gmail.com wrote:
Dear R experts,
I would like to substitute the values of a vector, say A, with the average
values taken over a factor
For example, lets assume
Thanks! That's what I need.
On Tue, Dec 29, 2009 at 1:33 PM, Dimitris Rizopoulos
d.rizopou...@erasmusmc.nl wrote:
have a look at ?ave(), e.g.,
A - 1:10
B - factor(c(1,2,2,4,4,4,7,7,7,7))
ave(A, B)
I hope it helps.
Best,
Dimitris
farida...@gmail.com wrote:
Dear R experts,
I
Hi all,
I have two columns of data, the first consists of parameter estimates and
the second probability estimates. I would like to extract the value from
the first column based on the max value of column two, but not sure how to
do this other than by extracting the value and then manually
On Tue, Dec 29, 2009 at 6:31 PM, Lisa lisa...@gmail.com wrote:
Thank you for your reply. But in the following case, “cat()” or “print()”
doesn’t work.
data.frame(cbind(variable 1, variable 2, cat(paste(variable, x), \n))),
where x is a random number generated by other R script.
Lisa
Yes,
Will something like this work for you:
x - cbind(runif(10), runif(10))
x
[,1] [,2]
[1,] 0.26550866 0.2059746
[2,] 0.37212390 0.1765568
[3,] 0.57285336 0.6870228
[4,] 0.90820779 0.3841037
[5,] 0.20168193 0.7698414
[6,] 0.89838968 0.4976992
[7,] 0.94467527 0.7176185
[8,]
Hi,
I wonder how to pass several functions and their arguments as arguments to
a function. For example, the main function is
f = function(X ) {
process(X)
...
process(X)
}
I have a few functions that operate on X, e.g. g1(X, par1), g2(X, par2),
g3(X, par3). par1, par2 and par3 are
Hi,
I wonder how to pass several functions and their arguments as arguments to
a function. For example, the main function is
f = function(X ) {
process1(X)
...
process2(X)
}
I have a few functions that operate on X, e.g. g1(X, par1), g2(X, par2),
g3(X, par3). par1, par2 and par3
I am trying to convert a string to a double using as.numeric
However, R is truncating the results!
Options(digits) is set to 7.
Can anyone shed some light on this?
Thanks!
b[1]
[1] 116.28125
summary(b[1])
Length Class Mode
1 character character
c - as.numeric(b[1])
ticspd wrote:
I am trying to convert a string to a double using as.numeric
However, R is truncating the results!
Options(digits) is set to 7.
Can anyone shed some light on this?
Thanks!
b[1]
[1] 116.28125
summary(b[1])
Length Class Mode
1 character
Hello,
How do I fit a mixed effect model with no intercept using lmer()?
Is the following syntax correct?
lmer(y ~ -1+x1+x2+(-1+x1+x2|id), data=dt)
Thanks,
Julia
--
View this message in context:
http://n4.nabble.com/lmer-how-to-lmer-with-no-intercept-tp990493p990493.html
Sent from the R
Hello,
I have practice data of motor action in the format:
S | Cond. | Time
+-+
01 | c | 1.23
01 | nc| 0.89
02 | c | 2.15
02 | nc| 1.80
.
I want to look at the learning curves graphically.
I will appreciate pointers to relevant functions /
Could you just transpose the matrix?
Otherwise you can write a simple function that should work.
Try this
-(mat1 - matrix(c(1, 2, 3, NA, 10, 2, NA, 8, 9, NA),nrow=2))
gl - function(x)length(x[!is.na(x)]
apply(mat1, 1, gl)
?get
something like:
v1 - c(0, 1, 0)
v2 - c(1, 1, 0)
v3 - c(2, 1, 2)
v4 - c(2, 2, 1)
v5 - c(0, 1, 1)
x - 5
cbind(v1, v2, get(paste('v', x, sep='')))
v1 v2
[1,] 0 1 0
[2,] 1 1 1
[3,] 0 0 1
On Tue, Dec 29, 2009 at 3:54 PM, Lisa lisa...@gmail.com wrote:
Thank you for your
Thanks. I also need column name v5.
Lisa
jholtman wrote:
?get
something like:
v1 - c(0, 1, 0)
v2 - c(1, 1, 0)
v3 - c(2, 1, 2)
v4 - c(2, 2, 1)
v5 - c(0, 1, 1)
x - 5
cbind(v1, v2, get(paste('v', x, sep='')))
v1 v2
[1,] 0 1 0
[2,] 1 1 1
[3,] 0 0 1
On Tue, Dec
It seems from your example that you're assuming all of the vectors have
the same length. If this is the case, then a data.frame might be your
friend.
df - data.frame(
v1 = c(0, 1, 0),
v2 = c(1, 1, 0),
v3 = c(2, 1, 2),
v4 = c(2, 2, 1),
v5 = c(0, 1, 1) )
x - 5
get.var -
I had an NA in one row of my data frame, so I called na.omit(). But I do
not understand where that row disappeared to.
fri=na.omit(fri)
fri
Date.OnlyDAY Hour Min15 Quarter Arrival.Val Arrival4
109/05/2008 Friday833 3 328
210/24/2008 Friday 21
Your script works well. Thank you very mcuh.
Lisa
Nutter, Benjamin wrote:
It seems from your example that you're assuming all of the vectors have
the same length. If this is the case, then a data.frame might be your
friend.
df - data.frame(
v1 = c(0, 1, 0),
v2 = c(1, 1, 0),
On 29-Dec-09 21:11:38, James Rome wrote:
I had an NA in one row of my data frame, so I called na.omit().
But I do not understand where that row disappeared to.
fri=na.omit(fri)
fri
Date.OnlyDAY Hour Min15 Quarter Arrival.Val Arrival4
109/05/2008 Friday833 3
ah perfect .. digits = 8 is actually 8 digits TOTAL, not 8 digits after the
decimal point!
setting the default to 15 works for me
thank you!
Berend Hasselman wrote:
ticspd wrote:
I am trying to convert a string to a double using as.numeric
However, R is truncating the results!
This came in at 11:21 today:
From 30/12-2009 at 15.30 untill 31/12-2009 at 08.00, ALL IT services at
SUND-IT will be shut down!
ALL forms of electronic communication will be affected at SUND; VPN,
Webmail, Mail, Print, Network drives, Servers ie.
Furthermore, Panum, CSS, Teilum, Museion and
Hi Magnus,
Magnus Torfason zulutime.net at gmail.com writes:
I just noticed (the hard way of course) that when a query returns 0
rows, the columns in the resulting data.frame get dropped as well. See
the following example code (where conn is an active connection to an
SQLite db):
On 29/12/2009 3:08 PM, Hao Cen wrote:
Hi,
I wonder how to pass several functions and their arguments as arguments to
a function. For example, the main function is
f = function(X ) {
process1(X)
...
process2(X)
}
I have a few functions that operate on X, e.g. g1(X, par1), g2(X,
Thanks. I don't know how many functions there will be.
What if g1, g2, g3 have variable number of parameters? say g1 has parg1a,
parg1b, and g3 has parg3a, parg3b, parg3c, parg4d.
f - function(X, fs, pars) {
process1(X)
for (i in 1:length(fs))
fs[[i]](pars[[i]]) # doesn't work
Hello,
I am not able to plot a circle of a given radius using symbols(). In
the example below, the circle appears too large:
plot(0, 0, xlim = c(-1, 1), ylim = c(-1, 1))
symbols(0, 0, circles = 1, inches = FALSE, add = TRUE)
What's happening?
Ery
On 29/12/2009 5:24 PM, Hao Cen wrote:
Thanks. I don't know how many functions there will be.
What if g1, g2, g3 have variable number of parameters? say g1 has parg1a,
parg1b, and g3 has parg3a, parg3b, parg3c, parg4d.
pars could be a list of lists of parameters instead of a list of
Hi r-users,
I would like to use prop.test code and I also calculate manually to test the
proportions for 2 groups. The problem is the answer for the p-value calculated
manually are different from prop.test. Here are the results:
## Manually
z value: z= (phat-p)/sqrt(pq/n) =
On Tue, Dec 29, 2009 at 11:46 AM, Julia7 liujul...@yahoo.com wrote:
Hello,
How do I fit a mixed effect model with no intercept using lmer()?
Is the following syntax correct?
lmer(y ~ -1+x1+x2+(-1+x1+x2|id), data=dt)
Yes. An alternative, which I prefer, is
lmer(y ~ 0 + x1 + x2 + (0 + x1 +
The circle is not too large; it just uses your x-units for the
radius and your plotting region isn't square.
Here are some things to try:
plot(0, 0, xlim = c(-1, 1), ylim = c(-1, 1))
symbols(0, 0, circles = 1, inches = FALSE, add = TRUE)
abline(v = c(-1,1)) # shows that x-units 'fit'
abline(h =
You need to learn how to use prop.test properly.
Is this a ONE-sample test or TWO-sample test?
Perhaps this will help:
#hand calculation
phat - 58/691
p - 56/691
q - 1-p
n - 691
z - (phat-p)/sqrt(p*q/n)
z
[1] 0.2787970
pnorm(-z)*2
[1] 0.7804006
#prop.test
prop.test(x=58, n=691, p=56/691,
ticspd wrote:
I am trying to convert a string to a double using as.numeric
However, R is truncating the results!
No it isn't! As someone phrased it recently, there's a difference
between an object and the display of an object. Ceci n'est pas une pipe.
Options(digits) is set to 7.
So
Hi Hao,
I suggest you try again, starting by read posting guide at footnote of this
email.
How about a title for the message? How about identify yourself?
bests
milton
On Tue, Dec 29, 2009 at 2:59 PM, Hao Cen h...@andrew.cmu.edu wrote:
Hi,
I wonder how to pass several functions and their
Hi Hao,
Ok.
Sorry for my last post.
bests
milton
On Tue, Dec 29, 2009 at 3:08 PM, Hao Cen h...@andrew.cmu.edu wrote:
Hi,
I wonder how to pass several functions and their arguments as arguments to
a function. For example, the main function is
f = function(X ) {
process1(X)
...
Hi steve,
Thank you so much for your reply.Im asking about the difference between two
cases:1) when I use svm in a regression system and 2) when I use svm in a
classification system. Is the code of using svm in these two cases the
same?This is the code for a regression system:
my_svm_model -
Hi everyone,
Thanks a lot for the explanationâ¦
I tried the following code to compute R2 for a regression system but it does
not work:
my_svm_model - function(myformula, mydata, mytestdata)
{
mymodel - svm(myformula, data=mydata)
k- summary(mymodel)
k$r.squared
}
Can anyone please tell
Dear useR's
I have a small basic problem which I am hoping to get some help with. I have
a data frame, testSeq_df, with 1 row and 500 columns. Each column is a
character (a,c,g or t). I want this sequence to have 4 factors (a,c,g,t).
When I try the following:
for(i in 1:500){
if
Please ignore my message. I figured it out myself.
The way to do it is:
levels(testSeq_df[,i]) - list(a=c('a'),c=c('c'),g=('g'),t=c('t'))
Thanks,
vishal
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Hi,
I seem to have to some problem using the rJava package. Any help is appreciated.
I have a sample java file Hello.java
public class Hello
{
public String sayHello()
{
String result=new String(Hell);
return result;
}
public static void main(String[] args)
{
}
}
and these are the
Hi Nancy,
2009/12/30 Nancy Adam nancyada...@hotmail.com:
Hi steve,
Thank you so much for your reply.
I’m asking about the difference between two cases:
1) when I use svm in a regression system and
2) when I use svm in a classification system.
Is the code of using svm in these two
thank you.
and for those who didn't get the reference:
http://en.wikipedia.org/wiki/Ren%C3%A9_Magritte
http://en.wikipedia.org/wiki/Ren%C3%A9_Magritte
Peter Dalgaard wrote:
ticspd wrote:
I am trying to convert a string to a double using as.numeric
However, R is truncating the
69 matches
Mail list logo
