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. Whe
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 si
I have the following code that when executed from the command line
works properly and produces a proper PDF. When the script is executed,
the PDF produced is considered corrupt. I am using R 2.7.2 on Mac OSX
10.5.4. Thank you in advance for the help!
library(lattice)
pdf(file="CLDiag2.pdf")
Ben Bolker wrote:
Martin Maechler stat.math.ethz.ch> writes:
"FEH" == Frank E Harrell 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! Let's hope
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 sefi
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 p
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 sa
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 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 neede
On Sep 6, 2008, at 6:07 PM, Ben Bolker wrote:
shalu 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, in
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 bea
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
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))
> apply(x,3,fun
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
bbayes[[i]]<-
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}
\b
shalu 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 through y
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.
Simil
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
https://stat
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 d
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
__
R-help@r-p
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 th
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
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 A
jpl math.unl.edu> writes:
>
>
> Hi,
>
> I am trying to fit a curve to data. My command line is:
>
>
model10=nls(offspring~((A*c^k)/gamma(k))*
((degdays-alpha)^(k-1))*exp(-c*(degdays-alpha)),
> start=list(A=30,k=2,c=.018,alpha=131))
>
try trace=TRUE to see where the function is try
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,
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 o
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 data.
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 dat
Hi Felix,
>> My problem is now, that I analyse data.frames with an unknown count of
>> columns. So to get rid of the first and last column for the "pred"
>> variable
>> and to select the last column for the "ref" variable, ...
Doubtless there are other routes. Generally I use ?length to get the
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 duplic
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 appreci
Hi Muechen,
OK !!! Thank you very much
Ricardo
Muenchen, Robert A (Bob) wrote:
>
> 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, SP
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
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"),heade
#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
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 = den
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
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, Septe
Martin Maechler stat.math.ethz.ch> writes:
>
> > "FEH" == Frank E Harrell 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 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)
return(c(unlist(x.freq.sel
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
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 the
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] sub
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
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/
http://zoonek2.
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 2
> "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 (Ad
> "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 w
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
"freq
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, b
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 versi
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_
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
Here is a start. You can delete the zeros:
> 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))
> apply(x,3,function(.mat){
+ rows
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 ba
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
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 linea
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 ti
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 ba
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
__
R
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
__
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
https://s
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, ca
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
s
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 how
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.
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,
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
https://stat.ethz.ch/mailm
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 so
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
h
69 matches
Mail list logo
