Here is one way of doing it:
x=c(1,2,3,1)
y=c(1,2,3,1)
ss=c(55,NA,55,88)
ss_byxy_test=tapply( ss, list( x, y), mean, na.rm=TRUE)
# use the 'reshape' package
ss_byxy_test
# now 'melt' the data to get it into a format for plotting
(ss_melt - melt(ss_byxy_test))
# create the plot area so you can
I left off the statement to load the reshape package. If you don't have it,
install it from CRAN:
x=c(1,2,3,1)
y=c(1,2,3,1)
ss=c(55,NA,55,88)
ss_byxy_test=tapply( ss, list( x, y), mean, na.rm=TRUE)
ss_byxy_test
# use the 'reshape' package
library(reshape)
# now 'melt' the data to get it into a
On 12/05/2009 11:45 PM, Tal Galili wrote:
Hi Romain,
First - I'd like to say that you where one of the bloggers I wanted to
e-mail and didn't get to it. I still remember our short session in useR
2009 where you showed me your plans for an R editor - I am still waiting
to see it once it comes
hi there!
how u doin?...hope fine.
look, i've a difficulty in simulation especially for the bayesian analysis
using R. Using the famous Markov Chain Monte Carlo (MCMC) method,
Methropolis hasting algorithm i want to simuate different distrbution like:
1/ Weibull-exponential
2/Gamma-Gamma
Hi, dear R users
I am a newbie in R and I wantto use the method of meximum likelihood
to fit a Weibull distribution to my survival data. I use optim as
follows:
optim(c(1, 0.25),weibull.like,mydata=mydata,method=L-BFGS-B,hessian
= TRUE)
My question is: how do I setup the constraints so that
Hi, dear R users
I am a newbie in R and I need to use the method of meximum likelihood to fit a
Weibull distribution to my survival data. I use optim as follows:
optim(c(1.14,0.25),weibull.like,mydata=mydata,method=L-BFGS-B,hessian = TRUE)
My question is: how do I setup the constraints that
Hello list,
Thanks in advance for whatever help you can give me. I have a data set that I
want to visualize that has four important dimensions (x,y,z,t).
What I would like to be able to do is plot the three spatial coordinates using,
say, scatterplot3d (or similar), and have the color of the
3rd and 4th dimension can easily be accommodated in a 2d plot by using
different color, shape of points etc. Therefore you can start with ggplot
package. You might also want to look into here :
http://had.co.nz/ggplot2/geom_point.html
Best,
Jared Nance wrote:
Hello list,
Thanks in advance
As per ?optim
Usage:
optim(par, fn, gr = NULL, ...,
method = c(Nelder-Mead, BFGS, CG, L-BFGS-B, SANN),
lower = -Inf, upper = Inf,
control = list(), hessian = FALSE)
Note that the optimx() function on R-forge
http://r-forge.r-project.org/R/?group_id=395
Thanks Dennis and Jorge!
On Sat, Dec 5, 2009 at 8:24 PM, Jorge Ivan Velez
jorgeivanve...@gmail.comwrote:
Dear Gary,
Here is a suggestion using read.table() and textConnection():
toread - id sex age inc r1 r2 r3
1 F 35 17 7 2 2
17 M 50 14 5 5 3
33 F 45 6 7 2 7
49 M 24 14 7
Sharpie wrote:
Tariq Perwez wrote:
Hi Everyone,
I have two very basic questions and would appreciate your help.
1. I would like to see/access the data that comes with a given R package
by
using a function like read.table(). For example, I just installed car
package from CRAN. I know that
Sarah wrote:
I am trying to fit a linear model with seasonal effects but keep getting
'Incorrect number of dimensions'. I have no idea what this means or how to
fix it. I am following instructions received from an instructor and it
worked for him, so I assume it has something to do with my
Please read ?optim and about its arguments
lower, upperBounds on the variables for the L-BFGS-B method.
Uwe Ligges
Liang Wang wrote:
Hi, dear R users
I am a newbie in R and I need to use the method of meximum likelihood to fit a Weibull
distribution to my survival data. I use optim as
hansoffate wrote:
I'm running an LSODA to generate some graphs, but I need to stop at a certain
point and use those values to generate another LSODA output. This is
working fine, but when I try to run the second LSODA, I get the Error in
eval(expr, envir, enclos) : object 'N' not found. Any
Peng Cai wrote:
Hi,
I'm plotting grouped barplot using following code and data. I need help with
re-ordering the labels.
1. On x-axis the factor AGE is grouped in order 0--4, 15--18, 5--14;
whereas I would like to have it in 0--4, 5--14, 15--18.
2. If I need to re-order RACE variable. How
Hi Romain.
1) Thank you for the link to how to add myself to PlanetR.
2) What I meant was that you don't need an online aggregator, since you know
how to use it with your own local software.
Cheers,
Tal
Contact
Details:---
Read the help page for `optim'. You already seem to be aware that L-BFGS-B
should be used. How much more harder is to read the help page further to
figure out how to supply the constraints?
There are also various other algorithms, but L-BFGS-B should be able to do
the job.
Ravi.
Tal,
For completeness:
On 5 December 2009 at 22:32, Tal Galili wrote:
| 1) Planet R is limited (for years) to 26 feeds only, and I don't remember
| seeing it evolve to include (or allow inclusion) of new R blogs that came
| around.
That was already addressed by Romain. It would appear that you
Dear list,
as described beforehand, i wanted to merge three single plots into one
JPEG-file. Since I use PNG now, code goes as follows:
png(...)
layout(...)
par(...)
plot(...)
par(...)
plot(...)
par(...)
plot(...)
which works out pretty well, but now a new problem occurred: R has
stopped
You need dev.off() I think.
-Ista
On Sun, Dec 6, 2009 at 10:53 AM, Walther, Alexander
awaltherm...@googlemail.com wrote:
Dear list,
as described beforehand, i wanted to merge three single plots into one
JPEG-file. Since I use PNG now, code goes as follows:
png(...)
layout(...)
Thanks Uwe, I got your suggestions part too.
2009/12/6 Uwe Ligges lig...@statistik.tu-dortmund.de
Peng Cai wrote:
Hi,
I'm plotting grouped barplot using following code and data. I need help
with
re-ordering the labels.
1. On x-axis the factor AGE is grouped in order 0--4, 15--18,
Thank you Uwe, Dennis, and Gary for your help. I have one more question:
I'm using pre-defined y-scales and trying to create grid lines.
As Female category has low sum value, its y-axis range from 0-150 whereas
Male ranges from 0-300. Is it possible to make them on same scale. Here is
the
@ Ista: no, R still struggles with the output. Thanks anyway. :) Is
there a way to merge multiple of these PNGs into a single PDF document
afterwards? I would like to refrain from bridging this gap with
additional software like f.i. Adobe Acrobat.
Cheers
Alex
Ista Zahn schrieb:
You need
Please ignore my last question. I found a way to handle that. One last
thing:
I'm defining my own y scales. In the process the bar starts from below the
y=0 line (or below the y-axis). Is there a way to get rid of it.
Here is the code and data I'm using
R Code (Data read in object dta)
dta$age
Although this is really a bibtex/Jabref question, I am hoping this
might still be the best place for help.
The bibtex entry in Jabref is like this
@MANUAL{R2009,
title = {R: A Language and Environment for Statistical Computing},
author = {{R Development Core Team}},
organization = {R
On Dec 6, 2009, at 11:47 AM, Walther, Alexander wrote:
@ Ista: no, R still struggles with the output. Thanks anyway. :) Is
there a way to merge multiple of these PNGs into a single PDF document
afterwards? I would like to refrain from bridging this gap with
additional software like f.i. Adobe
Oh yes, that's true, but I don't understand why the same exact model works
fine the first time I run LSODA, then running the model again with different
PARMS causes it to error with this error. I've been trying to debug or
find other relevant posts with this same error, but so far I have had no
Hi All,
I'm fitting an intercept ANOVA model using R. I'm using following contrast
statement, so my first level is considered as base level.
options(contrasts = c(contr.treatment, contr.treatment))
My data has two factors are Month and WeekDay. Due to default alphabetic
ordering, the
Hello there,
How would I go about creating a function that will read a vector containing
names of variables and then will output tables of proportions for all the
categorical variables and summary statistics appropriate for all quantitative
variables for all variables read by the function.
In order to get user input when running a batch file with Rterm in Windows
XP and Windows Server 2003, I have been using the R script below with
R-2.8.1 and earlier. This opens a command line window with the statement
Please enter the ID: and feeds the users response into the variable 'id'.
I have a one-variable data set in R.
The plot of histogram of my numerical variable suggests an inverse
gaussian distribution.
How can I obtain best estimation for the two parameters of inverse
gaussian based on my data?
Thanks.
--
View this message in context:
On Dec 6, 2009, at 7:19 AM, Kristen Mark wrote:
Hello there,
How would I go about creating a function that will read a vector
containing names of variables and then will output tables of
proportions for all the categorical variables and summary statistics
appropriate for all
Hi guys
I have data that contains the variances vt of the yields of 1, 2, 3, 4,
5,10, 20 year bonds. Assuming the Hull-White model for the yield of a t-year
zero-coupon bond, I have to estimate the σ of the Hull-White model using
nonlinear least squares and give a 95% confidence interval for each
hansoffate wrote:
Oh yes, that's true, but I don't understand why the same exact model works
fine the first time I run LSODA, then running the model again with different
PARMS causes it to error with this error. I've been trying to debug or
find other relevant posts with this same error, but
On Dec 6, 2009, at 1:14 PM, Xin Ge wrote:
Hi All,
I'm fitting an intercept ANOVA model using R. I'm using following
contrast
statement, so my first level is considered as base level.
options(contrasts = c(contr.treatment, contr.treatment))
My data has two factors are Month and WeekDay.
On Dec 6, 2009, at 12:48 PM, noclue_ wrote:
I have a one-variable data set in R.
The plot of histogram of my numerical variable suggests an inverse
gaussian distribution.
How can I obtain best estimation for the two parameters of inverse
gaussian based on my data?
Hi Tim,
Take a look at
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html
HTH,
Jorge'
On Sun, Dec 6, 2009 at 12:48 PM, noclue_ wrote:
I have a one-variable data set in R.
The plot of histogram of my numerical variable suggests an inverse
gaussian distribution.
Hello,
I want to write a function calculating posterior probability.
Can anyone help have a look at the following program? Here function f is
what I eventually
need. It is related to function f1, f2 and f3. The argument x is a vector.
I get evaluation errors here.
Thank you very much!
f1
On Dec 6, 2009, at 1:08 PM, king123 wrote:
Hi guys
I have data that contains the variances vt of the yields of 1, 2, 3,
4,
5,10, 20 year bonds. Assuming the Hull-White model for the yield of
a t-year
zero-coupon bond, I have to estimate the σ of the Hull-White model
using
nonlinear
I used to have the following in my .Rprofile:
if (length(.libPaths())==1)
.libPaths(paste(Sys.getenv(HOME),/Library/R/,paste(R.version$major,as.integer(R.version$minor),sep='.'),/library,sep=''))
-- and it added my user-defined library directory. Then I installed
packages there, so during an
On Sat, Dec 05, 2009 at 01:39:03PM -0500, David Winsemius wrote:
The question lacks so many specifics that would be needed to address
problem that I suspect persons with sufficient knowledge to address it
(and I am probably not among that group) simply passed it over and
waited for you
I have the .libPaths() defined in my .Rprofile, and it stopped having any
effect. Even if I try to do it in the REPL, as:
x -
paste(Sys.getenv(HOME),/Library/R/,paste(R.version$major,as.integer(R.version$minor),sep='.'),/library,sep='')
x
[1] /Users/alexyk/Library/R/2.10/library
On 06/12/2009 2:53 PM, Alexy Khrabrov wrote:
I have the .libPaths() defined in my .Rprofile, and it stopped having any
effect. Even if I try to do it in the REPL, as:
x -
paste(Sys.getenv(HOME),/Library/R/,paste(R.version$major,as.integer(R.version$minor),sep='.'),/library,sep='')
x
[1]
Hi,
I'm trying to create a new variable using the following code:
A$B2-ifelse(A$B %in% EDUCATION,EDUCATION,NA)
The problem is, I want the new variable to be labeled EDUCATION if the
original *contains* EDUCATION, not of it is an exact match. I've done
a lot of searching on wild cards and
On 06/12/2009 3:05 PM, Jennifer Brea wrote:
Hi,
I'm trying to create a new variable using the following code:
A$B2-ifelse(A$B %in% EDUCATION,EDUCATION,NA)
The problem is, I want the new variable to be labeled EDUCATION if the
original *contains* EDUCATION, not of it is an exact match. I've
On Sun, Dec 6, 2009 at 3:04 PM, Duncan Murdoch murd...@stats.uwo.ca wrote:
Does that directory exist? .libPaths(foobar) silently does nothing,
because (on my system) foobar is not a directory. (This is documented
behaviour, though it would perhaps be friendlier if it gave a warning when
it
On 7/12/2009, at 8:32 AM, Alexy Khrabrov wrote:
I used to have the following in my .Rprofile:
if (length(.libPaths())==1)
.libPaths(paste(Sys.getenv(HOME),/Library/R/,paste(R.version
$major,as.integer(R.version$minor),sep='.'),/library,sep=''))
-- and it added my user-defined library
Hi
Sebastien Bihorel wrote:
Dear R-users,
I would like to know how to pass arguments to gpar() without hard-coding
them. I tried to store my arguments in a list and passed this list to
gpar(), but it did find the way to do it properly. Any help would be
appreciated.
a-
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of William Dunlap
Sent: Friday, December 04, 2009 3:29 PM
To: Dennis Fisher; r-h...@stat.math.ethz.ch
Subject: Re: [R] Use of apply rather than a loop
You could try using merge:
On 4/12/2009, at 10:10 PM, Gough Lauren wrote:
Hi,
I have a numeric pixel image which I would like to divide into factors
for analysis in Spatstat. I have found that I can use cut.im()
function
to divide the range of pixel values into a series of equal length
intervals (e.g. if my
Hi,
I have 49 pairs in my data.frame 'data'
x y
76 80
138 143
67 67
29 50
381 464
23 48
37 63
120 115
... ...
how do I get a bootstrap sample of size n=50?
i have tried
sample(data,size=50,replace=TRUE)
and
sample(data,replace=TRUE)
both didnt
Jorge Ivan Velez wrote:
Hi Tim,
Take a look at
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html
HTH,
Jorge'
i don't think this is what the OP wanted. The nig is the normal inverse
Gaussian distribution, which is a special case of the generalized
inverse
Gary,
if i were you, i would use scan().
here is a piece of code.
# DO DATA INPUT IN R CONSOLE WITH SCAN() #
#--#
# COMPARABLE SAS CODE: #
#--#
# data
casperyc casperyc at hotmail.co.uk writes:
I have 49 pairs in my data.frame 'data'
x y
7680
... ...
how do I get a bootstrap sample of size n=50?
data[sample(nrow(data),size=50,replace=TRUE),]
__
R-help@r-project.org mailing list
Hi there,
i think that's not what i was aiming for...
i was aked to
Generate 50 Bootstrap samples and corresponding estimates
if i do data[sample(nrow(data),size=50,replace=TRUE),]
it will give me a table of 50 rows ( 50 sets of x and y)
then how do i estimate the mean? ( mean was esitmated
Hi,
I have a glm gives summary as follows,
Estimate Std. Errorz valuePr(|z|)
(Intercept) -2.03693352 1.449574526 -1.405194 0.159963578
A0.01093048 0.006446256 1.695633 0.089955471
N0.41060119 0.224860819 1.826024
data-data.frame(x=rnorm(49), y=rnorm(49))
t(sapply(1:50,function(i){colMeans(data[sample(nrow(data),size=nrow(data),replace=TRUE),])}))
r-help@r-project.org
2009/12/7 casperyc caspe...@hotmail.co.uk:
Hi there,
i think that's not what i was aiming for...
i was aked to
Generate 50
Hi everyone,
I have successfully created and analyzed my network data. I am new to R, and
Network Analysis too, but I want to color my vertex based on some of the
centrality measures calculated.
Can someone point me in the right direction? Again, I am new to R, but given
how powerful R
David Scott wrote:
Jorge Ivan Velez wrote:
Hi Tim,
Take a look at
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html
HTH,
Jorge'
i don't think this is what the OP wanted. The nig is the normal inverse
Gaussian distribution, which is a special case of the
Hello Everyone.
I¹m on a Mac, using rgl. Thanks to all the developers and maintainers on
it!
If I drag the mouse, I can rotate the view, if I hold down the option key or
the ctrl key and drag, I get scaling (though at times I seem to alter the
perspective?). Adding the shift key doesn't seem
Hi, Prof Nash
Thanks for your comment!
I modified my code to be (added an extra parametr):
optim(c(1.14,0.25,0.06), weibull.like, mydata=mydata, method=L-BFGS-
B, hessian = TRUE, lower = c(0, 0, 0), upper = c(Inf, Inf, 1))
But I had the following error:
Error in optim(c(1.14, 0.25, 0.06),
A new version of the ggplot2 web application is available at
http://www.yeroon.net/ggplot2. New features include 1D geom’s
(histogram, density, freqpoly), syntax mode (by clicking the tiny
arrow at the bottom), and some additional facet options. Furthermore
some minor improvements and fixes, most
Hello,
Is there a function that calculates MFE(Mean Forecast Error) and MSFE(Mean
Square Forecast Error) in R?
Cant find anything!
Thanks Tobias
--
View this message in context:
http://n4.nabble.com/Function-for-MFE-and-MSFE-tp950368p950368.html
Sent from the R help mailing list archive at
Hi，
I have a dataframe with x and y in pairs, ( 49 pairs i.e. 49 rows )
x y
76 80
138 143
67 67
29 50
381 464
23 48
37 63
120 115
how do i generate 50 Bootstrap samples and corresponding estimates?
i have tried
sample( data, size=50,replace=T)
does not seem to be right
Thanks.
How would I verify that one variable is the summation of two other variables?
Thank you!
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
I have a huge college dataset and I want to be able to find which colleges
appear more than once in the dataset. How would I go about doing this? The
variable college is a factor. If I just do table(college) it is impossible to
find where there are repetitions. Any help would be greatly
noclue_ wrote:
I have a one-variable data set in R.
The plot of histogram of my numerical variable suggests an inverse
gaussian distribution.
How can I obtain best estimation for the two parameters of inverse
gaussian based on my data?
Thanks.
See the examples for function pelp in
Hi,
How do I generate bootstrap p-values for goodness of fit test.
Thank you in advance for your help,
Juliet
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
college=c(Columbia,Columbia,Harvard,Harvard,NYU)
unique(college)
table(college)
unique(college)[table(college)1]
HTH
Daniel
-
cuncta stricte discussurus
-
-Ursprüngliche Nachricht-
Von: r-help-boun...@r-project.org
Tena koe Kristen
all.equal(theOneVariable, theSum)
HTH
Peter Alspach
PS It is advisable to use a more informative subject line.
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Kristen Mark
Sent: Monday, 7 December 2009
Hi there,
This seems to 'simple' for that. ( I mean the codes are too good)
Can I ask for a bit explaination? or some simplier (maybe few more line
codes) to archieve that goal?
And the codes you gave here,
first line give a table with x and y , random normal values,
but the second line, what
Here you have an example:
colleges - factor( c(a, a, b, c, c ))
colleges
[1] a a b c c
Levels: a b c
table( colleges )
colleges
a b c
2 1 2
colleges.count - table( colleges )[colleges]
colleges.count
colleges
a a b c c
2 2 1 2 2
which( colleges.count 1 )
a a c c
1 2 4 5
Kristen Mark
Bryan Hanson wrote:
Hello Everyone.
I¹m on a Mac, using rgl. Thanks to all the developers and maintainers on
it!
If I drag the mouse, I can rotate the view, if I hold down the option key or
the ctrl key and drag, I get scaling (though at times I seem to alter the
perspective?). Adding the
Thanks Duncan... I see how to do it. It appears that all 3 of my mouse
buttons are actually set to the trackpad button on my MacBook Pro. So if I
assign any button via pan3d I get the same behavior, and can't get back to
the default. I'll play with it!
Bryan
On 12/6/09 9:14 PM, Duncan
I am using a plot to inspect data points, and I would like to identify each
point with respect to an ID. At issue is that I am producing a faceted plot
with many IDs (96) and the key is far to large to accurately identify points
by color.
1) Can you direct me on labeling or printing data points
Hi Jared,
google R graph library 4d or
http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=90
bests
milton
On Sun, Dec 6, 2009 at 7:31 AM, Jared Nance nanc...@phys.washington.eduwrote:
Hello list,
Thanks in advance for whatever help you can give me. I have a data set
that I want
Hi Dave,
May be a good idea you (1) use subjects avoiding help, Hi there,
problem, etc. The list is archived for future search, and nobody with
problems with MCMC, Metropholis Hasting algorithm will search Hi there;
(2) how about minimum reproducible problem? If you are needing use this kind
of
Sorry. My previous code didn't quite display my problem correctly. To make
my point more clear, I want to use my groups as the label for the point,
or rather the point itself. The key becomes illegible with many groups.
x-as.factor(1:5)
y-rnorm(1:10)
f1-c(a,b)
f2-c(x,y)
dat-expand.grid(x=x,
Hi Tom,
The problem was that the bars have width 22.5, so setting the limits
to c(0, 360) chopped off the last bar. To get around that, you needed
to add half a bar width (11.25) to the start and end limits, and then
rotate the polar coordinate system so that 0 lines up again.
ggplot(wind.data,
I've just developed edtdbg, a small package that integrates R's debug()
with one's text editor.
Excerpt from the README file:
Goal
The debug() function in R is primitive. My goal was to make it more
usable by integrating it with one's text editor. Hence I wrote the
package
Hello
On 12/7/09, Norm Matloff matl...@cs.ucdavis.edu wrote:
Please give me your comments. I probably will upload to CRAN at some
point, possibly after an ESS guru contributes the ESS code. :-)
Do you think that the package can be made to work with Geany?
Thank you
Liviu
81 matches
Mail list logo
