That is the thing. As a new comer to 'R' I don't understand how to write a
formula when all I have is a time series. I don't know how to express the
independent and dependent variables in a formula when the object is a time
series. So please just solve this simple example and I will extrapolate
See ?comment
On Fri, 5 Sep 2008, Alexy Khrabrov wrote:
Is there a way to associate descriptions with the objects in the workspace,
and later retrieve them to know what the object was created for?
Thanks,
Alexy
__
R-help@r-project.org mailing list
On Fri, 5 Sep 2008, [EMAIL PROTECTED] wrote:
That is the thing. As a new comer to 'R' I don't understand how to write
a formula when all I have is a time series. I don't know how to express
the independent and dependent variables in a formula when the object is
a time series. So please just
hi
I am expecting the mails that will flow in R- help after subscribing this but
could not. Can you provide me the right form.
Thanks.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Dear R community
Hope every one be in best of his/her health. I have a situation in which there
are s-sectors. Each sector is further divided into r-rows and c-columns. All it
makes an array having dimension (r,c,s). e.g.
Hi im trying to add error bars to my barplots, there very basic, i have a few
grapghs where the y variable is different but on all the X variable is Age
(Adult and Juvenile) however this is split into two levels so i have males and
females, so my graph basically has four bars on it.
I know
Just another remakr on this thread.
I you have time series and think about its fourier transform
(EE language) then you should know that the statistical language
of that is spectral analysis or maybe
frequency domain time-series analysis
and the R function to consider should definitely be
New major versions of the caret packages (caret 3.37, caretLSF 1.23 and
caretNWS 0.23) have been uploaded to CRAN.
caret is a package for building and evaluating a wide variety of predictive
models. There are functions for pre-processing, tuning models using
resampling, visualizing the results,
hi im starting with R.have no idea to start...plz help
--
View this message in context:
http://www.nabble.com/new-to-R-tp19342903p19342903.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
On Sat, Sep 6, 2008 at 5:10 AM, sudeshna [EMAIL PROTECTED] wrote:
hi im starting with R.have no idea to start...plz help
Search the Internet for online tutorials and/or read an introductory
book (search for them, e.g., on Amazon.Com).
Good luck,
Paul
sudeshna wrote:
hi im starting with R.have no idea to start...plz help
Hi sudeshna,
There are several beginner's guides on the CRAN website. Go to:
http://cran.r-project.org
and select Contributed (second last option on the left).
Jim
__
Brown, Heidi wrote:
Hi im trying to add error bars to my barplots, there very basic, i have a few
grapghs where the y variable is different but on all the X variable is Age
(Adult and Juvenile) however this is split into two levels so i have males and
females, so my graph basically has four
Hello Martin,
Zitat von Martin Maechler [EMAIL PROTECTED]:
Just another remakr on this thread.
I you have time series and think about its fourier transform
(EE language) then you should know that the statistical language
of that is spectral analysis or maybe
frequency domain time-series
On Thu, 4 Sep 2008, Grant Gillis wrote:
Hello Professor Ripely,
Sorry for not being clear. I posted after a long day of struggling. Also
my toy distance matrix should have been symmetrical.
Simply put I have spatially autocorrelated data collected from many points.
I would like to do a
Hi R users,
Is is possible for me to use the try function with boot? I would to do
the bootstraping with a nonlinear model(it works well when R 1000).
But it does not work very well (when R is large) thus I try to use
try to resolve. I put the try function in two cases:
case1: put the
Brown, Heidi wrote:
Hi im trying to add error bars to my barplots, there very basic, i have a few
grapghs where the y variable is different but on all the X variable is Age
(Adult and Juvenile) however this is split into two levels so i have males and
females, so my graph basically has four
Dear R-Users,
I am currently looking for a way to test the equality of two correlations
that are related in a very special way. Let me describe the situation with
an example.
- There are 100 respondents, and there are 2 points in time, t=1 and t=2.
- For each of the respondents and at each of
I have to calculate a formula that gives me a ten components vector. I want
to see how the components behave at varying the variable i.
But when i run the following function:
bbayes-c()
for(i in 1:100) {
+ xx-t(X1_10)%*%X1_10
+ xxmeno1-solve(xx)
+ V-xxmeno1*i
+ Vmeno1-solve(V)
+ tx-t(X1_10)
+
you are subscribed to the r-help list.
On Sat, Sep 6, 2008 at 4:48 AM, Ram Kumar Basnet [EMAIL PROTECTED] wrote:
hi
I am expecting the mails that will flow in R- help after subscribing this but
could not. Can you provide me the right form.
Thanks.
[[alternative HTML version
Hi Ricardo,
You can search for comparisons by entering the packages that interest
you at:
http://finzi.psych.upenn.edu/search.html
Michael Mitchell wrote an interesting comparison of SAS, SPSS, Stata and
R at:
http://www.ats.ucla.edu/stat/technicalreports/
That report says little about R,
Oliver Bandel wrote:
Hello Martin,
Zitat von Martin Maechler [EMAIL PROTECTED]:
Just another remakr on this thread.
I you have time series and think about its fourier transform
(EE language) then you should know that the statistical language
of that is spectral analysis or maybe
ss == stephen sefick [EMAIL PROTECTED]
on Sat, 6 Sep 2008 09:40:50 -0400 writes:
ss you are subscribed to the r-help list.
Huh??? How would *you* know that?
The list of subscribers to R-help is not at all public,
and the fact that Ram Kumar's posting went through may well be
the
FEH == Frank E Harrell [EMAIL PROTECTED]
on Sat, 06 Sep 2008 07:19:33 -0500 writes:
FEH Brown, Heidi wrote:
Hi im trying to add error bars to my barplots, there very basic, i have
a few grapghs where the y variable is different but on all the X variable is
Age (Adult and
I apologize, my naive understanding that if an email was sent to the
list then the sender is subscribed. However, my hypothesis is
falseifiable.
On Sat, Sep 6, 2008 at 10:33 AM, Martin Maechler
[EMAIL PROTECTED] wrote:
ss == stephen sefick [EMAIL PROTECTED]
on Sat, 6 Sep 2008 09:40:50
On 9/6/08, sudeshna [EMAIL PROTECTED] wrote:
hi im starting with R.have no idea to start...plz help
--
http://www.math.ilstu.edu/dhkim/Rstuff/Rtutor.html
http://www.statmethods.net/index.html
http://rforsasandspssusers.com/
Rcmdr
rattle
http://www.sciviews.org/_rgui/
Dear R community
Initially i thought my problem has been solved but one thing which i found e.g.
if
1. All the elements of a sector are zero e.g
, , 7
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]000000000 0
[2,]00000
Hello
How can I change the function to get the rows with the sum (x+y+z) = 10?
Thank you very much
Joseph
- Original Message
From: Marc Schwartz [EMAIL PROTECTED]
To: joseph [EMAIL PROTECTED]
Cc: r-help@r-project.org
Sent: Wednesday, September 3, 2008 3:24:58 PM
Subject: Re: [R]
Hello everyone,
I would like to know if there is any function to calculate the mode value, or I
have to build one to do it.
Thanks so much
Carlos
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
look here:
http://www.nabble.com/how-to-calculate-the-mode-of-a-continuous-variable-td19214243.html#a19214243
On Sat, Sep 6, 2008 at 1:24 PM, Carlos Morales
[EMAIL PROTECTED] wrote:
Hello everyone,
I would like to know if there is any function to calculate the mode value, or
I have to build
Hi Carolos,
I know that it is not a elegant soluction, but may work. Almost for integer
values. Take care with float values.
modevalue-function(x)
{
x.freq-data.frame(table(x))
x.freq.max-max(x.freq$Freq)
x.freq.selected-subset(x.freq, x.freq$Freq==x.freq.max)
Martin Maechler maechler at stat.math.ethz.ch writes:
FEH == Frank E Harrell f.harrell at vanderbilt.edu
on Sat, 06 Sep 2008 07:19:33 -0500 writes:
FEH See http://biostat.mc.vanderbilt.edu/DynamitePlots for many reasons
not
FEH to use dynamite plots.
Ah! Very good!
Hi Jorge
I got the rows where V3 looks like this 10:10:10; Ithe sum here is 30 and not
10.
I want the rows where the sum is 10 for exaple 5:5:0 and 2:2:6
thanks
Joseph
- Original Message
From: Jorge Ivan Velez [EMAIL PROTECTED]
To: joseph [EMAIL PROTECTED]
Sent: Saturday,
On 09/06/08 17:24, Carlos Morales wrote:
Hello everyone,
I would like to know if there is any function to calculate the mode value, or
I have
to build one to do it.
I just did this the other day. Funny you should ask. It finds the
mode of each row of a matrix called Pbest.
Pbest.p -
Hi Jinsong,
I try to put the try in the c1.fun but it still does not work. Do
you mean this?
c1.fun-try(function(data,i){
+ d-data
+ d$density-d$fitted+d$res[i]
+ coef(update(c1.nmf,data=d))
+ } ,silent=T)
c1.try-boot(c1data, statistic = c1.fun, R=5000)
Error in nls(formula = density
#something like this?
V1 - c(1:10)
V2 - c(0,5,7,8,1,6,5,13,7,0)
V3 - c(9,5,6,8,1,7,5,33,88,0)
z - cbind(V1,V2,V3)
row.sums - rowSums(z)
d - cbind(z, row.sums)
subset(d, row.sums==10)
On Sat, Sep 6, 2008 at 2:25 PM, joseph [EMAIL PROTECTED] wrote:
Hi Jorge
I got the rows where V3 looks like this
Sorry I didn't read the problem carefully.
On Sat, Sep 6, 2008 at 2:53 PM, Jorge Ivan Velez
[EMAIL PROTECTED] wrote:
Hi Joseph,
Try this:
# Data set
DF=read.table(textConnection(V1 V2 V3
ab0:1:12
df1:2:1
cd1:0:9
be2:2:6
fc5:5:0),header=TRUE)
Dear R-colleagues,
another question from a newbie: I am creating a lot of simple
pivot-charts from my raw data using the reshape-package. In these charts
we have medical doctors judging videos in the columns and the videos
they judge in the rows. Simple example of chart/data.frame input with
two
Hi,
I just installed R, I'm work in UBUNTU and I don't have idea about how to
run a r-code written in emacs
into the shell.
Well I am in a shell, and obviously I can run simple commands over there,
Must I compile the program? if yes, How must I do that?
what is the extension?
I really
Not sure where your input came from. It's not in a format I would
have expected of an R object and the first line is not in a form that
would be particularly easy to read into a valid R object. Numbers are
no legitimate object names. It's also not clear what you want to do
with the
I fit a random effects linear model to data, and then tried to use it to
predict, but I got this error:
predict(lmeObject, newdata, level=0)
Error in eval(mCall$fixed)[-2] : object is not subsettable
This is a new error for me. It still occurs if I change the level to 1 or if I
change the
Hello Mr. Burns, Hello Mr. Dwinseminus
thank you very much for your incredible quick and efficient reply! I was
completely successful with the following command:
pairs-data.frame(pred=factor(unlist(input[,-c(1,ncol(input))])),ref=factor(input[,ncol(input)]))
In case of the input example
did you install ess for emacs?
On Sat, Sep 6, 2008 at 11:52 AM, Luisa [EMAIL PROTECTED] wrote:
Hi,
I just installed R, I'm work in UBUNTU and I don't have idea about how to
run a r-code written in emacs
into the shell.
Well I am in a shell, and obviously I can run simple commands over
Hi Joseph,
Try this:
# Data set
DF=read.table(textConnection(V1 V2 V3
ab0:1:12
df1:2:1
cd1:0:9
be2:2:6
fc5:5:0),header=TRUE)
closeAllConnections()
target=10
DF[sapply(strsplit(as.character(DF$V3), :), function(x)
sum(as.numeric(x))== target), ]
HTH,
i have a list of 6 each containing a dataframe of 96 observations as a
zoo object. Is there a way to plot these in one frame
par(mfrow=c(3,2))
this is what I tried
lapply(d, FUN=plot)
I can provide data, list is large.
thanks
--
Stephen Sefick
Research Scientist
Southeastern Natural Sciences
On Saturday, 06 September 2008, 08:52 (UTC-0700), Luisa wrote:
Hi,
I just installed R, I'm work in UBUNTU and I don't have idea about how to
run a r-code written in emacs
into the shell.
Well I am in a shell, and obviously I can run simple commands over there,
Must I compile the
Please tell me how to format data in a data frame so when currency amount is
displayed in a chart the axis tick labels contain leading $ signs.
Please also tell me if it is possible to rotate x axis labels using ggplot2.
Thank you,
Kurt
__
Hi Ralph,
I had the same problem you do a few months ago, and realized that
the question I had (does time show a different effect for X than Y) was not
best modeled as differences between correlations across individuals, but as
whether time interacts with condition.
I answered
2008/9/6 Luisa [EMAIL PROTECTED]:
Hi,
I just installed R, I'm work in UBUNTU and I don't have idea about how to
run a r-code written in emacs
into the shell.
Well I am in a shell, and obviously I can run simple commands over there,
Must I compile the program? if yes, How must I do that?
BIG THANKS!!! ))
--
View this message in context:
http://www.nabble.com/Area-of-density-tp19338958p19351299.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
I am trying to define 25 vectors of varying lengths, say y1 to y25 in a loop,
and then store the results of some computations in them. My problem is about
using some sort of concatenation for names. For example, instead of
initializing each of y1 through y25, I would like to do it in a loop.
shalu shahlar at hotmail.com writes:
I am trying to define 25 vectors of varying lengths, say y1 to y25 in a loop,
and then store the results of some computations in them. My problem is about
using some sort of concatenation for names. For example, instead of
initializing each of y1
Dear Friends,
I not sure whether this is an Sweave or a beamer problem.
The Rnw file:
\documentclass[compress,smaller]{beamer}
%\documentclass{article}
%\usepackage{beamerarticle}
\usepackage{Sweave}
\title{Psychophysics II}
\date{September 9, 2008}
\begin{document}
\frame{
\begin{Schunk}
I would suggest that you use a 'list' since it seems that the result
is more than a single value:
bbayes-list()
for(i in 1:100) {
xx-t(X1_10)%*%X1_10
xxmeno1-solve(xx)
V-xxmeno1*i
Vmeno1-solve(V)
tx-t(X1_10)
prpar-solve(Vmeno1+xx)
snpar-tx%*%y
This may come closer since it removes the zeros before comparison:
x=c(1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,2,2,2,3,3,0,0,0,0,0,0,0,0,0,0,1,1,1,2,2,3,3,3,4,4,4,0,0,0,0,0,0,0,1,1,1,2,2,2,3,3,3,4,4,4,
+ 0,0,0,0,0,0,1,2,2,2,2,2,0,3,3,0,4,4,0,0,0,0,0,0)
x=array(x,dim=c(3,6,5))
I would suggest that you use a list to store the values since it is
easier to create and reference:
output - list()
for (i in 1:10) output[[i]] - seq(i)
output
[[1]]
[1] 1
[[2]]
[1] 1 2
[[3]]
[1] 1 2 3
[[4]]
[1] 1 2 3 4
[[5]]
[1] 1 2 3 4 5
[[6]]
[1] 1 2 3 4 5 6
[[7]]
[1] 1 2 3 4 5 6 7
Hi Michael,
I think that beamer needs that frames that contain any verbatim text
or R output to be declared as fragile.
Try
\begin{frame}[fragile]
...
\end{frame}
On Sat, Sep 06, 2008 at 06:22:41PM -0400, Michael Kubovy wrote:
Dear Friends,
I not sure whether this is an Sweave or a
On Sep 6, 2008, at 6:07 PM, Ben Bolker wrote:
shalu shahlar at hotmail.com writes:
I am trying to define 25 vectors of varying lengths, say y1 to y25
in a loop,
and then store the results of some computations in them. My problem
is about
using some sort of concatenation for names. For
On Sep 6, 2008, at 4:24 PM, drflxms wrote:
Hello Mr. Burns, Hello Mr. Dwinseminus
snip
David, the input data.frame is the result of the reshape-command I
performed. I just copied it from the R-console into the e-mail. In
fact
the first column video is not part of the data, but needed
Please read the last line to every message to r-help. In particular
this question needs to include a cut down version of the data.
I'll take a guess at what it looks like:
library(zoo)
L - list(a = zoo(1:3), b = zoo(4:5))
plot(do.call(merge, L))
On Sat, Sep 6, 2008 at 4:47 PM, stephen sefick
On 06/09/2008 6:38 PM, Andrew Robinson wrote:
Hi Michael,
I think that beamer needs that frames that contain any verbatim text
or R output to be declared as fragile.
Try
\begin{frame}[fragile]
...
\end{frame}
The declaration I usually use is containsverbatim, but it probably
does the
the plot(do.call(merge, z.l)) works on the following data well. Is
there a way to get control of xlim so that it plots each individual
graph shows only the one day (figures the x axis on the range of the
data for each plot individually) and control labeling? Thanks in
advance, and sorry for not
Try xyplot.zoo with scale = list(relation = free)
specifying xlim as shown below:
library(zoo)
library(lattice)
zm - do.call(merge, z.l)
xlim - lapply(zm, function(x) range(time(na.omit(x
xyplot(zm, xlim = xlim, scale = list(relation = free))
On Sat, Sep 6, 2008 at 7:49 PM, stephen sefick
Ben Bolker wrote:
Martin Maechler maechler at stat.math.ethz.ch writes:
FEH == Frank E Harrell f.harrell at vanderbilt.edu
on Sat, 06 Sep 2008 07:19:33 -0500 writes:
FEH See http://biostat.mc.vanderbilt.edu/DynamitePlots for many reasons not
FEH to use dynamite plots.
Ah! Very
Hi,
I'm wanting to do a PCA on some data which is comprised of two different
groups (to see how well the groups are discriminated). Is there a way to
change the colour of the datapoints in a biplot so that I can easily see
which group is which (eg objects 1-100, red, 101-200, black).
Might be
Its a FAQ:
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-not-work_003f
On Sat, Sep 6, 2008 at 9:47 PM, Nathan Teuscher [EMAIL PROTECTED] wrote:
I have the following code that when executed from the command line
works properly and produces a proper PDF. When
65 matches
Mail list logo