> David Winsemius
> on Thu, 20 Oct 2011 01:51:28 -0400 writes:
> On Oct 19, 2011, at 11:29 PM, Alyse wrote:
>> Hello,
>>
>> I have a column in a data frame that need to be 10 digits long. As
>> such:
>>
>> Decimal.Year
>> 1 1994.25997
>> 2 1
Dear Folks--
All this seems to me to behave the way you expect, recognising that column b
is a factor:
> AA <- data.frame(a=3:4, b=c('x', 'y'))
> AA[1,]
a b
1 3 x
> as.numeric(AA[1,])
[1] 3 1
> AA[,2]
[1] x y
Levels: x y
> as.numeric(AA[,2])
[1] 1 2
> as.character(AA[,2])
[1] "x" "y"
But this se
late to the game but maybe this will help:
https://docs.google.com/viewer?a=v&q=cache:8PhCkZxP9zQJ:www.quantoid.net/Effects_package_4up.pdf+plot+effects+of+variables+in+R+nonlinear+arrows&hl=en&gl=us&pid=bl&srcid=ADGEEShUKOEcifzuWGxWvakh0yD4KtnLgBhLFvX5cCAkwewyQ75uznTw1OYybx6vrGuJflgMw6QYKGwuXQViN
Hi, I'm trying to overlay a geographical map with a heat map by following the
directions on
http://pages.stern.nyu.edu/~achinco/programming_examples/Example__PlotGeographicDensity.html.
However, the smaller my zoom level (the farther I zoom out), the more white
horizontal lines I have interspersed
let me clarify, i understand what differing x, y lengths mean. i understand
the concept of average effects, etc. i just don't understand how one would
fix it.
thanks.
--
View this message in context:
http://r.789695.n4.nabble.com/plotting-average-effects-tp3923982p3924003.html
Sent from the R h
hi... i am a phd student using r. i am having difficulty plotting average
effects. admittedly, i am not really understanding what each of the
commands mean so when i get the error i am not sure where the issue is.
here is my code... i will include the points at which there are errors
> da
So i have a simple function:
bmass=function(y){
weight=y$WT*y$MSTR
return(bio)
}
And want to apply to a whole bunch of rows in my data.frame:
final1=apply(final,1,yldbu)
BUT...recieve the following error:
"Error in y$WT : $ operator is invalid for atomic vectors"
However when i try:
> final[
I created a dendrogram (ddg0) using hclust in the usual way. I want to
collapse UP the tree in various ways, that is, from the leaves up to the
root. Optimally, I would give the id of a member of a final split in ddg0,
and return a new ddg1 with that split collapsed. Alt, I could give a depth
to co
Dear Jari,
Thank you for your quick reply and the time you have spent assisting with this
problem. Indeed the alias tool identifies one variable that when removed from
the capscale model solves the problem.
Once again greatly appreciate your assistance.
Regards
Steve Pawson
Scientist (Entomol
check this out
http://www.r-bloggers.com/pivot-tables-in-r/
--
View this message in context:
http://r.789695.n4.nabble.com/Aggregating-data-help-tp3923138p3923397.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mai
If it would help get any assistance with my issue, here's another
method I'm trying (using R sample data):
ggplot(mtcars, aes(disp)) +
geom_point(aes(y = mpg, colour = qsec))+
scale_colour_gradient(low="yellow", high="green")+
geom_point(aes(y = cyl, colour = qsec))+
scale_colour_gradient(low=
There's a great tutorial online that helped me out a lot - Lattice and Other
Graphics in R, by J H Maindonald at the Centre for Mathematics and Its
Applications at Australian National University.
http://maths.anu.edu.au/~johnm/r-book/2edn/xtras/rgraphics.pdf
I gave my lattice object a name, and t
Dear list members, dear Jay,
Well, I personally do not care about Revolutions Analytics selling their
products as this is also included into the idea of many open source
licences. Especially as Revolutions provide their packages to the
community and its is everybodies personal choice to buy th
Hello,
I am constructing a nomogram using cph and nomogram commands in Dr.
Harrell's Design/RMS package. The HR that I obtain for dichotomous and
categorical variables are identical to those that I obtain using STATA
stcox. However, the inter-quartile HR I obtain for continuous variables is
obviou
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 wan
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 dissertati
On Thu, Oct 20, 2011 at 3:42 PM, peter dalgaard 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 something.
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 = Da
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 houri
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 li
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 <
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 s
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')
# Us
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 wrote:
> 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
>
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
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 point
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
(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
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), dat
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, yma
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", "100
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
incl
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\00
(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 comma
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
>
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
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
extr
On Thu, Oct 20, 2011 at 7:29 PM, Bogaso Christofer
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 do analysis R using brow
Hi Uwe!
2011/10/20 Uwe Ligges :
> 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 I have so far:
> par(mar=
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 shape>1} I got my desired
result but when I take weibull(scale=1,shape=0.5) { i.e shape<1} it give
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)
Try Googling "R Portable"
Ken Hutchison
On Oct 20, 2554 BE, at 2:13 PM, jim holtman wrote:
> It runs fine off a flash drive.
>
> On Thu, Oct 20, 2011 at 1:29 PM, Bogaso Christofer
> wrote:
>> Dear all, the company I work for has Matlab installed for
>> statistical/mathematical calculations a
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 wrote:
> Dear Xu Jun,
>
> On Thu, 20 Oct 2011 14:41:30 -0400
>
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
Petar Milin 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 (http://
Dear Xu Jun,
On Thu, 20 Oct 2011 14:41:30 -0400
Xu Jun 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 polr
> funct
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 <- as.factor(ordwarm2
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 w
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 original
Hi,
On Thu, Oct 20, 2011 at 2:04 PM, Cem Girit 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 standard
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-
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 o
It runs fine off a flash drive.
On Thu, Oct 20, 2011 at 1:29 PM, Bogaso Christofer
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 possi
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
workin
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 Interv
Hi Uwe!
2011/10/20 Uwe Ligges :
> names.arg = rep(NA, nrow(file.codes))
Yes, that works beautifully. Danke schoen!
Henri-Paul
--
Henri-Paul Indiogine
Curriculum & Instruction
Texas A&M University
TutorFind Learning Centre
Email: hindiog...@gmail.com
Skype: hindiogine
Website: http://people
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 doe
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 y
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, x
I found it.
It is "ctr shift c"
Rui
On Thu, Oct 20, 2011 at 7:22 PM, Rui Esteves 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 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 wrote:
> Hi Rui,
>
> In the R terminal ctrl-c cancels the function, not the session.
>
> Cheers,
>
> Tsjerk
>
> On Oc
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" 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 that is taking ages.
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 de
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
nam
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 sugge
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 = FAL
On Thu, Oct 20, 2011 at 10:37 AM, peter dalgaard 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 multiplication
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
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
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 coul
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
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 f
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 s
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 Soc
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?
>
> > (sum(1000:1205))*(sum(1000
On Thu, Oct 20, 2011 at 7:50 AM, David Winsemius 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 * offe
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 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, jus
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 cor
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 think),
bbolker wrote:
>
> emmarosenfeld 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 resp
I got it! Thank you so much Jim!
On 20 October 2011 15:06, jim holtman wrote:
> try this: You were redfining 'transition' within the loop
>
> x <-
>
> c('A','D','F','H','N','A','C','H','F','D','F','F','H','K','G','D','D','B','N','K','O','V','S','S','F','H','J','U','K','T','H','S','R','G','S','R
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:
#
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
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 criterio
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 wrote:
> On Wed, Oct 19, 2011 at 4:08 PM, Sébastien Bihorel
> wrote:
> > Dear R-users,
> >
> > I would need some advices on the proper
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] "A
try this: You were redfining 'transition' within the loop
x <-
c('A','D','F','H','N','A','C','H','F','D','F','F','H','K','G','D','D','B','N','K','O','V','S','S','F','H','J','U','K','T','H','S','R','G','S','R','R','G','D','B','G','F','G','N','H','H','K','L','B','X','C','V','S','F','X','G','T','H',
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
hel
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
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 lis
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 aga
Thanks a lot for the answer, it works and really improves the speed of the
loop.
Best,
Sergio René
El , Eik Vettorazzi 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:
> >
> > I w
On Wed, Oct 19, 2011 at 4:08 PM, 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
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
>
>
> On 19.10.2011 22:08, Sébastien Bihorel wrote:
>
>> Dear R-users,
>>
>> I would nee
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 mail
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 = outlier(x1,l
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)
pos=which.
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
http://bioconductor.org/packages/2.8/bioc/vign
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
chan
1 - 100 of 113 matches
Mail list logo
