Dear R helpers,
I try to find the model parameters using mle2 (bbmle package). As I try to
optimize the likelihood function the following error message occurs:
Error in grad.default(objectivefunction, coef) :
function returns NA at
Hi,
How can I add to LDA ggplot axes the Percentages of variance explained?
Script:
/require(MASS)
require(ggplot2)
iris.lda-lda(Species ~ Sepal.Length + Sepal.Width + Petal.Length +
Petal.Width, data = iris)
datPred-data.frame(Species=predict(iris.lda)$class,predict(iris.lda)$x)
With R-3.0.1
Loading required package: Rmpi
Failed with error: ‘package ‘Rmpi’ was built before R 3.0.0: please re-install
it’
And when I try to reinstall Rmpi, I get this
after a whole bunch of 'yes's
checking mpi.h usability... no
checking mpi.h presence... no
checking for mpi.h...
Hello,
Maybe this link might help:
http://www.stats.uwo.ca/faculty/yu/Rmpi/install.htm
Regards,
Pascal
2013/8/13 Patrick Connolly p_conno...@slingshot.co.nz
With R-3.0.1
Loading required package: Rmpi
Failed with error: package Rmpi was built before R 3.0.0: please
re-install it
Dear
R-users,
Iâve been
looking at the lmer function (lme4 package) in order to set up a mixed linear
model
and something about the syntax of the random effects eludes me. Iâd like a
hand
with understanding a specific point, if someone does master this functionâ¦
Letâs say
that I have
Hi Gergely!
This sounds so exciting - so R is turning 20 years old? How did you set up your
R users' group? What are the best practices in going about to set one up?
I would be keen on establishing one here in Windhoek, Namibia.
Pancho Mulongeni
Research Assistant
PharmAccess Foundation
1
G'day
I try do compute some F-statistics of a singular spectrum analysis of a
timeseries sv
I run:
require(Rssa)
s - ssa(sv)
summary(sv)
Min. 1st Qu. MedianMean 3rd Qu.Max.
-4.238 2.761 6.594 6.324 10.410 15.180
r1 - reconstruct(s,groups = list(1:5))
r2 -
Hi Pancho,
there are already a bunch of R User Groups around the world:
http://rwiki.sciviews.org/doku.php?id=rugs:r_user_groups
The Revolution Analytics guys has posted some tips on how to found one (
http://www.revolutionanalytics.com/news-events/r-user-group/how-to-start-r-user-group.php)
and
Hello,
r1$df and r2$df don't exist.
Regards,
Pascal
2013/8/13 Ingo Wardinski i...@gfz-potsdam.de
G'day
I try do compute some F-statistics of a singular spectrum analysis of a
timeseries sv
I run:
require(Rssa)
s - ssa(sv)
summary(sv)
Min. 1st Qu. MedianMean 3rd Qu.Max.
1) Why use Nelder-Mead with optimx when it is an optim() function. You
are going from New York to Philadelphia via Beijing because of the extra
overhead. The NM method is there for convenience in comparisons.
2) NM cannot work with NA when it wants to compute the centroid of
points and search
Hi all,
Any insight into the code below would be appreciated - I don't
understand why two methods which I think should have equal distance
from the call don't.
Thanks!
Hadley
# Create simple class hierarchy
setClass(A, NULL)
setClass(B, A)
a - new(A)
b - new(B)
setGeneric(f, function(x, y)
That's the primary reason for the plot: so that you can look and think.
The test statistic is based on whether a LS line fit to the plot has zero slope. For
larger data sets you can sometimes have a significant p-value but good agreement with
proportional hazards. It's much like an example
Hello!
I have a very large experimental design space (all possible combinations of
all possible levels of several factors). For example, 'allcand' below has
1,875,000 possible combinations of 9 factors.
allcand-expand.grid(a1=as.factor(1:5),a2=as.factor(1:5),a3=as.factor(1:5),a4=as.factor(1:5),
Hadley,
The class AB inherits from A and from B, but B already inherits from class A.
So actually you only have an object of class B in your object of class AB. When
you call the function f R looks for a method f for AB objects. It does not find
such a method and looks for a method of the
Robert U tacsunday at yahoo.fr writes:
Dear
R-users,
[snip]
This question probably belongs on r-sig-mixed-mod...@r-project.org .
Followups there, please.
Let's say that I have 2 random effects, A (e.g. species, k=2) and B
(e.g. individuals, n=100). I made some research about model
Try this
S2 - data.frame(Group=rep(S, length(S)), Cat=factor(S))
B2 - data.frame(Group=rep(B, length(B)), Cat=factor(B))
V2 - data.frame(Group=rep(V, length(V)), Cat=factor(V))
table(rbind(S2, B2, V2))
-
David L Carlson
Associate Professor of Anthropology
The class AB inherits from A and from B, but B already inherits from class A.
So actually you only have an object of class B in your object of class AB.
When you call the function f R looks for a method f for AB objects. It does
not find such a method and looks for a method of the object
Thank you for your response, Terry.
To put the discussion into perspective, my data set is quite large with
over 160,000 samples and 38 variables. The event is true for all samples in
this dataset. The distribution is zero-inflated (i.e. most events occur at
time = 0).
The result of the cox.zph
Oops! Ok So I have this file:
SampleName Individual Age Gender
1 4 80 M
2 15 56 F
3 1 75 F
4 15 56 F
5 2 58 F
6 4 80 M
And I want to pull out paired samples, so the resulting file would look
something like this:
SampleName Individual Age Gender
1 4 80 M
2 15 56 F
4 15 56 F
6 4 80 M
.kripa
Hi experts,
I have a dataframe with 100k+ records. it has a key/id column and 25 code
columns. I would like to restructure it having a row for each code column.
I have a structure like this (used dput):
structure(list(DSYSRTKY = structure(c(1L, 2L, 3L, 3L, 4L, 4L), .Names =
c(1,
2, 3, 4, 5, 6),
Hi,
Yes, I have searched stack overflow.
My issue is to simply change coloring in boxes and legend in my bwplot. I have
done this many times in lattice, but now I have been tweaking the plot somewhat
and I can no longer apply the color changes.
I would really appreciate some help.
A. Zakrisson
I am teaching a summercourse this and the next week where the students are
using R.
We have downloaded and use the new version R.3.0. It has worked perfectly until
today where some of the basic functions have started NOT to work.
Examples are sd() and lm ()
The message we get is Error:
Thank you David.
I had to sort the data afterwards for it to work:
S-1:86
B-1:15
V-1:45
S2 - data.frame(Group=rep(S, length(S)), Cat=factor(S))
B2 - data.frame(Group=rep(B, length(B)), Cat=factor(B))
V2 - data.frame(Group=rep(V, length(V)), Cat=factor(V))
table(rbind(S2, B2, V2))
If you take an example which works with slots,
setClass(A, representation(a = numeric)
setClass(B, contains = c(A), representation(b = numeric))
a - new(A, a = 2)
b - new(B, a = 3, b = 2)
setClass(AB, contains = c(A, B))
new(AB, a = 2, b = 3)
You see, that there is only one @a slot, the one
Hi,
The conditions are still not clear.
dat2- dat1[dat1$Individual%in% reps,]
dat2
# SameName Individual Age Gender
#1 1 4 80 M
#2 2 15 56 F
#4 4 15 56 F
#6 6 4 80 M
A.K.
I manipulated the code you sent and it works perfectly, thanks!
.kripa
Date: Tue, 13 Aug 2013 08:10:53 -0700
From: smartpink...@yahoo.com
Subject: Re: [R] pulling out pairs from data frame
To: kripa...@hotmail.com
CC: r-help@r-project.org
Hi,
The conditions are still not clear.
Hi again,
I have a question on some outliers that I have in my response variable (wich
are bird counts). At the beginning I did not drop them
out because they are part of the normal counts and I considered them
ecologically correct.
However, I
tried some of the same models without
HI,
Your desired output is not clear. May be this helps:
#dat1 is the dataset
dat1$ID- 1:nrow(dat1)
library(reshape2)
res1-melt(dat1,id.vars=c(ID,DSYSRTKY))
res1$value-res1$value!=
res1[,2]- as.integer(as.character(res1[,2]))
res1[,3]-as.character(res1[,3])
colnames(res1)[3:4]-c(CODE,PRIMARY)
Thanks for your interest and prompt answer!
What I try to estimate is the correlation of one bird species counts with a set
of environmental parameters. The count data are zero-inflated and
overdispersed. I am modeling with hurdle-negative binomial-mixed effects.
The results are very difficult
I do not know what you are exactly estimating, but if it is about count models
and the model fit gets better when you drop the outliers, it does not say, that
the model is now more correct. It just says, if the data were without the
outliers, this model would fit good.
Overdispersion in count
?duplicated
yourframe[!duplicated(yourframe)$Individual,]
-- Bert
On Tue, Aug 13, 2013 at 8:12 AM, Kripa R kripa...@hotmail.com wrote:
I manipulated the code you sent and it works perfectly, thanks!
.kripa
Date: Tue, 13 Aug 2013 08:10:53 -0700
From: smartpink...@yahoo.com
Subject: Re:
The central question is: What caused the 3 unusual values? What is
their scientific relevance? Only you can answer that, not us.
-- Bert
On Tue, Aug 13, 2013 at 8:51 AM, Marta Lomas lomasv...@hotmail.com wrote:
Thanks for your interest and prompt answer!
What I try to estimate is the
Sorry. Typo. Corrected version is:
yourframe[!duplicated(yourframe$Individual),]
-- Bert
On Tue, Aug 13, 2013 at 9:05 AM, Bert Gunter bgun...@gene.com wrote:
?duplicated
yourframe[!duplicated(yourframe)$Individual,]
-- Bert
On Tue, Aug 13, 2013 at 8:12 AM, Kripa R kripa...@hotmail.com
Dear R users,
occasionally I find .Rhistory and/or .RData files cluttered around in my
file structure. Is there a way to tell R not to save such files? Or to
use one central location where to save them (if they are of any use)? I
have looked through options() to no avail.
Cheers
Jannis
I don't see a question in what you wrote. Your graph has some similarities
to some of my examples.
Please look at the demo in the HH package
## install.packages(HH) ## if necessary
library(HH)
demo(bwplot.examples, package=HH)
Rich
On Tue, Aug 13, 2013 at 10:00 AM, Anna Zakrisson Braeunlich
I have to process a set of survey data with questions that are formatted
like this;
1) Pick your top three breeds (pick 3)
1 Rottweiler
2 Pit Bull
3 German Shepard
4 Poodle
5 Border Collie
6 Dalmation
7 Mixed Breed
and the answers are formatted like this:
Respondent, Question1
I think I understand your question. You need to make sure that you are
setting the right parameters in your theme. Use trellis.par.get() to have
a look at the MANY possible settings. For example, in your case, to have
the boxplots and rectangles be the same color:
my.theme - list(
Thanks Bert!
I think they are relatively important. What I am doing is comparing 2003 with
2013 distribution and use of this species in an specific sampled area. They are
currently way lower numbers than in 2003, however in both years the data are
zero inflated. Most of the outliers are in 2003
R-3.0.1 (use all digits for describing an R version) is not the problem.
Most likely you are masking a function
or something like that. When you started the R session, did you get a
message about restoring a previous session?
If so, then close R, find the directory in which you were working and
Dear R users,
I've just updated the WriteXLS package (on R 3.0.1) and I now have an error
when exporting a data.frame with the argument Encoding=latin1. For
example, these two lines work:
library(WriteXLS)
WriteXLS(iris, iris.xls)
whereas these ones don't work:
library(WriteXLS)
Bert,
dat1-structure(list(SameName = 1:6, Individual = c(4L, 15L, 1L, 15L,
2L, 4L), Age = c(80L, 56L, 75L, 56L, 58L, 80L), Gender = c(M,
F, F, F, F, M)), .Names = c(SameName, Individual,
Age, Gender), class = data.frame, row.names = c(NA, -6L
))
Your solution gives:
Dear r users,
what could cause such an error:
internal error -3 in R_decompress1
unfortunately the error kills all my usual error catching mechanisms an
appears on a remote cluster so I can not really tell you which command
etc is causing it.
Thanks for any hints on where to dig for the
I've written the following function
CoursePrep - function (Source, SaveName) {
Clean$TERM - as.factor(Clean$TERM)
Clean$INST_NUM - as.factor(Clean$INST_NUM)
Clean$zGrade - with(Clean, ave(GRADE., list(TERM, INST_NUM), FUN =
scale))
write.csv(Clean,paste(SaveName, csv, sep =.), row.names
On 13/08/2013 18:47, Jannis wrote:
Dear r users,
what could cause such an error:
internal error -3 in R_decompress1
unfortunately the error kills all my usual error catching mechanisms an
appears on a remote cluster so I can not really tell you which command
etc is causing it.
It is a
Hi,
You could try:
lapply(split(Clean,list(Clean$TERM,Clean$INST_NUM)),function(x)
shapiro.test(x$GRADE))
A.K.
- Original Message -
From: Robert Lynch robert.b.ly...@gmail.com
To: r-help@r-project.org
Cc:
Sent: Tuesday, August 13, 2013 1:46 PM
Subject: [R] ave function
I've written
Yes, you're right.
So I guess you should match on duplicated values, something like (untested)
with(dat1, dat1[individual %in% individual[duplicated(individual)],]
which is presumably essentially what you gave.
-- Bert
On Tue, Aug 13, 2013 at 10:41 AM, arun smartpink...@yahoo.com wrote:
Hi,
My desired output for my sample!! using dput():
structure(list(ID = c(1, 2, 3, 4, 5, 6, 7, 8,
9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30,
31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41,
42, 43, 44, 45, 46, 47, 48), DSYSRTKY = c(10005,
10005,
Dear All,
I am running a negative binomial model in R using the package pscl in oder to
estimate bed sediment movements versus river discharge. Currently we have
deployed 4 different plates to test if a combination of more than one plate
would better describe the sediment movements when the
Hi,
Many apologies for the simplicity (hopefully!) of this request - I can't
find it on the forum, but it may have been asked in the past.
I have a data frame consisting of ~2000 rows. I simply want to take the
average of the first 6, then the next 6, then the next 6 until the end of
the table.
sir i have installed plant breeding library well. But when i import the
file in R and give command
* data(fulldial)
Warning message:
In data(fulldial) : data set fulldial not found*
above warning message is found please guide me. my data is under
MALE FEMALE YIELD 1 1 53.333 1 2 52.333 1
The following should help:
What does R ask you each time you quit R? Answer no.
Start R with
R --no-save
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 8/13/13 9:15 AM, Jannis bt_jan...@yahoo.de wrote:
Dear R
On 08/13/2013 12:17 PM, Walter Anderson wrote:
I have to process a set of survey data with questions that are formatted
like this;
1) Pick your top three breeds (pick 3)
1 Rottweiler
2 Pit Bull
3 German Shepard
4 Poodle
5 Border Collie
6 Dalmation
7 Mixed Breed
and the
According to your first post,
NewDataFrame - data.frame(ID=integer(), DSYSRTKY=integer(), CODE=character(),
PRIMAIRY=logical())
The new output dataset: Out1
str(Out1)
'data.frame': 48 obs. of 4 variables:
$ ID : chr 1 2 3 4 ...
$ DSYSRTKY: chr 10005 10005 10005
What about something like this:
tmp - data.frame(var1 = rnorm(36), ind = gl(6,6))
with(tmp, tapply(var1, ind, mean))
You can see that your version of
mean(tmp[1:6,c(var1)])
gives the same as mine for the first 6 rows.
-Original Message-
From: r-help-boun...@r-project.org
From: Kevin E. Thorpe [mailto:kevin.tho...@utoronto.ca]
Sent: Monday, August 12, 2013 11:00 AM
Subject: Re: [R] Memory limit on Linux?
What does ulimit -a report on both of these machines?
Greetings,
Sorry for the delay. Other fires demanded more attention...
For the system in which
I have a dataFrame
sID - c(a, 1,2,3, b, 4,5,6)
rID - c(shr1125, bwr331, bwr330, vjhr1022)
tmp - data.frame(cbind(sID,rID))
but I need to split tmp$sID into three different columns, filling locations
where tmp$sID has only one value with NA.
I can split tmp$sID by the comma
tmp.1 -
On 08/13/2013 03:06 PM, Stackpole, Chris wrote:
From: Kevin E. Thorpe [mailto:kevin.tho...@utoronto.ca]
Sent: Monday, August 12, 2013 11:00 AM
Subject: Re: [R] Memory limit on Linux?
What does ulimit -a report on both of these machines?
Greetings,
Sorry for the delay. Other fires demanded
Lauria: For historical reasons the logistic regression (binomial with
logit link) model portion of a zero-inflated count model is usually
structured to predict the probability of the 0 counts rather than the
nonzero (=1) counts so the coefficients will be the negative of what you
expect based on
You could also try:
##Out1 is the output dataset
Out1$PRIMAIRY-as.logical(Out1$PRIMAIRY) #changing the class
#dat1 input dataset
vec1- paste(dat1[,1],dat1[,2],colnames(dat1)[2],sep=.)
res2-reshape(dat1,idvar=newCol,varying=list(2:26),direction=long)
res3-res2[order(res2[,4]),]
res4-
Try,
sID - c(a, 1,2,3, b, 4,5,6)
tmp1 - strsplit(sID,',')
tmp2 - lapply(tmp1,
function(x) if (length(x)==1) c('','',x) else x )
tmp3 - matrix(unlist(tmp2),ncol=3, byrow=TRUE)
rID - c(shr1125, bwr331, bwr330, vjhr1022)
newdf - data.frame(cbind(tmp3,rID))
You'll need to name
Hello.
I am using the aovp() function from the library lmPerm with one factor
(group: 3 levels) controlling for 2 covariates. I now want to conduct a
post-hoc test using the same model. Unfortunately, I did not find an
appropriate test which works with 2 covariates. I would be grateful for any
On 13-08-2013, at 18:46, jsf1982 jamie.free...@ucl.ac.uk wrote:
Hi,
Many apologies for the simplicity (hopefully!) of this request - I can't
find it on the forum, but it may have been asked in the past.
I have a data frame consisting of ~2000 rows. I simply want to take the
average of the
Hi,
You could try:
set.seed(24)
dat1- as.data.frame(matrix(sample(1:50,29*6,replace=TRUE),ncol=6))
((seq_len(nrow(dat1))-1)%/%6)+1
# [1] 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5
#For a particular column:
aggregate(dat1[,5],list(((seq_len(nrow(dat1))-1)%/%6)+1),FUN=mean)
#
Hi,
You could try:
tmp[,1]- as.character(tmp[,1])
tmp[,1][-grep(,,tmp[,1])]-paste0(,,,tmp[,1][-grep(,,tmp[,1])])
tmp2-data.frame(read.table(text=tmp[,1],sep=,,header=FALSE,stringsAsFactors=FALSE),rID=tmp[,2],stringsAsFactors=FALSE)
colnames(tmp2)[1:3]-paste(sID,letters[1:3],sep=.)
tmp2
#
Hi everyone:
I am running a code in R and I get the following message after using
large files (files larger than 2GB):
Runtime error!
this application has requested the Runtime to terminate it in an usual way.
Please contact the application's support team for more information
Another person
I have a set of survey data where I have answers to identify preference
of three categories using three questions
1) a or b?
2) b or c?
3) a or c?
and want to obtain weights for each of the preferences
something like X(a) + Y(b) + Z(c) = 100%
I am at a loss how how to calculate this from
It would seem that in going to a 64 bit architecture you have not escaped your
memory problems. Such problems are highly varied in details, so you would need
to be much more specific about how you are encountering this problem before
anyone could help. Read the Posting Guide and make a
On 08/13/2013 11:41 AM, Siraaj Khandkar wrote:
On 08/13/2013 12:17 PM, Walter Anderson wrote:
I have to process a set of survey data with questions that are formatted
like this;
1) Pick your top three breeds (pick 3)
1 Rottweiler
2 Pit Bull
3 German Shepard
4 Poodle
5 Border
Hi,
You could try:
dat2- read.table(text='
Respondent, Question1
1, 1,4,7
2, 2,7,5
3, 6,3,5
4,
',sep=,,header=TRUE,stringsAsFactors=FALSE)
library(stringr)
dat2New-cbind(dat2,do.call(rbind,lapply(
str_split(str_trim(dat2[,2]),,),as.numeric)))
colnames(dat2New)[3:5]-
Hi Hadley,
I suspect that the dispatch algorithm doesn't realize that selection
is ambiguous in your example. For 2 reasons:
(1) When it does realize it, it notifies the user:
setClass(A, NULL)
setGeneric(f, function(x, y) standardGeneric(f))
setMethod(f, signature(A,
I had a similar problem and found when looking
inside one of the lattice functions that the
legend colours are controlled by the superpose
series eg superpose.line, superpose.polygon etc
in trellis.par.set/get or par.settings
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
Among the many solutions, here is the one using phenology package:
library(phenology)
plot_errbar(1:100, rnorm(100, 1, 2),
xlab=axe x, ylab=axe y, bty=n, xlim=c(1,100),
errbar.x=2, errbar.y=rnorm(100, 1, 0.1))
or
x - 1:100
plot_errbar(x=1:100, rnorm(100, 1,
Hi Listers
I have been trying to import a .ps graphic file into R using the grImport
package but I keep getting the following error message
Error in PostScriptTrace(fish.ps) :
status 127 in running command 'gswin32c.exe -q -dBATCH -dNOPAUSE
-sDEVICE=pswrite
Hello,
What is the result of sessionInfo()?
Regards,
Pascal
2013/8/14 Andrew Halford andrew.half...@gmail.com
Hi Listers
I have been trying to import a .ps graphic file into R using the grImport
package but I keep getting the following error message
Error in PostScriptTrace(fish.ps) :
Dear R Helpers,
I would like to pull the CIK number from the web page
http://www.sec.gov/cgi-bin/browse-edgar?CIK=MSFTFind=Searchowner=excludeaction=getcompany
If you put this web page into your browser you will see the CIK number in
red on the left side of the page near the top.
When I try
75 matches
Mail list logo