On Wed, 4 Mar 2009, David Winsemius wrote:
In what ways is rnorm not a satisfactory answer?
My guess was that CB wants to generate a finite population whose mean and
variance are specified, which would involve rnorm() followed by centering and
scaling.
-thomas
--
David Winsemius
Hello,
I have a data.frame called spec containing data about samples. But I don't
have these data for all my samples.
So if I have data (that is code of the sample is in spec$Code), I would like
to assign data1 to the variable m.
If I don't have this data, I would like to assign 1 to m.
I tried
glaporta glaporta at freeweb.org writes:
Dear list,
I read the SOM package manual but I don't understand how to perform (for
example) 1) the SOM analysis on Iris data 2) with a visualization similar to
that of figure 7 in
http://www.cis.hut.fi/projects/somtoolbox/package/papers/techrep.pdf
plot(x,rho,pch=id)
-Mensaje original-
De: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]en
nombre de Dipankar Basu
Enviado el: 03 March 2009 20:11
Para: r-help@r-project.org
Asunto: [R] scatter plot question
Hi R Users,
I have a dataframe like this:
id x rho
A
Hi
r-help-boun...@r-project.org napsal dne 04.03.2009 09:11:06:
Hello,
I have a data.frame called spec containing data about samples. But I
don't
have these data for all my samples.
So if I have data (that is code of the sample is in spec$Code), I would
like
to assign data1 to the
Dear Mohammad,
have a look at the finance task view on CRAN:
http://cran.at.r-project.org/web/views/Finance.html
(Dirk has nicely updated this page recently).
In addition, Patrick Burns provides a recipe for PC-GARCH models on his
web-site:
Hi
Read what barplot does and look to your plot.
If you want each row to be plotted as stacked bar with names uder each bar
taken from peaople variable then you need to
transpose your matrix - barplot takes columns
plot without names - you do not want them really plotted
add names under each
Dear Lists,
This is an off topic (OT ).
I recently took Anne Milley's interview .In Part 1 of the interview , Anne
talks about SAS, WPS, other softwares she studied like SPSS,.She also talks
about the difference between small and big companies , what sets SAS apart
and the famous licensing model
Dear Carlos and Kjetil,
Thanks for your answer.
I do not think that is the way to go. If you believe that your algorithm
is better than the existing one, talk to the author of the package and
discuss the improvement. The whole community will benefit.
I should be able to *easily* modify it and
ARDIA David david.ardia at unifr.ch writes:
Dear all,
I very often use the R function apply, for speedup purposes. I am now also
using MATLAB, and would like to use the same kind of function.
I have already asked MATLAB people, and the answer is : vectorize... but of
course, this is not
Vedula, Satyanarayana svedula at jhsph.edu writes:
I need to select one row per patient i in clinic j. The data is organized
similar to that shown below.
...
If patient has outcome recorded at visit 2, then outcome = outcome
columns at visit 2
If patient does not have visit 2, then
Dear Dave,
Another option would be to manually create the breaks of your plot.
library(ggplot2)
myBreaks - sort(c(pretty(mtcars$mpg), min(mtcars$mpg)))
ggplot(aes(x = wt, y = mpg), data = mtcars) + geom_point() +
geom_hline(yintercept = min(mtcars$mpg)) + scale_y_continuous(breaks = myBreaks)
matel.rtest accepts both Euclidean and non Euclidean matrices. If both
distance matrices are Euclidean, then a faster algorithm is used to
speed up the permutation procedure.
You can also use mantel.randtest which is exactly similar to rtest but
faster as computation are in C.
Sarah Goslee
Hi Guys,
I'm a numbie to R.
I searched the forum for plotting pdf's of functions but couldn't find one
that explained my question.
I have been asked to plot the pdf fX(x) = 1/pi(1+x^2). -infinity x
infinity.
I would be greatfull if someone might point me to a webpage that might give
me a
On 04-Mar-09 07:55:11, Bob Roberts wrote:
Hi,
I'm trying to write a regular expression that captures numbers in the
form 9,007,653,372,262.48 but does not capture dates in the form
09/30/2005
I have tried numerous expressions, but they all seem to return the
dates as well. Thanks.
Testing
Not clear what else is allowed to match and what not but this will
match strings of numbers commas and decimal points but not
match if there is anything else in the string like a slash:
regexpr(^[0-9.,]+$, 9,007,653,372,262.48)
On Wed, Mar 4, 2009 at 2:55 AM, Bob Roberts quagmire54...@yahoo.com
have a look at ?curve(), e.g.,
curve((1 + x^2) / pi, -15, 15)
I hope it helps.
Best,
Dimitris
Beetle wrote:
Hi Guys,
I'm a numbie to R.
I searched the forum for plotting pdf's of functions but couldn't find one
that explained my question.
I have been asked to plot the pdf fX(x) =
2009/3/4 Daniel Nordlund djnordl...@verizon.net:
Something like this may help get you started.
std.pop - function(x,mu,stdev){
((x-mean(x))/sd(x)*stdev)+mu
}
population - std.pop(rnorm(1000),10,5)
Hope this is helpful,
Dan
Very helpful. I hadn't thought of a simple roll-your-own
2009/3/4 Thomas Lumley tlum...@u.washington.edu:
On Wed, 4 Mar 2009, David Winsemius wrote:
In what ways is rnorm not a satisfactory answer?
My guess was that CB wants to generate a finite population whose mean and
variance are specified, which would involve rnorm() followed by centering
but this is not a probability density function, since lim(abs(x)-
infty)=infty. But there is just a pair of brackets missing in Beetles post.
btw. homework?
Dimitris Rizopoulos schrieb:
have a look at ?curve(), e.g.,
curve((1 + x^2) / pi, -15, 15)
I hope it helps.
Best,
Dimitris
Beetle
I feel really dumb for having to ask this, but here I go anyway...
I perform structural equation modeling on a survey, using about 25
variables that create a total of 5 latent variables (factors).
Applying sem (using the sem package) was a piece of cake, even for an
(SEM) layman, thanks for THE
Thank you.
I do not know where the source code is available.I never uploaded anything to
CRAN. I am just an R user so far.
Can you please tell me ?
I did notify the package maintainer mentioned in the package on-line
documentation, namely
Eric Aldrich eldr...@gmail.com, twice.
I suspect hic
Hello !!
II'm trying to test for my fixed effects using an lmer with quasipoisson
errors.
Since my lmer model is corrected for overdispersion using this kind of
errors, I should use during model simplification in my Anovas *F test *and
not *Chi square test* to compare two models. So I write:
Matthieu Stigler wrote:
Hello
I looked on the archives but did not find answer for that...
Running R CMD check for a package, i get an error:
Error in parse(n = -1, file = file) : unexpected symbol at
3341: }
But how can I find which file is guilty? What is this 3342 referring to?
Given
Please read the posting guide which asks you to answer basic questions
such as:
Which R / lattice versions are we talking about?
Which is the any format?
Are you using the Devices directly or are you using some other way to
copy contents of one device into another device?
What is the exact,
Hi, take a look on pls package and it's documentation, there are
examples also for NIR data.
http://mevik.net/work/software/pls.html
Article form Journal of Statistical Software
http://www.jstatsoft.org/v18/i02
Also Caret package can be used to evaluate pls and other regreesion models:
Matthieu Stigler wrote:
Hello
I looked on the archives but did not find answer for that...
Running R CMD check for a package, i get an error:
Error in parse(n = -1, file = file) : unexpected symbol at
3341: }
But how can I find which file is guilty? What is this 3342 referring to?
Finally
On Wed, Mar 4, 2009 at 2:48 AM, Daniel Nordlund djnordl...@verizon.net wrote:
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of CB
Sent: Tuesday, March 03, 2009 10:05 PM
To: David Winsemius
Cc: r-help@r-project.org
Subject: Re:
Dear all,
I wrote a function test1 in test1.R.
Right, I am writing another function test2 on test2.R and trying to use
test1 function.
How can I do?
Is there any similar way like including test1.R in test2.R file?
Thank You Very Much.
Jia-Ming
$BD%2HLC(B
Dear Sir or Madam,
I've been using R for one year and I really appreciate it.
I would like to know if a version performing parallel computations on multicore
computers and computer clusters exists.
Thank you very much.
Yours sincerely.
Enrico Giorgi
Jim Lemon a écrit :
Matthieu Stigler wrote:
Hello
I looked on the archives but did not find answer for that...
Running R CMD check for a package, i get an error:
Error in parse(n = -1, file = file) : unexpected symbol at
3341: }
But how can I find which file is guilty? What is this 3342
Many, many thanks for all the answers!
Notepad++ looks very promising although it does not have a project file
management facility. But it has a very clean appearance. I'll have to look
into SciTE which also sounds quite good. There seem to be some good
alternatives.
Meanwhile, I found a
To calculate Percentile for a set of observations Excel has percentile()
function. R function quantile() does the same thing. Is there any
significant difference btw percentile and quantile?
Regrads,
--
View this message in context:
That's great - thanks!
On 4 Mar 2009, at 08:59, ONKELINX, Thierry wrote:
Dear Dave,
Another option would be to manually create the breaks of your plot.
library(ggplot2)
myBreaks - sort(c(pretty(mtcars$mpg), min(mtcars$mpg)))
ggplot(aes(x = wt, y = mpg), data = mtcars) + geom_point() +
Hi,
I think you just have to call it with the correct arguments inside test2
function. So the arguments for test1 functions should be consider to be also
arguments of test2 function.
test2 - function( test1arguments, test2arguments ){
...
value - test1( test1arguments )
...
}
Paulo Ricardo Gherardi Hein phein1...@gmail.com writes:
I am new here (since jan2009) and up to now, I not seen anyone commenting
about principal component analysis and regression PLS to analyze spectral
information in R system. Sorry, I am a R starter...
Anybody have any package, or trick
Hi Paulo,
You might also want to look at something like the glmnet package (Friedman,
Hastie, and Tibshirani). This carries out penalized regression, is designed
to work with high numbers of predictors/inputs/columns and relatively few
samples/obervations/rows, and is very fast.
See:
Dear All,
1. Does anyone have experience of running R on a server inside a Citrix
shell - I'd like to get R onto the server and would be greatful for any
tips or direction on the matter.
2. This may seem like a silly question so forgive my ignornace.
Most of the data I currently work with is
Hi,
You can use the source() function to load functions from other files, like:
source(test1.R)
cheers,
Paul
Chang Jia-Ming wrote:
Dear all,
I wrote a function test1 in test1.R.
Right, I am writing another function test2 on test2.R and trying to use
test1 function.
How can I do?
Is
Hello
I am am new to R and any help with the following would be appreciated:
I have a list (example attached) and I would like to create a new list
which is a filtered version of this list. I.e I would like a list that
only contains elements with this value:
Chr 10 : 21853562 - 21855482
Any
Matthieu Stigler wrote:
Jim Lemon a écrit :
Matthieu Stigler wrote:
Hello
I looked on the archives but did not find answer for that...
Running R CMD check for a package, i get an error:
Error in parse(n = -1, file = file) : unexpected symbol at
3341: }
But how can I find which file is
On Tue, 03 Mar 2009 22:10:42 +0100, David Winsemius
dwinsem...@comcast.net wrote:
That is what I thought to be the critical paragraph. The variance is
assumed to be = 1 when you use family=gaussian rather than the default
of family=qgauss. You give it a vector, 1000*rnorm(100), that ranges
Hi
r-help-boun...@r-project.org napsal dne 04.03.2009 13:18:30:
Hello
I am am new to R and any help with the following would be appreciated:
I have a list (example attached) and I would like to create a new list
which is a filtered version of this list. I.e I would like a list that
only
John Fox wrote:
Dear Sebastian,
What you're looking for are factor-score coefficients, which would allow you
to estimate the values of the factors from the observed variables. Then,
given the original dataset from which the input-covariance matrix to sem()
was computed, you could find the
Dear Sebastian,
What you're looking for are factor-score coefficients, which would allow you
to estimate the values of the factors from the observed variables. Then,
given the original dataset from which the input-covariance matrix to sem()
was computed, you could find the factor scores. The sem
In lme4 I wamnt to analyse a dataset in which the random effects are
somewhat intricate. The nearest approximation of what I want to know is
lmer(hr~hpos+((bin==b)|VP)+(1|(VP:exp)), data)
In this design VP denotes the subject variable, and exp the experimental
unit (which is obviously nested in
Dear Nikol,
Try this:
do.call(c,lapply(yourlist,function(x) x[2] ))
HTH,
Jorge
On Wed, Mar 4, 2009 at 7:18 AM, Nikol Simecek ns...@cam.ac.uk wrote:
Hello
I am am new to R and any help with the following would be appreciated:
I have a list (example attached) and I would like to create a
Hello,
If test2 and test1 functions are in two files like test2.R and test1.R,
how could I do to let test2 know test1 function?
Cheers,
Jia Ming
2009/3/4 Usuario R r.user.sp...@gmail.com
Hi,
I think you just have to call it with the correct arguments inside test2
function. So the
I usually use the ?source function.
so at the top of your test2.R file, you would put a line like:
source('C:\\test1.R')
but with the correct path to where ever you have stored your test1.R
file.
Hope that helps a little,
Tony
On 4 Mar, 10:29, Chang Jia-Ming chang.jiam...@crg.es wrote:
Dear
Hello,
I need to insert a line after every eigth row (group/suset) which should
contains the following:
(an incremented ID),0,1,1,1,1
I have two problems with constructing my code:
1. I am getting NAs in three of the columns and
2) I can not find a way to write the ID (group name) at the
Dear list,
when running a mood.test() (part of package stats) on slightly longer
vectors (than the example from the help-page) we get the error-message
shown below : once both vectors tested are of length 50 this error oocurs.
Note, that this problem didn't occur with R-2.7.x (or even older
Hello
My question might sound awkward, but I am looking for a way to somehow
convert a plot in R into a base64 string.
Here's an idea, but it is not at all satisfying.
1. write the plot to the harddisk:
---
png(toto.png)
plot(c(1,2,3))
dev.off()
Nikol Simecek wrote:
Hello
I am am new to R and any help with the following would be appreciated:
I have a list (example attached) and I would like to create a new list
which is a filtered version of this list. I.e I would like a list that
only contains elements with this value:
Chr 10 :
Petr PIKAL wrote:
Hi
r-help-boun...@r-project.org napsal dne 04.03.2009 13:18:30:
Hello
I am am new to R and any help with the following would be appreciated:
I have a list (example attached) and I would like to create a new list
which is a filtered version of this list. I.e I would
megh megh74 at yahoo.com writes:
To calculate Percentile for a set of observations Excel has percentile()
function. R function quantile() does the same thing. Is there any
significant difference btw percentile and quantile?
If you check the documentation of quantile, you will note that
download the source tarball - everything should be in there.
On Wed, Mar 4, 2009 at 5:15 AM, mau...@alice.it wrote:
Thank you.
I do not know where the source code is available.I never uploaded anything to
CRAN. I am just an R user so far.
Can you please tell me ?
I did notify the package
Frank E Harrell Jr f.harrell at vanderbilt.edu writes:
For Linux I recommend the engauge-digitizer package.
It's Java and works under Windows.
Dieter
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hi,I would like to fill the dots in this graph. I would appreciate help to
do that if possible. If fill is not possible, can I make it bright? I do NOT
want to increase the size.
__
dotplot(hu.event~fitted.adj.cat,data=risk.benefit.cast,groups=treat,
auto.key=list(space=right))
On Wed, Mar 4, 2009 at 12:09 AM, Vedula, Satyanarayana
sved...@jhsph.edu wrote:
Hi,
Could someone help with coding this in R?
I need to select one row per patient i in clinic j. The data is organized
similar to that shown below.
Two columns - patient i in column j identify each unique
Dear R-Experts,
I have an empirical dataset with 150 subjects for 24 observations.
In each observation, each subject can have a score in the range 0:3.
I made then a simple index making the sum of the values in each row,
so each subject have a score between 0 and 72.
I was thinking about what
Hi,
I think you need to set the parameter pch to 19.
- pch=19: solid circle,
- pch=20: bullet (smaller circle),
- pch=21: circle,
- pch=22: square,
- pch=23: diamond,
- pch=24: triangle point-up,
- pch=25: triangle point down.
Regards
2009/3/4 Veerappa Chetty
Dear all,
I have a problem with the integration of the following function.
Could you please give some suggestions?
Thank you very much!!!
y-rnorm(n=100)
f-function(x,xi,h){
n-length(xi)
Ke-c()
for(t in 1:n) {
Ke[t]-dnorm((x-xi[t])/h)
}
fke-sum(Ke)*(1/(n*h))
fke-fke^2
return(fke)
}
Many thanks for all your suggestions - much appreciated!
Nikol
__
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,
Hello!
I would like to ask whether the seasonality implemented in arima() is additive
or multiplicative? I searched a lot, but I could not find an answer to that
question, although it has been asked other times too.
Thank you very much for your attention.
Regards,
Martin
Dear All,
1. Does anyone have experience of running R on a server inside a Citrix
shell - I'd like to get R onto the server and would be greatful for any
tips or direction on the matter.
2. This may seem like a silly question so forgive my ignornace.
Most of the data I currently work with is
andrea.tor...@apat.it wrote:
Dear all,
I have a problem with the integration of the following function.
Could you please give some suggestions?
Thank you very much!!!
y-rnorm(n=100)
f-function(x,xi,h){
n-length(xi)
Ke-c()
for(t in 1:n) {
Ke[t]-dnorm((x-xi[t])/h)
}
To see the code (in this case at any rate) you just type the name of
the function. I produced a snippet from what appeared on my screen. I
then copy-pasted to a text-editor, copy-pasted the missing section to
the d8 segment and assigned the to to wt.filter, which effectively
overwrote the
Rolf == Rolf Turner r.tur...@auckland.ac.nz writes:
On 4/03/2009, at 11:50 AM, Michael A. Miller wrote:
Sports scores are not statistics, they are measurements
(counts) of the number of times each team scores. There
is no sampling and vanishingly small possibility of
Dear Ajay,
just to deny the implicit statement 'corporate user'='moron' surfacing
here and there in this interesting thread :^). This might be a
statistical regularity but should by no means be considered a theorem,
as there are counter-examples available. You can find people willing to
learn
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: Wednesday, March 04, 2009 3:17 AM
To: Daniel Nordlund
Cc: r-help@r-project.org
Subject: Re: [R] How to generate fake population (ie. not
sample) data?
On Wed, Mar 4, 2009 at 2:48 AM, Daniel
Yes, I aware of those definitions. However I wanted to know the difference
btw the words Percentile and quantile, if any. Secondly your link
navigates to some non-english site, which I could not understand.
Dieter Menne wrote:
megh megh74 at yahoo.com writes:
To calculate
#Hi,
#
#somebody knows how to remove the outer box around a wireframe and
reduce the height
#
#
test = data.frame(expand.grid(c(1:10), c(1:10)))
z = test[,1] + test[,2]
test = cbind(test, z)
names(test) = c(x, y, z)
require(lattice)
wireframe(z ~ x*y, data = test, par.box = c(col =
On 04-Mar-09 16:10:29, megh wrote:
Yes, I aware of those definitions. However I wanted to know the
difference btw the words Percentile and quantile, if any.
Secondly your link navigates to some non-english site, which I could
not understand.
Percentile and quantile are in effect the same
Hi guys I have been using R for a few months now and have come across an
error that I have been trying to fix for a week or so now.I am trying to
build a classifer that will classify the wine dataset using Naive Bayes.
My code is as follows
library (e1071)
wine-
Hello - and sorry for what might look like a simple graphics question.
I am building an interaction plot for d:
d=data.frame(xx=c(3,3,2,2,1,1),yy=c(4,3,4,3,4,3),zz=c(5.1,4.4,3.5,3.3,-1.1,-1.3))
d[[1]]-as.factor(d[[1]])
d[[2]]-as.factor(d[[2]])
print(d)
interaction.plot(d$xx, d$yy, d$zz,
(Ted Harding) wrote:
snip
So, with reference to your original question
Excel has percentile() function. R function quantile() does the
same thing. Is there any significant difference btw percentile
and quantile?
the answer is that they in effect give the same results, though
differ
Hi All,
This is a slightly arcane question, but I'm wondering if anyone else
uses vi mode with R? On my platform, across several versions, there is
some broken behaviour. When executing commands like 'df)' (to delete
up to the next bracket) the cursor moves to the next ), but nothing is
(Sorry for the repeat. Forgot to copy R-help)
Try,
test = data.frame(expand.grid(c(1:10), c(1:10)))
z = test[,1] + test[,2]
test = cbind(test, z)
names(test) = c(x, y, z)
require(lattice)
wireframe(z ~ x*y, data = test,
par.settings = list(axis.line = list(col = transparent)),
par.box = c(col
Thank you for finding this. Yes in some cases the parameter settings need to
be updated by a call to plot.new for the calculations to be correct (if you
carried out your example 2 more times you would see that the 3rd plot is also
incorrect since it is still using the dimensions of the 2nd
On 04-Mar-09 16:56:14, Wacek Kusnierczyk wrote:
(Ted Harding) wrote:
snip
So, with reference to your original question
Excel has percentile() function. R function quantile() does the
same thing. Is there any significant difference btw percentile
and quantile?
the answer is that they in
Here is another approach that still uses strspit if you want to stay with that:
tmp - '(-0.791,-0.263].(-38,-1.24].(0.96,2.43]'
strsplit(tmp, '\\.(?=\\()', perl=TRUE)
[[1]]
[1] (-0.791,-0.263] (-38,-1.24] (0.96,2.43]
This uses the Perl 'look-ahead' indicator to say only match on a
:-) works!
Sundar Dorai-Raj schrieb:
(Sorry for the repeat. Forgot to copy R-help)
Try,
test = data.frame(expand.grid(c(1:10), c(1:10)))
z = test[,1] + test[,2]
test = cbind(test, z)
names(test) = c(x, y, z)
require(lattice)
wireframe(z ~ x*y, data = test,
par.settings = list(axis.line =
Hello,
I have two time series objects, 1 is yearly (population) and the other is
quarterly (bankruptcy statistics). I would like to produce a quarterly time
series object that consists of bankruptcy/population. Is there a pre-built
function to intelligently divide these time series.
The series I
Greg Snow wrote:
Here is another approach that still uses strspit if you want to stay with
that:
tmp - '(-0.791,-0.263].(-38,-1.24].(0.96,2.43]'
strsplit(tmp, '\\.(?=\\()', perl=TRUE)
[[1]]
[1] (-0.791,-0.263] (-38,-1.24] (0.96,2.43]
This uses the Perl 'look-ahead'
For the record
residuals(model)
1 2 3 4 5
5.55860143 -0.00073852 2.49255235 -1.41987341 -0.00042425
6 7 8
-0.94389158 2.72987046 -1.15760836
residuals(model, pearson)
1 2 3
Hi
I'm trying to do a GAM analysis and have the following codes entered
into R (density is = sample number, alive are the successes)
density-as.real(density)
y-cbind(alive,density-alive)
library(mgcv)
m1-gam(y~s(density),binomial)
at which point I get the following error message
Error in
Hi Jim,
How are you? I saw your posting. I am trying to do clustering for co
authorship.What I have is undirected graph .I want to have clusters for 393
nodes.
I am attaching the file along with this mail.If you move to the section
Cluster I am looking to do something like that.Is it
Dear R Community,
I am plotting this simple x-y plot (raw data plot attached).
I cant fit a linear regression line to it. I have to figure out what is the
best fit for this graph. Is there a way to tell which regression to use for
this kind of data?
Also, after selecting the best fit model, I
Dear R-experts,
recently, I started to discover the world of R. I came across a problem,
that I was unable to solve by myself (including searches in R-help, etc.)
I have a flat table similar to
key1key2value1
abcd_1 BP 10
abcd_1 BSMP1A
abcd_1 PD 25
abcd_2 BP 20
Hi,
I am looking for a colormap (in color) that look like a gradient in gray
scale. It is to allow people without color printer to print the color graph
and have something meaningful in gray scale.
It can be something like this
plot(1:6,col=c(1,7,5,3,2,4),pch=c(1,20,20,20,20,20))
but with
On Wed, Mar 4, 2009 at 10:52 AM, Dimitri Liakhovitski ld7...@gmail.com wrote:
Hello - and sorry for what might look like a simple graphics question.
I am building an interaction plot for d:
d=data.frame(xx=c(3,3,2,2,1,1),yy=c(4,3,4,3,4,3),zz=c(5.1,4.4,3.5,3.3,-1.1,-1.3))
See ?reshape
Uwe Ligges
Christian Pilger wrote:
Dear R-experts,
recently, I started to discover the world of R. I came across a problem,
that I was unable to solve by myself (including searches in R-help, etc.)
I have a flat table similar to
key1key2value1
abcd_1 BP 10
abcd_1
At 12:19 AM 3/4/2009, you wrote:
plot(x,rho,pch=id)
Or this.
Tim
dat
id x rho
1 A 1 0.1
2 B 20 0.5
3 C 2 0.9
labels-dat$id
labels
[1] A B C
plot(dat$x,dat$rho,pch=labels)
__
R-help@r-project.org mailing list
On Wed, 4 Mar 2009, thibert wrote:
Hi,
I am looking for a colormap (in color) that look like a gradient in gray
scale. It is to allow people without color printer to print the color graph
and have something meaningful in gray scale.
It can be something like this
anujgoel wrote:
Dear R Community,
I am plotting this simple x-y plot (raw data plot attached).
I cant fit a linear regression line to it. I have to figure out what is the
best fit for this graph. Is there a way to tell which regression to use for
this kind of data?
Also, after selecting the
Hello,
I have two time series objects, 1 is yearly (population) and the other is
quarterly (bankruptcy statistics). I would like to produce a quarterly time
series object that consists of bankruptcy/population. Is there a pre-built
function to intelligently divide these time series:
br.ts =
Hello!
I asked in this forum about what kind of seasonality the function arima() from
stats implements. Now that I have been answered that it implements the
Box-Jenkins multiplicative seasonality, I would like to ask whether there is in
R possibility to model ARIMA with additive seasonality. I
On Mar 4, 2009, at 1:22 PM, anujgoel wrote:
Dear R Community,
I am plotting this simple x-y plot (raw data plot attached).
I cant fit a linear regression line to it. I have to figure out what
is the
best fit for this graph.
That is virtually impossible to define rigorously. The best fit
Pls forgive me heavy-handed data generation -- newby ;-)
### start ###
# example data
g - rep.int(c(A, B, C, D), 125)
t - rnorm(5000)
a - sample(t, 500, replace=TRUE)
b - sample(t, 500, replace=TRUE)
# what I actually want to have:
boxplot(a | b ~ g)
# but that does obviously not produce what
Excel 2003's help for percentile just says it interpolates
between the quantiles in the data:
Array is the array or range of data that defines relative standing.
K is the percentile value in the range 0..1, inclusive.
If array is empty or contains more than 8,191 data points,
1 - 100 of 165 matches
Mail list logo
