Thanks Josh,
I appreciate your comments. Merge is quite easy and I guess I should just go
that route. I don't have so many data points so the speed isn't really much
of an issue. I was just curious whether there was anything more elegant.
Yes, the function is not clear or ideal for
check the help archives.
hclust with method=ward might be what you are looking for
On Wed, Oct 19, 2011 at 2:43 PM, Jörg Reuter jo...@reuter.at wrote:
I am very new to R, so sorry that I ask stupid things.
I want compare a Matrix row by row and at the end I want to a Matrix with
the
My big problem is not ward, it is to calculate the Levenshtein-Distance.
check the help archives.
hclust with method=ward might be what you are looking for
On Wed, Oct 19, 2011 at 2:43 PM, Jörg Reuter jo...@reuter.at wrote:
I am very new to R, so sorry that I ask stupid things.
I want
Dear All:
I would like to know if there is plausible way to say to R to remove all
elements in the memory but the sequence. I have a code which makes a loop,
and what I want is after the programme has performed all the operation over
every ith element, to remove all the objects, expect the
Hi Sergio,
how about this:
rm(list=setdiff(ls(),z))
cheers.
Am 20.10.2011 11:00, schrieb Sergio René Araujo Enciso:
Dear All:
I would like to know if there is plausible way to say to R to remove all
elements in the memory but the sequence. I have a code which makes a loop,
and what I want
On 19.10.2011 22:08, Sébastien Bihorel wrote:
Dear R-users,
I would need some advices on the proper way to call a particular function.
This function is called scope.char and it is embedded in the step.gam
function from the gam package. I am trying to call scope.char directly in a
script but I
Steve Pawson Steve.Pawson at scionresearch.com writes:
My apologies for the delay in responding to your request for further
information I have been travelling for
work since you replied and have only just returned to email contact.
The output from the traceback is as follows
# This is the
Hello ALL!
I am running Linux, Fedora 15 64-bits, and R on it. I need to use
WinBUGS and R2WinBUGS, but as far as I read, WinBUGS is closed project,
to be continued with/as OpenBUGS. Thus, I have found R2OpenBUGS on
OpenBUGS Contributed Code (http://openbugs.info/w/UserContributedCode),
not on
Am 19.10.2011 17:33, schrieb Fong Chun Chan:
Hi,
Has anyone else run into this weird behaviour where the text in the plots
created using Cairo are always italicized. For example,
library(Cairo)
Cairo(file='cairo_created', type='pdf', dpi=100)
plot(1:10)
dev.off()
This produces the
i am new to the quantmod package . so if the answer is trivial please excuse
me. i want to study stock values within a day. i get current stock updates
using getQuotes and then want to produce usual quantmod graphs with that
values. also the graph should be able of adding technical indicators.
Dear R users:
I want to assign a global font family to a levelplot such as
par(family=serif) in traditional plotting. However, I can only change part of
the text using
panel.text(x,y,A,fontfamily=serif)
or scales=list(fontfamily=serif) to change the axis number font family.
How can I change
On Oct 20, 2011; 10:07am Jörg Reuter wrote:
I want compare a Matrix row by row and at the end I want to a Matrix with
the Levenshtein-Distance.
Jörg,
To begin with, try the following at the command prompt:
##
RSiteSearch(Levenshtein)
Shows, amongst other hits, that package vwr has a
Dear R community,
I tried to post a query about the use of the rose.diag function from the
CircStats package to the r-help forum a few days ago. It gave all appearances
of having been successful. It would have been my first post; I am correctly
recieving all the R-help emails, but but the
Hi!
Is There a way to manually create an affybatch object from qPCR array data?
--
View this message in context:
http://r.789695.n4.nabble.com/Creating-affybatch-objects-from-matrix-data-from-qPCR-array-tp3921559p3921559.html
Sent from the R help mailing list archive at Nabble.com.
Yes, I see many package. But the Problem is, the Package compare strings.
But there is a diffrent between the caracter 1 and 2 and the number 12.
All package I see compare every letter, but a number have many digits but
the digits make only sense if the function see them together. I tough to
On 10/20/2011 04:14 AM, ali_protocol wrote:
Hi!
Is There a way to manually create an affybatch object from qPCR array data?
Hi --
affybatch is for Affy probe information only; perhaps you mean an
ExpressionSet and the Bioconductor mailing list
Hi all,
I am doing analyses with haplo.stats in R. I know that when doing analysis
with haplo.glm, I have to define my outcome variable, i.e. if it is
binomial or gaussian.
When looking at interactions , eg. gender * geno or life events * geno do I
have to specify whether the independent factor
Hi,
I was looking for a solution to the same problem. I just found something and
thought it might be useful to share.
Suppose you look for the positions of the minimum value in matrix D. Using a
little euclide division solves the problem easily. Here is my code.
min(D)
Dv=as.vector(D)
Hi All,
I am working on the dataset in which some of the variables have more than
one observations with outliers .
I am using below mentioned sample script
library(outliers)
x1 - c(10, 10, 11, 12, 13, 14, 14, 10, 11, 13, 12, 13, 10, 19, 18, 17,
10099, 10099, 10098)
outlier_tf1 =
Can anyone give me links to reviews/comparisons of R with Amos for SEM? I
have found some but they are a little old (2009).
Ravi
--
View this message in context:
http://r.789695.n4.nabble.com/Structural-equation-modelling-in-R-compared-with-Amos-tp3921654p3921654.html
Sent from the R help
Thanks Uwe for reply.
The license is actually a problem: my package is distributed with a LGPL3
license, which is incompatible with GPL-2.
I'll try to find a work-around.
Sebastien
2011/10/20 Uwe Ligges lig...@statistik.tu-dortmund.de
On 19.10.2011 22:08, Sébastien Bihorel wrote:
Dear
On Wed, Oct 19, 2011 at 4:08 PM, Sébastien Bihorel pomc...@free.fr wrote:
Dear R-users,
I would need some advices on the proper way to call a particular function.
This function is called scope.char and it is embedded in the step.gam
function from the gam package. I am trying to call
Thanks a lot for the answer, it works and really improves the speed of the
loop.
Best,
Sergio René
El , Eik Vettorazzi e.vettora...@uke.de escribió:
Hi Sergio,
how about this:
rm(list=setdiff(ls(),z))
cheers.
Am 20.10.2011 11:00, schrieb Sergio René Araujo Enciso:
Dear All:
Hi Ravi,
Look at openmx, it uses R to do matrix optimization especially for SEM (though
more general too). It does not have the draw paths interactively feature, but
it is extremely powerful and flexible. I have used Amos, EQS, Mplus, Lisrel,
and OpenMx and I believe OpenMx is competitive
I see, I struggle too with SPSS. But I use this company for help:
http://www.ivoryresearch.com/custom-statistical-services-spss.php they are
really helpful
--
View this message in context:
http://r.789695.n4.nabble.com/need-help-on-read-spss-tp3893430p3921825.html
Sent from the R help mailing
Hi all,
I'd like to thank those who helped me with my previous loop function
question with agents/events. I have solved the problem with the advice from
this community.
I have now moved on to the next step, which requires me to find all the
transition pair within an event. A sample data and the
Have you considered simply changing your numbers into strings with
as.character()?
E.g.,
library(vwr)
levenshtein.distance(aba,cda)
num1 - 121; num2 - 341
levenshtein.distance(as.character(num1),as.character(num2))
I find that last line a little verbose to type, so I'd write a little
helper
try this: You were redfining 'transition' within the loop
x -
Here is another way of getting the lists
newtrans - lapply(split(mydata, mydata$agents), function(.agent){
+ .all - paste(.agent$actions, collapse = '')
+ .indx - embed(seq(nchar(.all)), 2)
+ substring(.all, .indx[, 2], .indx[, 1])
+ })
newtrans
$`2857`
[1] LB
$`293`
[1] AD DF FH
Thanks for the hack Gabor! That works just fine... (I guess until the
step.gam function is re-structured :D)
On Thu, Oct 20, 2011 at 9:20 AM, Gabor Grothendieck ggrothendi...@gmail.com
wrote:
On Wed, Oct 19, 2011 at 4:08 PM, Sébastien Bihorel pomc...@free.fr
wrote:
Dear R-users,
I would
Did you read the documentation for ?outlier. It clearly states that it
removes the single (possibly repeated) value with the largest distance
from the mean. That's only 10099 hereyou could perhaps apply the
function more than once or write your own outlier removal script using
whatever
On Oct 20, 2011, at 7:53 AM, Jörg Reuter wrote:
Yes, I see many package. But the Problem is, the Package compare
strings.
But there is a diffrent between the caracter 1 and 2 and the number
12.
All package I see compare every letter, but a number have many
digits but
the digits make only
Dear r-help listers,
I am using effects to produce an effect plot after the proportional
odds logistic regression model. There is no problem for me to estimate
the model, but when it comes to the graphing, I was stuck. see the
codes below:
I got it! Thank you so much Jim!
On 20 October 2011 15:06, jim holtman jholt...@gmail.com wrote:
try this: You were redfining 'transition' within the loop
x -
bbolker wrote:
emmarosenfeld emmarosenfeld at hotmail.co.uk writes:
Dear all,
I am having some problems trying to run a GLMM model with zero-inflation
using the alpha version of glmmADMB (0.6.4) using R (2.13.1) in Windows
and
I would greatly appreciate some help.
My count
Dear Ravi,
I would also suggest you to have a look at 'lavaan' (www.lavaan.org). It
has an extremely straightforward yet flexible model syntax and very good
documentation. Although its statistical capabilities are not comparable to
those of Mplus (which is the most powerful SEM-software I
On Oct 20, 2011, at 6:35 AM, richard hewitt wrote:
Dear R community,
I tried to post a query about the use of the rose.diag function from
the CircStats package to the r-help forum a few days ago. It gave
all appearances of having been successful. It would have been my
first post; I am
On Oct 20, 2011, at 1:45 AM, Lei Jiang wrote:
When using power operator ^, both the base and index are coerced to
type
real. The range of real is larger than integer.
However, an ordinary multiplication operator * offers it to two
integers (as
always).
So, to avoid the warning, just add
Try the memisc package which seems to have an independent
implementation from PSPP (at least looks that way when compiling on
linux) plus the ability to select variables before reading in the
entire dataset.
Paul Bivand
-
Paul Bivand
Head of
On Thu, Oct 20, 2011 at 7:50 AM, David Winsemius dwinsem...@comcast.net wrote:
On Oct 20, 2011, at 1:45 AM, Lei Jiang wrote:
When using power operator ^, both the base and index are coerced to type
real. The range of real is larger than integer.
However, an ordinary multiplication operator
On Oct 20, 2011, at 16:50 , David Winsemius wrote:
That does seem to be an infelicity that ought to be fixed. Using the infix
addition operator does that same sort of coercions, so why would one expect
the infix multiplication operator to refuse to do it?
Dear Xu Jun,
I'm not sure whether this is the source of the error, but it may help to spell
the xlevels argument correctly (it is not xlevles).
I hope this helps,
John
John Fox
Sen. William McMaster Prof. of Social Statistics
Department of
Lavaan is a nice package too, thanks for mentioning it. Given sufficient
understanding of matrices, I think OpenMx is arguably as or more powerful than
Mplus---it runs circles around Mplus in terms of flexibility, although Mplus
does what it does well, with (IMO) fairly sensible defaults, and
Hi!
I have 2 problems in drawing a stacked bar plot:
(1) This is a stacked bar plot with more than 100 bars next to each
other. So there should not be names at the bottom of the bars because
the bars are too narrow.I tried arg.names=NULL but that does not
work because R uses the row names
On 20.10.2011 18:03, Henri-Paul Indiogine wrote:
Hi!
I have 2 problems in drawing a stacked bar plot:
(1) This is a stacked bar plot with more than 100 bars next to each
other. So there should not be names at the bottom of the bars because
the bars are too narrow.I tried arg.names=NULL
Moohbear,
I may be able to help, in that I've had some graduate course work in plant
breeding and some experience with R.
However, I can't tell for certain with the information provided. What papers
that use SAS or SPSS are you citing?
The question I would ask, in order to determine if I
Need help with finding out approx over each row of a matrix. Here is the
setup:
years - matrix(c(1,2,3,1,2,3),nrow=2, ncol=3,byrow=TRUE)
rates - matrix(c(1,2,3,11,12,13), nrow = 2, ncol=3, byrow=TRUE)
points - matrix(c(1.5, 1.5), nrow=2, ncol=1, byrow=TRUE)
so basically i have above three
I am new to R, Sweave, and LaTeX. I am running StatET with Eclipse 3.7
on a Windows
7 x64 OS. StatET also includes the Sweave add-on, which I use to try and
compile a PDF document.
I also installed MiKTeX 2.9, because, previously, R could not execute texi2dvi
(presumably, because there was no TeX
On Thu, Oct 20, 2011 at 10:37 AM, peter dalgaard pda...@gmail.com wrote:
On Oct 20, 2011, at 16:50 , David Winsemius wrote:
That does seem to be an infelicity that ought to be fixed. Using the
infix addition operator does that same sort of coercions, so why would one
expect the infix
Hi Uwe! Thanks for the feedback. However, now the graph does not
draw the y-axis nor labels and gives the following error message:
Error in barplot.default(t(file.codes), beside = FALSE, names.arg = NA, :
incorrect number of names
This is the R code now:
barplot(t(file.codes), beside =
Dear all, the company I work for has Matlab installed for
statistical/mathematical calculations and really not ready to go with R
(even installing exe file for R). Therefore I was wondering is it possible
to do analysis R using browser like IE, without installing anything?
Thanks for your
On 20.10.2011 18:58, Henri-Paul Indiogine wrote:
Hi Uwe! Thanks for the feedback. However, now the graph does not
draw the y-axis nor labels and gives the following error message:
Error in barplot.default(t(file.codes), beside = FALSE, names.arg = NA, :
incorrect number of names
Hi,
This question seems very basic but I cannot find an answer on google.
I have a R session on a linux command line.
I called a function that is taking ages.
I want to cancel the function but without killing the R session.
What is the shortcut?
Thanks,
Rui
[[alternative HTML version
Hi Rui,
In the R terminal ctrl-c cancels the function, not the session.
Cheers,
Tsjerk
On Oct 20, 2011 7:16 PM, Rui Esteves ruimax...@gmail.com wrote:
Hi,
This question seems very basic but I cannot find an answer on google.
I have a R session on a linux command line.
I called a function
Hi Tsjerk,
In my command line it does not.
Maybe it because I am using linux.
That is my problem.
Thank you for answering,
Rui
On Thu, Oct 20, 2011 at 7:18 PM, Tsjerk Wassenaar tsje...@gmail.com wrote:
Hi Rui,
In the R terminal ctrl-c cancels the function, not the session.
Cheers,
I found it.
It is ctr shift c
Rui
On Thu, Oct 20, 2011 at 7:22 PM, Rui Esteves ruimax...@gmail.com wrote:
Hi Tsjerk,
In my command line it does not.
Maybe it because I am using linux.
That is my problem.
Thank you for answering,
Rui
On Thu, Oct 20, 2011 at 7:18 PM, Tsjerk Wassenaar
Dear Professor Fox,
Thanks a lot! That is an embarrassing error. After I cleaned up and
simplified my codes, and ran the following two lines:
myologit - polr(factor(warm) ~ yr89 + male + white + age + ed + prst,
+ data=ordwarm2, method=c(logistic))
plot(effect(age, myologit,
On 20.10.2011 19:22, Rui Esteves wrote:
Hi Tsjerk,
In my command line it does not.
Perhaps you arfe using a function that calls compiled code that does not
look for user interrupts?
try it on
repeat{1+1}
...
Uwe Ligges
Maybe it because I am using linux.
That is my problem.
Thank
I have been searching the web for an answer on whether there is a package for
group-based trajectory modeling in R. Something along the lines of what PROC
TRAJ (http://www.andrew.cmu.edu/user/bjones/index.htm) accomplishes in SAS.
Does anyone have any experience working with a package that
Hi Uwe!
2011/10/20 Uwe Ligges lig...@statistik.tu-dortmund.de:
names.arg = rep(NA, nrow(file.codes))
Yes, that works beautifully. Danke schoen!
Henri-Paul
--
Henri-Paul Indiogine
Curriculum Instruction
Texas AM University
TutorFind Learning Centre
Email: hindiog...@gmail.com
Skype:
Hello,
I need some results from the survival analysis of my data
that I do not know whether exist in Survival Package or how to obtain if
they do:
1. The Mean survival time
2. The standard error of the mean
3. Point and 95% Lower Upper Confidence
Dear Professor Fox,
I didn't include codes in between these two commands (polr and
plot(effect()). I am still trying to work it out. Maybe it's related
to how I coded some of the factor variables as I got some successful
runs after I changed the coding for some factor variables. I will keep
It runs fine off a flash drive.
On Thu, Oct 20, 2011 at 1:29 PM, Bogaso Christofer
bogaso.christo...@gmail.com wrote:
Dear all, the company I work for has Matlab installed for
statistical/mathematical calculations and really not ready to go with R
(even installing exe file for R). Therefore I
Hi-
I'm wondering if anyone can help me with my code. I'm coming up dry
when I try to get a p-value from the following code. If I make a
histogram of my resampled distribution, I find the difference between
by groups to be significant. I've ranked the data since I have
outliers in one
A bit new to R and I'm working on something a bit more challenging than I am
used to- so here's whats going on:
Data inputs: 9 different shapefiles (.shp) of different point locations
(lat, long)
shapefile bounding box (lat/long corner points (14, 1)
(15,1) (14, 2) (15,2))
Hi,
On Thu, Oct 20, 2011 at 2:04 PM, Cem Girit gi...@biopticon.com wrote:
Hello,
I need some results from the survival analysis of my data
that I do not know whether exist in Survival Package or how to obtain if
they do:
1. The Mean survival time
2. The
Dear Xu Jun,
It's really not possible for me to know the source of the error without a
complete, reproducible example, but if I had to guess, I'd guess that there's a
scoping problem of some sort, with the value of yr89 coming from somewhere that
you don't expect. How about fitting the
There have been posts in the past regarding similar questions, but many of
them looked dated. I am using xyplot to show variability within my
replicates, and so far so good, but I would like to refine the plot. I've
used the code below to graph my response variable against year (coded as a
factor
Dear Professor Fox,
Now I got it. It all comes from my unfamiliarity with the effect
function. I forgot the c part in the given.values option, plus it
looks like plot(effect()) does not like factor(warm) in the polr
function. So here are the two working lines:
ordwarm2$warm2 -
Dear Xu Jun,
On Thu, 20 Oct 2011 14:41:30 -0400
Xu Jun junx...@gmail.com wrote:
Dear Professor Fox,
Now I got it. It all comes from my unfamiliarity with the effect
function. I forgot the c part in the given.values option, plus it
looks like plot(effect()) does not like factor(warm) in the
Petar Milin pmilin at ff.uns.ac.rs writes:
Hello ALL!
I am running Linux, Fedora 15 64-bits, and R on it. I need to use
WinBUGS and R2WinBUGS, but as far as I read, WinBUGS is closed project,
to be continued with/as OpenBUGS. Thus, I have found R2OpenBUGS on
OpenBUGS Contributed Code
Hello,
I have a dataset with student performance on a math test. There are
multiple cases for each student (identified by id) and the concept as a
variable.
rtest
id test subject gradeconcept correct tested per_corr
year
11 83 Mathema 8 8.2.D 1
Dear Prof. Fox,
I just picked up R not long ago, and apologize that I am not that
familiar with some basics. I am trying to replicate what I can do with
Stata in R. Thanks for all your help!
On Thu, Oct 20, 2011 at 2:47 PM, John Fox j...@mcmaster.ca wrote:
Dear Xu Jun,
On Thu, 20 Oct 2011
Try Googling R Portable
Ken Hutchison
On Oct 20, 2554 BE, at 2:13 PM, jim holtman jholt...@gmail.com wrote:
It runs fine off a flash drive.
On Thu, Oct 20, 2011 at 1:29 PM, Bogaso Christofer
bogaso.christo...@gmail.com wrote:
Dear all, the company I work for has Matlab installed for
Can someone please help me out with this? The ggplot2 suggestion works
great but I've spent a few days trying to figure out how to plot 2
variables with it and I'm stuck. Here's my example code:
library(ggplot2)
#Here's the 1st plot
x-rnorm(100)
y-rnorm(100)
z-rnorm(100)
d - data.frame(x,y,z)
Dear Sir/madam,
I'm getting a problem with a R-code which calculate Fisher Information
Matrix for Hybrid Censored Weibull Distribution. My problem is that:
when I take weibull(scale=1,shape=2) { i.e shape1} I got my desired
result but when I take weibull(scale=1,shape=0.5) { i.e shape1} it gives
Hi Uwe!
2011/10/20 Uwe Ligges lig...@statistik.tu-dortmund.de:
arrange it outside by, e.g. increasing the size of margins (see argument
mar in ?par) and place a separate legend (see ?legend) into the margins
(see xps argument in ?par).
I could not find 'xps', do you mean 'xpd'?
This is what
On Thu, Oct 20, 2011 at 7:29 PM, Bogaso Christofer
bogaso.christo...@gmail.com wrote:
Dear all, the company I work for has Matlab installed for
statistical/mathematical calculations and really not ready to go with R
(even installing exe file for R). Therefore I was wondering is it possible
to
Hello all,
I'm struggling to grasp how I might use lapply() instead of looping to
run approx() on a list consisting of multiple arrays - each of equal
dimension.
But simpler than that, I haven't been able to successfully apply
approx() to an array, unless I loop through the third dimension and
Generally getting intraday/real-time data requires some sort of (paid)
source, but if you are willing to just strip free quotes off the
internet non-stop, perhaps something like this:
while(TRUE) {
cat( file = FileName.txt, c(getQuote(YHOO), recursive = T), \n,
append = T)
}
# You could get
Hi,
Please send your information to the r-help list, not just to me, but do
note that the list is plain-text only.
But surely all you are looking for is:
dt-c(37,41,40,38,38,37,44,45,48,43,48,46,54,60,32,45,55,62,42,62,62,62,47,42,59,43,60,60,51,43,50,51,47,42,47,51)
mean(dt)
[1] 48.16667
(Iske - matrix(c(1, 1, 1, 2, 2, 2, 10, 1, 1, 5,
1,2,2,2,1,1,1,4,4,4,4,3,3,3,30,3,3,3,3,3,3,3,5,4,4,4,4,4,44,2,2,2,2,4,5,5,5,5,5,5,5,5,95,1,2),
ncol = 5))
numtochar - function(Zahl){
text - rawToChar(as.raw(Iske[1,]))}
(Iske.char[]-apply(Iske,1,numtochar))
I have a little problem with the command
Does this work?
apply(Iske, 1, function(x) rawToChar(as.raw(x)))
apply(Iske, 1, function(x) rawToChar(as.raw(x)))
[1] \001\002\003\004\005 \001\002\003\004\005 \001\002\036\004\005
[4] \002\001\003\004\005 \002\001\003\004\005 \002\001\003,\005
[7] \n\004\003\002\005 \001\004\003\002\005
Hi,
It looks like you are trying to manually bootstrap. Take a look at:
require(boot)
?boot
as an added advantage of using boot instead of trying to do it
manually, you can easily parallelize. In fact, if you are using one
of the pre-release versions of 2.14.0, the new parallel package is
I have a table that looks like this:
structure(list(speed = c(3,9,4,8,7,6), C = c(0.697, 0.011, 0.015, 0.012, 0.018,
0.019), house = c(1,
1, 1, 1, 1, 1), date = c(719, 1027, 1027, 1027, 1030, 1030),
hour = c(18, 8, 8, 8, 11, 11), id = c(1000, 1,
10001, 10002, 10003, 10004)),
Not doing much work with spacial stats or shapefiles, I can't help in
too much detail, but here are some R commands that might help for each
part:
a.
# This will help you pick random points within your bounded box
runif2d - function(n, xmin, xmax, ymin, ymax){
stopifnot(all(xmax xmin, ymax
I'd dont quite get what you are asking, but here's my best guess and
you can tweak it go get what you need.
I'd do it in two passes since you want two rather unrelated objects.
structure(list(speed = c(3,9,4,8,7,6), C = c(0.697, 0.011, 0.015,
0.012, 0.018, 0.019), house = c(1,
1, 1, 1, 1, 1),
On Oct 20, 2011, at 5:05 PM, Sarah Goslee wrote:
Hi,
Please send your information to the r-help list, not just to me, but
do
note that the list is plain-text only.
But surely all you are looking for is:
dt-
c
On 21/10/11 06:32, magono wrote:
A bit new to R and I'm working on something a bit more challenging than I am
used to- so here's whats going on:
Data inputs: 9 different shapefiles (.shp) of different point locations
(lat, long)
shapefile bounding box (lat/long corner
On Oct 20, 2011, at 23:48 , Joshua Wiley wrote:
Hi,
It looks like you are trying to manually bootstrap.
Nope. It's a manually performed approximate Wilcoxon test. Which is fair enough
if the object is to learn something. (Notice, however, that the ExactRankTests
package eats this sort
Hi,
Just include pch=c(1,2,3), col=c('red','blue','yellow') or your
choices in xyplot
Weidong Gu
On Thu, Oct 20, 2011 at 2:15 PM, wisc_maier cmai...@wisc.edu wrote:
There have been posts in the past regarding similar questions, but many of
them looked dated. I am using xyplot to show
Hi:
Here's a way using the reshape2 package.
library('reshape2')
rsub - subset(rtest, concept %in% c('8.2.D', '8.3.A', '8.3.B'))
# want year ahead of concept in the variable list
rsub - rsub[, c(1:4, 9, 5:8)]
cast(rsub, id + test + subject + grade + year ~ concept, value_var = 'per_corr')
#
I'm not seeing how to plot the quantities associated with two values of a
factor by reading the ?xyplot help page and Deepayan's book. Perhaps I need
to generate many different subsets from the data frame, but that would
require _many_ new data frames.
The structure of the data frame (for a
AFAIK, you can't 'add' two ggplot2 graphs together; the problem in
this case is that the two color scales would clash. If you're willing
to discretize the z values, then you could pull it off. Here's an
example:
d - data.frame(x = rnorm(100), y = rnorm(100), z = factor(1 +
(rnorm(100) 0)))
d1 -
Jannis,
I'm not complete sure I understand your first point, but maybe someone
from REvolution will weigh in. Nobody is forcing anyone to purchase
any products, and there are attractive alternatives such as the CRAN R
and R Studio (to name two). This issue has arisen many times of the
various
Hi:
Here's one way with the plyr package. Using ds as the name of your
data frame (thank you for the dput and clear description of what you
wanted, BTW),
library('plyr')
ddply(ds, .(date), mutate, minspd = min(speed), Cmin =
C[which.min(speed)], diff = C - Cmin)
speed C house date hour
Hi Rich
Without a dataset I am not sure what you need.
Further down the track you will need to plot what you finally need
and this may help you decide what you need now
for finals try
library(lattice)
library(latticeExtra)
useOuterStrips(
xyplot(variable1 + variable2 ~ Date|Site, data =
On Thu, Oct 20, 2011 at 3:42 PM, peter dalgaard pda...@gmail.com wrote:
On Oct 20, 2011, at 23:48 , Joshua Wiley wrote:
Hi,
It looks like you are trying to manually bootstrap.
Nope. It's a manually performed approximate Wilcoxon test. Which is fair
enough if the object is to learn
Hi!
I am trying to use ggplot2 to create a stacked bar plot. Previously I
tried using barplot() but gave up because of problems with the
positioning of the legend and other appearance problems. I am now
trying to learn ggplot2 and use it for all the plots that I need to
create for my
Hi,
I am very happy. My problems are solved without one little thing:
(Iske - matrix(c(1, 1, 1, 2, 2, 2, 1, 1, 1, 5, 1, 2, 2, 2, 1, 1, 1,
4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 5, 4, 4, 4, 4, 4, 4, 2,
2, 2, 2, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 1, 2), ncol = 5)) #My Matrix
Iske- Iske+33 #I want
100 matches
Mail list logo