Hi there
I have an Rscript and I am looking for a way to install a package
non-interactively. In Rscript {Utils}, I saw an example which does
something like this, however this does not seem to work for my particular
example. I am trying to install the following package in an Rscript
(without
Hi,
I am trying to find a simple way to numerically solve a system of two
equations equal to zero with two unknowns (x_loc and y_loc). Here is a mock
data set and below it, the equations I need to solve.
theta-c(180,135,90)/(2*pi)
x-c(0,0,15)
y-c(20,0,0)
0 =
On 10/09/12 01:52, Fred wrote:
Thanks Arun,
I can manage something with that, just need then to delete the first raw
with photoshop !
Huh? WTF?
If you're doing that then you're doing something completely wrong
and not understanding how R works.
You didn't include the context of your
Hi,
I have two matrices. (same size, same number of rows,
same names for the rows - one has data in sorted
order, another one master data in unsorted order. I need to
combine both).
Rows names of Matrix A
P, Q, R, S, T, U, V ...(some order) - each row has 4 values
Row names of
Hi
On 07/09/12 09:35, Benjamin Tyner wrote:
Update: seems one way to skin this cat is to add
gp = gpar(fontfamily=mono)
to the viewport() call itself. If anyone has any suggestions for a
robust way to extract this piece information from the key grob itself
(it's nested several levels
On 12-09-10 2:03 AM, Diviya Smith wrote:
Hi there
I have an Rscript and I am looking for a way to install a package
non-interactively. In Rscript {Utils}, I saw an example which does
something like this, however this does not seem to work for my particular
example. I am trying to install the
Dear fellow R users.
I have released version 1.2.8 of the psych package. This replaces 1.2.4 from
the April release.
The more important changes
o Modified score.irt so that subjects who miss all items or pass all
items are given an estimate based upon the (product) of the
On 10-09-2012, at 01:22, berg1546 wrote:
Hi,
I am trying to find a simple way to numerically solve a system of two
equations equal to zero with two unknowns (x_loc and y_loc). Here is a mock
data set and below it, the equations I need to solve.
theta-c(180,135,90)/(2*pi)
x-c(0,0,15)
Le dimanche 09 septembre 2012 à 19:29 -0500, Sridhar Iyer a écrit :
Hi,
I have two matrices. (same size, same number of rows,
same names for the rows - one has data in sorted
order, another one master data in unsorted order. I need to
combine both).
Rows names of Matrix A
P, Q, R, S,
Sorry to bother everyone. What has changed in Sweave since 2009?
Sweave(Crs.Rnw)
Error: c('Crs.HdrFct.Rnw', 'Crs.Rnw') is not ASCII and does not
declare an encoding
No error in 2009 with the same files.
heers
Christian
--
Christian W. Hoffmann,
CH - 8915 Hausen am Albis, Switzerland
Dear all,
I am including in a plot 6 different lines (?lines) with 6 different line types.
The problem is that I have so dense information that the line types are not
visible any more.
In the code below
myLength-length(currentSet)
plot(seq(from=1,to=myMax,length.out=myLength),currentSet,
Thanks everyone for your comments and help.
Actually, I found my mistake (probably stupid for y'all). I was drawing
first the heatmap and then putting the
png(heatmap.2.png).
What you have to do is :
a- read.csv ( .csv)
a_matrix - data.matrix (a)
row.names(a_matrix) -a$Name # In my case I
Hi,
i would like to plot a few hundred .csv files.
Each file contains one curve with x,y values to plot.
I have been searching for gnu r read many files qplot
and similar words. I found for loops that use assign to generate
one variable containing a dataframe.
When i uesed the classic plot'
Hi
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Fred
Sent: Sunday, September 09, 2012 7:45 PM
To: r-help@r-project.org
Subject: Re: [R] how to save a heatmap.2 in png /jpeg /tiff
Hey Jeff, sorry for the it just doesn't
As usual, Michael was faster than I in responding. Let me add a few thoughts of
my own. See comments below in text.
Best,
Wolfgang
--
Wolfgang Viechtbauer, Ph.D., Statistician
Department of Psychiatry and Psychology
School for Mental Health and Neuroscience
Faculty of Health,
On 09/10/2012 06:14 PM, Alaios wrote:
Dear all,
I am including in a plot 6 different lines (?lines) with 6 different line types.
The problem is that I have so dense information that the line types are not
visible any more.
In the code below
myLength-length(currentSet)
Dear all,
I am including in a plot 6 different lines (?lines) with 6
different line types.
The problem is that I have so dense information that the
line types are not visible any more.
Perhaps plot a lowess or loess curve to show trend instead of a point to point
line?
S Ellison
Hello,
Something like this?
step - myLength %/% 500
ix500 - seq(1, myLength, by = step)
curSubset - currentSet[ix500]
plot(seq(from = 1, to = myMax, length.out = 500), curSubset, ...etc...
Hope this helps,
Rui Barradas
Em 10-09-2012 09:14, Alaios escreveu:
Dear all,
I am including in a
On 12-09-10 3:28 AM, Christian Hoffmann wrote:
Sorry to bother everyone. What has changed in Sweave since 2009?
Sweave(Crs.Rnw)
Error: c('Crs.HdrFct.Rnw', 'Crs.Rnw') is not ASCII and does not
declare an encoding
No error in 2009 with the same files.
heers
Christian
vignettes are now
On 10/09/2012 11:49, Duncan Murdoch wrote:
On 12-09-10 3:28 AM, Christian Hoffmann wrote:
Sorry to bother everyone. What has changed in Sweave since 2009?
Sweave(Crs.Rnw)
Error: c('Crs.HdrFct.Rnw', 'Crs.Rnw') is not ASCII and does not
declare an encoding
No error in 2009 with the same files.
Hi everyone,
Here is a summary of what I am trying to do. Please take a look and let me
know how can i complete this below mentioned objective using R for a large
dataset.
**
- *Objective:*
- To calculate the number of points that fall inside each cell of a
given raster
-
Hi Rolf,
I suggested that in matrix, duplicate row names are possible.
set.seed(1)
dat1-data.frame(a=sample(1:50,5,replace=TRUE),b=rnorm(5,15))
dat2-as.matrix(dat1)
row.names(dat2)-rep(a,5)
A.K.
- Original Message -
From: Rolf Turner rolf.tur...@xtra.co.nz
To: Fred
HI,
You can also use ?merge()
Same example data from Milan:
A - matrix(runif(4*10), 10)
rownames(A) - LETTERS[1:10]
B - matrix(runif(10), 10)
rownames(B) - sample(LETTERS[1:10])
#Either, use
C-as.matrix(merge(A,B[order(rownames(B)),],by=row.names))
#or
HI Sridhar,
Try this:
#I saved the file as .csv.
testMaster-read.csv(test1Master.csv)
test2Sorted-read.csv(test2Sorted.csv)
testMaster1-testMaster[order(testMaster[,1]),]
test2Sorted1-test2Sorted[order(test2Sorted[,1]),]
Combinedat-as.matrix(merge(testMaster1,test2Sorted1,by=X))
Jonas Stein news at jonasstein.de writes:
Hi,
i would like to plot a few hundred .csv files.
Each file contains one curve with x,y values to plot.
I have been searching for gnu r read many files qplot
and similar words. I found for loops that use assign to generate
one variable
Hello,
I´m really sorry to bother you with this problem but it took me a whole day
searching for a solution.
I want the y axis labels positioned closer to the y axis, by default they are
partly outside the window.
Example:
x - as.zoo(EuStockMarkets)
par(las=1)
plot.zoo(x, type=l, main=,
Dear all,
As David mentioned, I used his R-code to try to see the dimension behind the
'LO2' variable. These are the results:
lapply (subsets, function (x) {table(x$LO2)})
[[1]]
nee geen atrofie ja atrofie aanweizg
173 0
[[2]]
nee geen atrofie ja
HI Sridhar,
You may not need to order the dataset by the X to do the merge.
It's already ordered.
Combinedat1-merge(testMaster1,test2Sorted1,by=X)
new1-merge(testMaster,test2Sorted,by=X)
identical(Combinedat1,new1)
#[1] TRUE
A.K.
- Original Message -
From: Sridhar Iyer
Hi ,
I am new to R . I am facing difficulty how to make pairwise comparisions.
For example. I have a file which looks like below
a b c d
x 3 6 7 6
y 7 8 6 5
z 5 4 7 8
Here I need to look for the each pairwise comparisions (ab,ac,ad,bc,bd,cd
for each row)
For instance ,looking at first row,
Thanks for this. A few things
First, yes, my lmer syntax was indeed bad - I was writing this as an example of
what my data code look like. Apologies. So, 1|studyID indeed. Also
1/variance.
I've also been wondering - I often have more than 2 drugs - so, T1, T2, T3,
etc. It depends on
Howdy,
Last week I got some great help on why I was getting an error code when trying
to run this model, thanks everyone! I was able to get the code up and running
beautifully for several data sets. Now I am getting different errors with this
new data set. I can't figure out why, I have
Thank you for your suggestion - my only question is what is the 'par'
variable?
--
View this message in context:
http://r.789695.n4.nabble.com/Solving-a-system-of-two-equations-tp4642626p4642657.html
Sent from the R help mailing list archive at Nabble.com.
Hello, I have a question regarding the predict command for ARIMA models. Why do
I not have to give any input except the number of steps for the forecast? Which
input value is used (in my case for y(t))? Will the last parameter of the
estimate vector be used (in my example 6)? x -
Hello,
I have a list with gene names, fold changes (=expression level) and
chromosomes.
Names fold change chromosome
hz 1.5 2
If I plot fold change versus chromosome (or vice versa):
plot (ch, fc)
I see only the chromosomes with numbers but not those with letter (x
I having issues with my memory size in R for Windows. I have implemented the
3gb switch for Windows 7 32 bit with 8gb of memory, only 4gb is available
for use due to 32 bit limitations. 512Mb is taken by the shared video card,
the memory limit is 2559mb yet I cannot use more than 1.5 gb of
Hi to all
I hope you can help me.
# I'm trying to plot the following ratio data with standard error bars
(horizontal).
ratioBiomass - c(1.327, 0.865, 1.900, 0.992, 1.469, 1.381, 1.230,1.269,
2.411, 1.288, 1.861, 0.714, 1.341, 1.362, 1.065, 2.374)
SEratio - c(19.28, 5.04, 0.01, 0.01, 0.90,
1. Be forewarned: There may be special approaches built into specific
packages that are better than generalist approaches like the one I
suggest here.
2. The general approach to this sort of thing is with ?outer, if I
understand you correctly.
3. Have you read an Introduction to R (ships with
On 10-09-2012, at 17:10, berg1546 wrote:
Thank you for your suggestion - my only question is what is the 'par'
variable?
par is the argument of the function f.
So when you do
a - c(2,3)
f(a)
then in the body of the function par will be equal to c(2,3).
I strongly advise you to read
On 10/09/2012 16:13, chrisk wrote:
I having issues with my memory size in R for Windows. I have implemented the
3gb switch for Windows 7 32 bit with 8gb of memory, only 4gb is available
for use due to 32 bit limitations. 512Mb is taken by the shared video card,
the memory limit is 2559mb yet I
You will have to **show** what you did; it isn't sufficient to
**describe** what you did.
The problem could be that somewhere in your program there is something
that only works in an interactive session.
The problem could be that your command-line attempt is using a different
working directory
The answer to your question probably lies in your data. If you follow the
instructions at the footer of every email on this list and provide an example
that we can run (reproducible), we will be able to do more than guess.
On Sep 10, 2012, at 9:14 AM, Hermann Norpois wrote:
Hello,
I have a list with gene names, fold changes (=expression level) and
chromosomes.
Names fold change chromosome
hz 1.5 2
If I plot fold change versus chromosome (or vice versa):
plot (ch, fc)
On Sep 10, 2012, at 7:00 AM, priya wrote:
Hi ,
I am new to R .
How new? Is this still in a file? ... or has it been imported to an R object.
If it really still is in a file, then Bert's suggestion of first reading
Intro to R rather than posting to R help is on point. If it is an R object,
On Sep 10, 2012, at 11:03 AM, David Winsemius wrote:
On Sep 10, 2012, at 7:00 AM, priya wrote:
Hi ,
I am new to R .
How new? Is this still in a file? ... or has it been imported to an R object.
If it really still is in a file, then Bert's suggestion of first reading
Intro to R
HI,
Sorry, there was a mistake. I forgot to add some code.
Try this:
dat1-read.table(text=
a b c d
x 3 6 7 6
y 7 8 6 5
z 5 4 7 8
,sep=,header=TRUE)
new1-sapply(combn(dat1[,1:4],2),`[`,1:3)
colnames(new1)-sapply(combn(colnames(dat1),2),`[`,1)
new2-as.matrix(new1)
dim(new2)-c(3,2,6)
Hi,
I am trying to identify duplicate values in a column in a date frame. The
duplicated function identifies the duplicate rows in the data frame but it
only does this for the second record, not both records. Is there a way to
mark both rows in the data frame as TRUE?
try this:
dfA$dups-duplicated(dfA$Value) | duplicated(dfA$Value, fromLast = TRUE)
On Mon, Sep 10, 2012 at 2:23 PM, kborgmann borgm...@email.arizona.edu wrote:
Hi,
I am trying to identify duplicate values in a column in a date frame. The
duplicated function identifies the duplicate rows in
On Mon, Sep 10, 2012 at 11:23:42AM -0700, kborgmann wrote:
Hi,
I am trying to identify duplicate values in a column in a date frame. The
duplicated function identifies the duplicate rows in the data frame but it
only does this for the second record, not both records. Is there a way to
mark
A conditioning factor for bwplot() is the stream name. There are 17
streams so each panel is comparatively small and the names in the panel
header are truncated at both ends. I would like to wrap the names on two
lines for each panel but do not see if this can be done when I look at the
online
Hello,
Try the following.
#-- this is new
colfun - colorRampPalette(c(black, white))
cols - colfun(nrow(dfInOutBiom2))
#-- end
plotBiomRatio2 - segplot(
reorder(factor(Species), ratioBiomass) ~ (-SEratio ) + SEratio,
data = dfInOutBiom2,
#-- this
Hi,
What could the cause of
history()
Error in savehistory(file) : no history available to save
savehistory(file=myhist)
Error in savehistory(file) : no history available to save
save.image()
be?.
Cheers
Christian
--
Christian W. Hoffmann,
CH - 8915 Hausen am Albis, Switzerland
On Sep 10, 2012, at 11:35 AM, Rich Shepard wrote:
A conditioning factor for bwplot() is the stream name. There are 17
streams so each panel is comparatively small and the names in the panel
header are truncated at both ends. I would like to wrap the names on two
lines for each panel but do
On Mon, 10 Sep 2012, David Winsemius wrote:
The text that appears in the panel header is called the strip. I
didn't see it in the help page for bwplot. You need to look at:
David,
That explains why I did not find it in the book's index.
?strip.custom.
And since you have the Lattice book
Hello,
Please use dput(dfA) to post your data examples. The following is it's
output. All one has to do is to copy and paste to an R session to get
the data example.
dfA -
structure(list(Site = c(929L, 929L, 930L, 930L, 931L, 932L),
State = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = VA,
Thank you for the help and fast response. I'm going to look into getting
some time on the unix machines.
On Sep 10, 2012 10:28 AM, Prof Brian Ripley [via R]
ml-node+s789695n4642670...@n4.nabble.com wrote:
On 10/09/2012 16:13, chrisk wrote:
I having issues with my memory size in R for
Dear R users,
I want to apply zero inflated models with continuous and categorical
variables and I used pscl package from R and the zeroinf() function. My
question are the follow:
a) The value of fitted.values is mu or (1-p)*mu? where p is the probability
of zero came form a zero point mass
b)
Thanks!!
That did the trick!!
--
View this message in context:
http://r.789695.n4.nabble.com/Identifying-duplicate-rows-tp4642679p4642683.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Dear all,
please do you have any recommendation about a more advanced function in R
for lowess/loess regression ?
the basics lowess() or loess() do not perform as well as I would expect.
thanks very much,
Bogdan
[[alternative HTML version deleted]]
In addition, try this:
x - rep(1:10,2); y - runif(20);f-factor(rep(paste(This
is\nLevel,1:2),e=10))
xyplot(y~x|f,
par.strip.text=list(lineheight=.8,lines=3))
I find the documentation here confusing and incomplete:
par.strip.text is listed as a parameter both in ?xyplot and
?strip.custom. The
On Sep 10, 2012, at 1:43 PM, Bogdan Tanasa wrote:
Dear all,
please do you have any recommendation about a more advanced function in R
for lowess/loess regression ?
the basics lowess() or loess() do not perform as well as I would expect.
thanks very much,
You should probably look at
Thanks, David. I am working on a biology-related project where I have a
matrix of gene-to-gene interactions,
( interactions * distances ) and depending on the basic R functions, I do
get vary different results ...
Bogdan
On Mon, Sep 10, 2012 at 1:52 PM, David Winsemius
On Sep 10, 2012, at 1:48 PM, Bert Gunter wrote:
In addition, try this:
x - rep(1:10,2); y - runif(20);f-factor(rep(paste(This
is\nLevel,1:2),e=10))
xyplot(y~x|f,
par.strip.text=list(lineheight=.8,lines=3))
Bert;
I had imagined that the factor level character strings could be
On Sep 10, 2012, at 8:24 AM, Nicolai Schneider wrote:
Hello,
I´m really sorry to bother you with this problem but it took me a whole day
searching for a solution.
I want the y axis labels positioned closer to the y axis, by default they are
partly outside the window.
Example:
Hi,
I just got started analyzing data for a charter school, and am trying to
open Riverside Publishing's file of ITBS scores. I just assume that
someone, somewhere has opened this file before, and I was wondering if this
person had the code right now. There are over 100 variables, all of
Dear R helpers,
I am trying to test for linkage disequilibrium (LD) in R between all pairs
of loci. Importantly, my loci are microsatellites and have up to 30
alleles per locus. Do you know of any R packages that can either (1)
calculate measures of LD directly or (2) test for the presence of
On Mon, 10 Sep 2012, Bert Gunter wrote:
In addition, try this:
x - rep(1:10,2); y - runif(20);f-factor(rep(paste(This
is\nLevel,1:2),e=10))
xyplot(y~x|f,
par.strip.text=list(lineheight=.8,lines=3))
Bert,
Thanks very much. If the docs are less than clear to you, they will
certainly be
Either way, there is a par() parameter to control the axis and label title
position.
Look at the entry for 'mgp' in ?par
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 9/10/12 2:17 PM, David Winsemius
A permissions problem in the current working directory?
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 9/10/12 1:29 PM, Christian Hoffmann c-w.hoffm...@sunrise.ch wrote:
Hi,
What could the cause of
history()
Error in
Hello R List,
I'm trying to run a linear mixed effects analysis on fairly large datasets with
lmer (from the lme4 package) on a 32-bit Windows machine running XP with 3 GB
of RAM. It's not working. (details below)
I've researched the ff and bigmemory packages, but it appears they won't handle
Is it somehow possible to use the the predict method of the raster package
with a prediction model obtained using caret's train()?
/
predict(rasterStack, carettrain, filename=..., progress='text',
format=BSQ,
datatype='INT1U', type='response', overwrite=TRUE) /
Thank you
--
View this
To exclude the color key use
plotBiomRatio2 - segplot(
reorder(factor(Species), ratioBiomass) ~ (-SEratio ) + SEratio,
data = dfInOutBiom2,
col = black,
level = NULL,
#-- this changed
colorkey = FALSE,
#-- end
draw.bands = FALSE,
Your problem seems to be exactly the problem I've been having for the past
week. So glad to find this thread.
However, when I follow the syntax you used to get around the issue, I still
didn't get the output to work in the way yours did. I don't believe my
factors are truly nested however (nor
On Sep 10, 2012, at 5:14 PM, semperparatus wrote:
Your problem seems to be exactly the problem I've been having for the past
week. So glad to find this thread.
However, when I follow the syntax you used to get around the issue, I still
didn't get the output to work in the way yours did. I
Hi all
I've been getting some emails from the R package maintainers that I need
to update some code in a CRAN packge that uses FORTRAN, to comply with
(not so recent) changes. I've been a little busy...
I'm having trouble adjusting my code. I hope someone can help.
The package was working
McCall, Ken (CMG-Dayton Ken.McCall at coxinc.com writes:
I'm trying to run a linear mixed effects analysis on fairly large
datasets with lmer (from the lme4 package) on a 32-bit Windows
machine running XP with 3 GB of RAM. It's not working. (details
below)
I've researched the ff and
I want to change it because I don't want to compare in this instance between
conditions, but I simply want to see the contrast t-statistic between
patient and control at every level of condition (1, 2, and 3).
From there I'd like to be able to plot the t-statistic for the contrast
between patient
I don't understand at all what the problem is, but a workaround
is simply not to specify the pdf argument in name=value form
but rather in positional form. I.e. just make it the first argument
to .Fortran(). As in:
tmp - .Fortran(pdf, other parameters here)
Works for me! :-)
cheers,
Hi Team,
I am not sure how to run some r queries in continuos time interval for say in
each 30 minutes.
Will be very grateful if some one can give me some heads up in good direction.
Unfortunately I dont have any query written or data which I can share and I
know most of the respected people
78 matches
Mail list logo