Hi Andrew,
Do you know which coefficients are missing (i.e. which terms have
missing coefficients)? This would help to narrow down the possibilities.
Is it certain that all levels of all factors occur in all replicates? If
not then you will get different numbers of coefficients (in the
parame
Hi!
I want to create a random matrix with 15 variables, each variable having
1000 observations.
Between each two variables, I want to define a specific (*not *random)
correlations between them, but still saving the "randomness" of each
variable (mean=zero, s.d=1).
How can I do this in R?
thanks,
Thank you all - these suggestions have been very helpful! I've got it doing
what I wanted know, and I appreciate the help!
Emily
On Jan 27, 2013, at 1:18 PM, Rui Barradas wrote:
> Hello,
>
> Try the following. Create a file called test.R with these instructions:
>
> cmd <- commandArgs(TRUE)
Hi,
I am intending to save a path-describing character object in a slot of a
class I'm working on. In order to have the option to use "system.file" etc
in such string-saved path definitions, I wrote this
ExpressionEvaluator <- function(x){
x <- tryCatch(
expr=base::parse(text=x),
erro
Thanks a million for all help provided!! I can do what I intend to using
the "for loop". However, I'm still eager to try the list.files approach.
Here is the error message that I got using Ivan's code:
> list_of_dataset <- do.call(read.table, file_names)
Error in do.call(read.table, file_names) :
On 28-01-2013, at 06:17, Purna chander wrote:
> Dear admin members,
>
>
> My Inbox is being flooded with the posts every time. As a reason, I
> wish to unsubscribe from this forum.
>
> Can you suggest me how to do that.
>
At the bottom of every message is a link to the R-help mailing list p
Dear List,
I'm using gam in a multiple imputation framework -- specifying the knot
locations, and saving the results of multiple models, each of which is
fit with slightly different data (because some of it is predicted when
missing). In MI, coefficients from multiple models are averaged, as a
Dear admin members,
My Inbox is being flooded with the posts every time. As a reason, I
wish to unsubscribe from this forum.
Can you suggest me how to do that.
Regards,
Purna
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/
I have only *two* datasets from normal and cancer samples.
CancerNormal
--
mRNA1 3049
mRNA2 199200
... ... ...
mRNA1000 1340
Each samples contain several thousan mRNA microarray express
On Jan 27, 2013, at 8:16 PM, Giuseppe Amatulli wrote:
Dear all,
thanks for your input.
Bert - yes you get the point, i would to like to draw an ellipse
contour
for a population quantile.
Indeed, as you mention data.ellipse() should draw that. In other
words if i
re-run my model for anothe
Hi Karl! Thanks for writing!
Doesn't this format require a column for a factor for every variable
present in any observation, whether or not that variable is present in the
observation in question? I think I end up with data that consists mainly of
columns of variables that are NAs for all but a f
Dear all,
thanks for your input.
Bert - yes you get the point, i would to like to draw an ellipse contour
for a population quantile.
Indeed, as you mention data.ellipse() should draw that. In other words if i
re-run my model for another prediction (getting a new vector b) i would
have the 95% prob
Hi
On 17/01/13 13:19, p_conno...@slingshot.co.nz wrote:
Paul Murell's article "What's in a Name" in The R Journal Vol 4/2
gives an interesting example of editing a stacked barplot of the barley
data. Using the method described in that article, it's easy to do
something along the lines of
grid.
Hi Fabrice,
The cross.area parameter gives the size of the intersection, which cannot be
larger than the size of either set 1 (area1 parameter) or set 2 (area2
parameter). You probably want:
venn.plot <- draw.pairwise.venn(
area1 = 3186 + 5880,
area2 = 325 + 5880,
cross.area = 5880);
Paul
--
How do you get those values from the example header file that you included?
On Sun, Jan 27, 2013 at 9:15 PM, Pascal Oettli wrote:
> ?extract
>
> HTH
> Pascal
>
>
> Le 27/01/2013 22:53, Jonsson a écrit :
>>
>> having 12 files with 12 hdrs for one year:these files are raster
>> (projected WGS8
?extract
HTH
Pascal
Le 27/01/2013 22:53, Jonsson a écrit :
having 12 files with 12 hdrs for one year:these files are raster
(projected WGS84,lat
long):https://echange-fichiers.inra.fr/get?k=rLSyoavrnifGyH5XrlO
samples = 1440
lines = 720
ba
Depends what algorithm you are using. Read ?set.seed
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Live:
Ok. Thanks!
On Sun, Jan 27, 2013 at 7:44 PM, Jeff Newmiller
wrote:
> ?cat
> ---
> Jeff NewmillerThe . . Go Live...
> DCN:Basics: ##.#. ##.#. Live Go...
>
?cat
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (
Dear R People:
Here is an unusual question: if you are using "source", it just
prints the output. If you use source with "echo", you get both
commands and output. Is there a way just to show the commands,
please?
Thanks,
Erin
--
Erin Hodgess
Associate Professor
Department of Computer and M
Thank you all.
Yes. I just miss understand this part. Now it is OK.
On Sun, Jan 27, 2013 at 7:41 PM, Paul Boutros wrote:
> Hi Fabrice,
>
> The cross.area parameter gives the size of the intersection, which cannot be
> larger than the size of either set 1 (area1 parameter) or set 2 (area2
> par
On 27-Jan-2013 23:50:57 Ben Bolker wrote:
> Fabrice Tourre gmail.com> writes:
>
>> Dear list,
>> When I use VennDiagram package, I got a error as follow:
>>
>> venn.plot <- draw.pairwise.venn(
>> area1 = 3186,
>> area2 = 325,
>> cross.area = 5880);
>>
>> Error in draw.pairwise.venn(area1 = 318
Fabrice Tourre gmail.com> writes:
>
> Dear list,
>
> When I use VennDiagram package, I got a error as follow:
>
> venn.plot <- draw.pairwise.venn(
> area1 = 3186,
> area2 = 325,
> cross.area = 5880);
>
> Error in draw.pairwise.venn(area1 = 3186, area2 = 325, cross.area = 588) :
> Impossible
All:
Aha! -- The light dawneth, methinks (maybe...)
Giuseppe: I re-read the SAS link you sent and if I have parsed it
correctly, what SAS chooses to call an ellipse for "prediction" -- a
rather idiosyncratic way to describe it, imo -- I believe the rest of
us would call a contour for a population
Hi,
I am learning R. I've been using set.seed() for a while, but without
actually understanding the significance of the "number" we put in the
brackets. e.g. set.seed(135) & set.seed(930).
Can anyone shed some light on this please?
--
View this message in context:
http://r.789695.n4.nabble.co
Thanks for your help! Your original solution would have worked just fine
too if I didn't have a number of other large data frames loaded. This is
my first time working with such large data sets, so I'd never previously
run out of available memory. I finished up the project, though I'm sure
all b
I am trying to download multiple files from a ftp. However, I get this error
after downloading several files: cannot open URL ''. The stopage occurs
randomly; can be after the 10th or 30th etc. The code I am using is:
#
#rm(list=ls(all=TRUE))
require(RCurl)
require("R.u
Hi,
I did some googling. I found multiple ODE solver packages. But I am
wondering if there is any terminal value ODE solver (linear case) for
multi-dimension case or not? I have not come across any. Any guidance would
be appreciated.
Robert
[[alternative HTML version deleted]]
_
Dear list,
When I use VennDiagram package, I got a error as follow:
venn.plot <- draw.pairwise.venn(
area1 = 3186,
area2 = 325,
cross.area = 5880);
Error in draw.pairwise.venn(area1 = 3186, area2 = 325, cross.area = 588) :
Impossible: cross section area too large.
Does anyone have suggestion
On Jan 27, 2013, at 4:10 PM, Dennis Fisher wrote:
R 2.15.1
OS X
Colleagues,
I encountered the following unexpected behavior today:
The following command yielded the expected result:
paste(c("TEXT1", "TEXT2"), collapse="|")
Result:
[1] "TEXT1|TEXT2"
However, abbreviating "col
Hi Dennis: One of function argument matching rules in R that arguments
after the dotdotdot have to be matched exactly ( see code for paste below )
so that's why your attempt doesn't work. But I would have been surprised
also so I'm not trying to imply that one should know which functions have
arg
R 2.15.1
OS X
Colleagues,
I encountered the following unexpected behavior today:
The following command yielded the expected result:
paste(c("TEXT1", "TEXT2"), collapse="|")
Result:
[1] "TEXT1|TEXT2"
However, abbreviating "collapse" by even one character:
paste(c("TEXT1",
Hi Michael!
Yes, I do use Cohens d. As a matter of fact my thesis supervisor told me to use
1 as the value for standard deviation for all of my studies.
Unfortunately I am not totally sure myself why to do this have you ever used
such an approach?
kind regards,
Alma
Dear Giuseppe and Bert,
I also didn't follow what's intended, more or less for the same reasons as
Bert mentioned, which is why I didn't reply to the initial posting. In the
car package, confidenceEllipse() draws confidence ellipses for a pair of
coefficients from a statistical model, and dataElli
Hello,
Try the following. Create a file called test.R with these instructions:
cmd <- commandArgs(TRUE)
if(any(grepl("--ext", cmd))){
suffix <- cmd[grep("--ext", cmd)]
suffix <- unlist(strsplit(suffix, ":"))[2]
fl <- paste0("test_", suffix, ".txt") # underscore included
}
Hi,
I tried with bigger dataset.
set.seed(25)
names <- sample(c("bob", "joe", "cr...@gmail.com", "emily",
"j...@yahoo.com"),5e6,replace=TRUE)
set.seed(1651)
emails
<- sample(c("b...@cup.com", "joesm...@gmail.com", "cr...@gmail.com",
"emi...@yahoo.com", "j...@yahoo.com"),5e6,replace=TRUE)
df
Thanks all. I will give them all a go and let you know the outcome.
kind regards
--
View this message in context:
http://r.789695.n4.nabble.com/Converting-column-of-strings-to-boolean-tp4656739p4656774.html
Sent from the R help mailing list archive at Nabble.com.
_
Hi,
It sounds like you just want to write a command line script using R
and you would pass the suffix/prefix as command line args, no?
Why not just go with what Peter has already suggested with
`commandArgs`, or if you want a more feature-rich command line arg
parser, you can try:
http://cran.r-
Hi Yasha,
I guess you got Uwe's response.
I created `df2` with the intention of getting the two results from the
original dataset.
For example, after you get the first result
df[,1][grep("@",df$names)]<- ""
#you can get the second result by:
df[df$names!="",]
# names emails
#1
Hi,
thanks for your replay.
My values of a and b are respectively:
a = observation of an event
b = prediction of a model.
Therefore i would like to draw the confidence region for predicting a new
observation, and according to this
http://v8doc.sas.com/sashtml/insight/chap40/sect35.htm the predicti
Hello all (again),
I received a very helpful answer to this question, and would like to pose one
more:
Right now I have this script, which is being called from the command line,
writing output to two generically named files ("pvalues" and "qvalues") that
are named in the script using the lin
Hi,
When I look at the summary of an rpart object run on my data, I get 7 nodes but
when I plot the rpart object, I get only 3 nodes. Should the number of nodes
not match in the results of the 2 functions (summary and plot) or it is not
always the same?
Look forward to your reply,
Carol
--
You two were 100% right, it was just a memory issue. This was part of a
bigger project where I had a number of data frames loaded, all with 1-5
million rows. Cleaned up my code to have less data frames loaded at once,
and everything is working great. Thanks for the help!
On Jan 27, 2013 9:46 AM,
You appear to be quite confused. I have no idea what your "a" and "b"
below mean. The SAS documentation that you quote is for prediction
from a bivariate normal. I don't know what that has to do with your
problem, nor with car's confidence ellipses for parameters for a
univariate regression. I woul
Hi,
You could use library(plyr) as well
library(plyr)
pnew<-colSums(aaply(laply(split(as.data.frame(p),((1:nrow(as.data.frame(p))-1)%/%
25)+1),as.matrix),c(2,3),function(x) x))
res<-rbind(t(pnew),colSums(p))
row.names(res)<-1:nrow(res)
res<- 100-100*abs(res/rowSums(res)-(1/3))
A.K.
- Orig
Hi all, I have a set of 54 files that I need to convert from ASCII grid
format to .shp files to .bnd files for BayesX.
I have the following R code to operate on those files:
library(maptools)
library(Grid2Polygons)
library(BayesX)
library(BayesXsrc)
library(R2BayesX)
readfunct <- function(x)
{
u
On 13-01-27 10:37 AM, Alexander Senger wrote:
Hello useRs,
I would like to draw a 3D-surface using rgl with a point-like
light-source within the scene, that is with finite distance of the
light-source to the surface to be lit.
The rgl package doesn't support that.
From the help to the 'li
Hello useRs,
I would like to draw a 3D-surface using rgl with a point-like
light-source within the scene, that is with finite distance of the
light-source to the surface to be lit.
>From the help to the 'light3d' command I read:
"They [the light-sources] are positioned either in world space or
Hi Ray!
I'm insisting with list.files...!
What about like this (untested)?
file_names <- list.files(path="C:/.../data", pattern=".dat$",
full.names=TRUE)
list_of_dataset <- do.call(read.table, file_names)
Let me know if this helps!
Ivan
--
Ivan CALANDRA
Université de Bourgogne
UMR CNRS/uB 6282
On 27.01.2013 12:50, Richard D. Morey wrote:
Dear Contributors,
I am asking help on the way how to solve a problem related to loops for
that I always get confused with.
I would like to perform the following procedure in a compact way.
Consider that p is a matrix composed of 100 rows and three
On 27.01.2013 07:11, ypodeswa wrote:
Actually, it worked perfectly for my sample data, but my actual data has
5.5 million rows, and grep doesn't seem to work with over a million rows.
Any idea on a workaround?
It is not a matter of grep() but of available memory, I guess.
Hence try to redu
having 12 files with 12 hdrs for one year:these files are raster
(projected WGS84,lat
long):https://echange-fichiers.inra.fr/get?k=rLSyoavrnifGyH5XrlO
samples = 1440
lines = 720
bands = 1
header offset = 0
file type
> Dear Contributors,
> I am asking help on the way how to solve a problem related to loops for
> that I always get confused with.
> I would like to perform the following procedure in a compact way.
>
> Consider that p is a matrix composed of 100 rows and three columns. I need
> to calculate the su
On Jan 27, 2013, at 08:33 , Emily Sessa wrote:
> Hi all,
>
> I am trying to use the scan function in an R script that I am calling from
> the command line on a Mac; at the shell prompt I type:
>
> $ Rscript get_q_values.R LRT_codeml_output
>
> in the hope that LRT_codeml_output will get pa
Hello,
Something like this?
x <- sample(c("red", "blue", "green", "yellow"), 100, replace = TRUE)
cnames <- unique(x)
sapply(cnames, function(.x) x == .x)
Hope this helps,
Rui Barradas
Em 26-01-2013 22:25, domcastro escreveu:
Hi
I'm trying to convert a column of strings (nominal types) to
Hello,
I think there is an error in the expression
100-(100*abs(fa1[i]/sum(fa1[i])-(1/3)))
Note that fa1[i]/sum(fa1[i]) is always 1. If it's fa1[i]/sum(fa1), try
the following, using lists to hold the results.
# Make up some data
set.seed(6628)
p <- matrix(runif(300), nrow = 100)
idx <- se
That said,
> wilcox_test(x ~ factor(y), distribution = "exact")
or the same with oneway_test, i.e would be ok?
2013/1/27 Achim Zeileis
> On Sun, 27 Jan 2013, Kay Cichini wrote:
>
> Thanks for the reply!
>>
>> Still, aren't there issues with 2-sample test vs y and excess zeroes
>> (->many tie
On Sun, 27 Jan 2013, Kay Cichini wrote:
Thanks for the reply!
Still, aren't there issues with 2-sample test vs y and excess zeroes
(->many ties), like for Mann-Whitney-U tests?
If you use the (approximate) exact distribution, that is no problem.
The problem with the Wilcoxon/Mann-Whitney tes
Thanks for the reply!
Still, aren't there issues with 2-sample test vs y and excess zeroes
(->many ties), like for Mann-Whitney-U tests?
Kind regards,
Kay
2013/1/26 Achim Zeileis
> On Fri, 25 Jan 2013, Kay Cichini wrote:
>
> Hello,
>>
>> I'm searching for a test that applies to a dataset (N=
Hi all,
I am trying to use the scan function in an R script that I am calling from the
command line on a Mac; at the shell prompt I type:
$ Rscript get_q_values.R LRT_codeml_output
in the hope that LRT_codeml_output will get passed to the get_q_values R
script. The first line of that script
Hi,
Anther possibility may be to use:
library(MatrixModels)
model.Matrix(~d-1,sparse=FALSE)
#7 x 4 Matrix of class "ddenseModelMatrix"
# dblue dgreen dred dyellow
#1 0 0 1 0
#2 0 1 0 0
#3 0 0 1 0
#4 1 0 0 0
#5 0
Hi
I'm trying to convert a column of strings (nominal types) to a set of
boolean / binary / logical values. For example, in the column there is red,
blue, green and yellow. There are 100 rows and each has a colour. I want to
convert the column to 4 columns: red, blue, green,yellow and then either
Thank you _very much_ Ilai for the rapid and accurate answer.
It works and indeed helps a lot. Both to solve the question and to help
me progress !
Possibly this will also help others.
Thanks again
Tito
Le 27/01/2013 01:11, ilai a écrit :
>
>
> On Sat, Jan 26, 2013 at 10:26 AM, Tito de Morais
Dear Contributors,
I am asking help on the way how to solve a problem related to loops for
that I always get confused with.
I would like to perform the following procedure in a compact way.
Consider that p is a matrix composed of 100 rows and three columns. I need
to calculate the sum over some ro
64 matches
Mail list logo