On 12/12/2011 02:48 PM, jalfaro wrote:
Hi all.
I am having a problem using color2D.matplot
I am trying to visualize several different matrices under two color ranges.
One color range corresponds to values less than 1 and the second color range
for values greater than 1. However, the minimum
Hi jalfaro,
Apologies for the typo in the last message, the first two lines in the
second example should read:
cellcolors[d1]-
color.scale(c(0,1,d[d1]),cs1=c(0,0,1),cs2=c(0,0,1),cs3=1)[-1:2]
That is, minus the first two values returned.
Jim
__
On 12/12/2011 06:00 PM, Nicola Van Wilgen wrote:
Dear fellow R users,
I would like to draw a sunflowerplot because I have data (decade by
month) that plots multiple times on the same x-y co-ordinates. Further I
would like to colour each of the points/sunflower leaves on the plot
according to
Hello to all,
I still have this doubt.
I'd like to colour the different labels of my dendrogram each one with a
different colour. How can I do? I guess I could do using *edgetext* and
then *t.col* or* lab.col* but I don't know how to add edgetext to my
dendrogram. Can you help me please?
Example:
Hello *,
I'm trying to round time stamps, consisting of a date and time part, to 10
minutes intervals
The timestamps have the following format: %d.%m.%Y %H:%M
*timeStamp - c(31.03.2011 09:30, 31.03.2011 09:41, 31.03.2011 10:04,
31.03.2011 10:10, 31.03.2011 10:28, 31.03.2011 10:35, 31.03.2011
You haven't specified how you want 00:05 rounded (it is neither nor
5). But here's one way
t2 - strptime(timeStamp, format=%d.%m.%Y %H:%M)
t2$min - round(t2$min, -1)
format(t2, %d.%m.%Y %H:%M)
[1] 31.03.2011 09:30 31.03.2011 09:40 31.03.2011 10:00 31.03.2011 10:10
[5] 31.03.2011 10:30
I would like to run a for loop with an index going from 0 to 499 but the
following seems to miss out the first value:
C - 499
for (i in 0:C)
The alternative is:
C - 500
for (i in 1:C)
{
#Then every time I use i, I replace it with i-1
}
Is this a good way to do it or is tere a better way?
Dear list
I have a matrix such as the following:
1 0,5
1 0.7
1 0.5
2 1
3 0
4 0.2
I'd like to plot the histogram of the first column (this is very easy) and
then for each bin of the obtained histogram to plot another histogram on
the z-axis of the frequencies of the second column for the values
On 2011-12-11 22:49, Bert Gunter wrote:
Inline below.
-- Bert
Graph the data sensibly to figure out what's going on. Statistical
machinationsand anova tables with P values alone are not sufficient
and can be opaque or misleading.
If you do not know what sensibly is (or even if you do),
Hi, what about this:
C - 0:499
for (i in C)
{
cat (i, )
}
Best regards,
Andris Jankevics
AIO
Groningen Bioinformatics Centre
Groningen Biomolecular Sciences and Biotechnology Institute
University of Groningen
Kerklaan 30, Haren, 9751 NN, The Netherlands
On Mon, Dec 12, 2011 at 11:44, Thomas
Hi,
On Mon, Dec 12, 2011 at 6:44 AM, Thomas Chesney
thomas.ches...@nottingham.ac.uk wrote:
I would like to run a for loop with an index going from 0 to 499 but the
following seems to miss out the first value:
C - 499
for (i in 0:C)
First off, you've named your variable for an existing
You may want to check how your loop misses out the zero:
C - 5
for (i in 0:C) print(i)
## gives me
[1] 0
[1] 1
[1] 2
[1] 3
[1] 4
[1] 5
Regards,
Enrico
Am 12.12.2011 12:44, schrieb Thomas Chesney:
I would like to run a for loop with an index going from 0 to 499 but the
following seems to
Did you try this:
http://r.789695.n4.nabble.com/coloring-leaves-in-a-hclust-or-dendrogram-plot-td795496.html
or this:
http://bioinformatics.mdanderson.org/Software/OOMPA/ClassDiscovery/html/plotColoredClusters.html
edgetext labels the segments, I think, and is not what I understand
you to want,
On 2011-12-12 1:48, Richard M. Heiberger wrote:
Thank you for you use of HH.
I think the right graph for this data is the much simpler ancova function
library(HH)
ancova(y ~ year * Trt, data=mydata)
where we see that the three treatments have totally different slopes.
Thank you very much for
Hello Everyone
(Quick) question: Does anyone know a package/command or simply a way of
creating a table of different panel data estimations (estimated using
/*plm()*/ ) just as *mtable()* does for models estimated with /*lm()*/?
It seems *mtable* (and *apsrtable* equally) only support /*lm*/ and
On Mon, 12 Dec 2011, ManuelS wrote:
Hello Everyone
(Quick) question: Does anyone know a package/command or simply a way of
creating a table of different panel data estimations (estimated using
/*plm()*/ ) just as *mtable()* does for models estimated with /*lm()*/?
It seems *mtable* (and
Hi!
In my study on kelp I want to see how wave exposure and current (both
continuous variables) affect kelp physical measures along the Norwegian
coast. To assure a balanced design, I stratified my study area into 9
different classes, defined by all combinations of three levels of wave
Hello,
I am trying to write a script to analyse R scripts, which contain only
functions. Therefor, I want to test the R script if it contains onyl correct
code. Is there a nicer and more efficent function than
test - try(source(someRscript))
if(class(test) == try-error){
print(Not executable)
Dear gls-experts,
while reading and testing some examples of the book
introductionary time series analysis with R,
I encountered the following fact which puzzles me.
Confidence intervals for global temperature time series (P99)
computed from general least squares (GLS) to fit the time series.
Dear all,
I'm using the function lm.gls but I'm puzzled by the very limited output. I
do get the regression coefficients but not their standard errors. Neither
summary(), anova() or vcov() works properly for a lm.gls object.
Preferably I look for a function like vcov that will give me the
i use library(ncdf) to read this file as follow
library(ncdf)
sst.nc = open.ncdf(title)
lonall = get.var.ncdf(sst.nc,'lon')
latall = get.var.ncdf(sst.nc,'lat')
precip = get.var.ncdf(sst.nc,'pre')
close(sst.nc)
if i use this method my pc freeze and not respond until i restart it is
there
--
View
Thanks for help. Ive already sorted this out.
robert
--
View this message in context:
http://r.789695.n4.nabble.com/split-date-nad-time-tp4164191p4185788.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing
On Dec 11, 2011, at 10:19 PM, jalfaro wrote:
Hello,
Were you ever able to find a solution to your problem?
I have been trying on and off to solve this problem for several
months now.
I would love to hear if you found a solution.
Neither the person who posted the problem a year ago or you
Dear,
I am using the solnp command (package Rsolnp) for a problem with
equality and inequality constraints.
I am not getting convergence for my problem but apart that I get 1
warning per iteration saying: ?In cbind(temp, funv) : number of rows
of result is not a multiple of vector length
Dear statistician experts,
Sorry if this is a trivial question, or the old same question (i don't know
what is the efficient key word for this issue).
In order to understand the calculation of parameter of logistic regression,
I did an exercise through spreadsheet following the procedural
Hi Jim,
I think this will work well for what I need it for - thanks a lot!
I would still be interested (for interest sake) if anyone knows how to
get this to work using sunflower plot.
Thanks again,
Nicola
-Original Message-
From: Jim Lemon [mailto:j...@bitwrit.com.au]
Sent: 12
Dear Mr. Ripley,
your answer perfectly fits my needs, thank you very much.
(And yes, rounding 5 to even is what I wanted).
Prof Brian Ripley wrote
t2 - strptime(timeStamp, format=%d.%m.%Y %H:%M)
t2$min - round(t2$min, -1)
format(t2, %d.%m.%Y %H:%M)
[1] 31.03.2011 09:30 31.03.2011 09:40
Dear R users, I got the following problem. Given that
data[3,2]
[1] 010252
Code:
intro - data.frame()
intro[1,1] - as.character(data[3,2])
write.csv(intro, file='intro.csv')
In 'intro.csv' file I am loosing the 0 in frot of 10252, which I need. Is
there a way to keep the full character saved?
PEL wrote on 12/07/2011 02:37:42 PM:
Hi all,
I have dataframe that was created from the fusion of two dataframes.
Both
spanned over the same time intervall but contained different
information.
When I put them together, the info overlapped since there is no holes in
the
time interval of
Dear R users,
I am using the code below to generate a fitted value of b. I have about 300
different values for for y (y1, y2, ...y300) which means I will have to write
the code below 300 times to generate the 300 different fitted values for y. Is
there a short way of doing that ?
Many thanks
?prop.test
The very first line in the help file on the function prop.test():
prop.test can be used for testing the null that the proportions
(probabilities of success) in several groups are the same ...
How to interpret the results of the prop.test() you ran? The proportion
preferring
This isn't a reproducible example, so I can't provide code,
but I would put all the y variables in a list and use lapply().
Sarah
On Mon, Dec 12, 2011 at 8:16 AM, Mintewab Bezabih
mintewab.beza...@economics.gu.se wrote:
Dear R users,
I am using the code below to generate a fitted value of b.
On 12/12/2011 01:16 PM, Mintewab Bezabih wrote:
Dear R users,
I am using the code below to generate a fitted value of b. I have about 300
different values for for y (y1, y2, ...y300) which means I will have to write
the code below 300 times to generate the 300 different fitted values for
On 12/12/2011 10:39 AM, Alexander wrote:
Hello,
I am trying to write a script to analyse R scripts, which contain only
functions. Therefor, I want to test the R script if it contains onyl correct
code. Is there a nicer and more efficent function than
test - try(source(someRscript))
threshold wrote on 12/12/2011 06:23:33 AM:
Dear R users, I got the following problem. Given that
data[3,2]
[1] 010252
Code:
intro - data.frame()
intro[1,1] - as.character(data[3,2])
write.csv(intro, file='intro.csv')
In 'intro.csv' file I am loosing the 0 in frot of 10252, which I
On 12/12/2011 10:39 AM, Alexander wrote:
Hello,
I am trying to write a script to analyse R scripts, which contain only
functions. Therefor, I want to test the R script if it contains onyl correct
code. Is there a nicer and more efficent function than
test - try(source(someRscript))
On 12/12/2011 10:19 AM, tony333 wrote:
i use library(ncdf) to read this file as follow
library(ncdf)
sst.nc = open.ncdf(title)
lonall = get.var.ncdf(sst.nc,'lon')
latall = get.var.ncdf(sst.nc,'lat')
precip = get.var.ncdf(sst.nc,'pre')
close(sst.nc)
if i use this method my pc freeze and
Hi
Dear R users,
I am using the code below to generate a fitted value of b. I have about
300 different values for for y (y1, y2, ...y300) which means I will have
to write the code below 300 times to generate the 300 different fitted
values for y. Is there a short way of doing that ?
Thank you! I'm going to look into that straight away.
Manu
--
View this message in context:
http://r.789695.n4.nabble.com/Package-command-for-creating-a-table-of-panel-models-tp4185744p4186077.html
Sent from the R help mailing list archive at Nabble.com.
I was wondering if there is way to place constraints upon the plinear
algorithm of nls, or rather is there a manner in which this can be achieved
because nls does not allow this to be done.
I only want to place constraints on one of the nonlinear parameters, a, such
that it is between 0 and
On Dec 12, 2011, at 7:23 AM, threshold wrote:
Dear R users, I got the following problem. Given that
data[3,2]
[1] 010252
Code:
intro - data.frame()
intro[1,1] - as.character(data[3,2])
write.csv(intro, file='intro.csv')
In 'intro.csv' file I am loosing the 0 in frot of 10252, which I
This question is perhaps best asked on the specialized Bioconductor mailing
list. The audience there is quite knowledgeable and will certainly be able to
help.
Michael
On Dec 11, 2011, at 8:08 AM, ali_protocol mohammadianalimohammad...@gmail.com
wrote:
Dear all,
Is there anyway too
Why can't you just replace temp/data with filename (no quotes)? I'm not sure
I get the question...
Michael
On Dec 12, 2011, at 5:26 AM, Paul Menzel paulepan...@users.sourceforge.net
wrote:
Dear R folks,
I have several data sets I want to process automatically using R. I
found littler
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained,
Thanks a lot. I didn't find that. I prefered the option with class
discovery package which we have to install from OOMPA project.
barbara
On 12 December 2011 11:58, Sarah Goslee sarah.gos...@gmail.com wrote:
Did you try this:
Please sign yourself out, using the linke that appears at the bottom of
each and every message from the list.
On Mon, Dec 12, 2011 at 9:48 AM, Kaigang Li likaig...@yahoo.com wrote:
[[alternative HTML version deleted]]
__
R-help@r-project.org
Dear R users,
I want to test for serial correlation in each equation of VAR.
The multivariate version of this test, for a vector X of class VAR for
instance is
var.ser-serial.test(X,lags.pt=16,type=PT.adjusted) ### in the vars
package
Would any of you suggest me how to modify this to test for
On Mon, Dec 12, 2011 at 2:19 AM, tony333 tony33...@hotmail.com wrote:
i use library(ncdf) to read this file as follow
library(ncdf)
sst.nc = open.ncdf(title)
lonall = get.var.ncdf(sst.nc,'lon')
latall = get.var.ncdf(sst.nc,'lat')
precip = get.var.ncdf(sst.nc,'pre')
close(sst.nc)
if i use
Hello R-Experts,
I've got a question, concerning the automation of a number of regressions
(lm) with the help of a loop (for i in ).
The situation is as follows (the code follows after that):
I have my data in an access database. I have historical data for 2000 parts,
for each of this parts
Hi everybody,
I have a matrix with 3 columns (Date, MeanArea and SdArea). I want to draw a
figure showing the variable MeanArea in terms of the Date. But instead to
use the variable SdArea as bar error, I want to use “polygon error”. I use
this code but the output does not seem good.
Hi Michael,
As far as I know, there is no Arima function in R that allows multiple
seasonal patterns. There's 2 approaches you can take:
1. Use a matrix of dummy variables (e.g. jan/feb/mar/...dec and
mon/tues/wed/...sun) as an xreg term in the arima function.
2. Use a matrix of Fourrier terms
Indeed in txt it looks fine. Anyway, I must stay without 0 because csv is THE
format.
I got another question. why for
Table - matrix(0,8,3)
day = Monday
Table[1,1]=day
all other elements become characters too?
Thanks, robert
--
View this message in context:
The dput() command is essential in doing so.
Michael
On Dec 11, 2011, at 8:50 AM, Tal Galili tal.gal...@gmail.com wrote:
Yes -
E-mail us reproducible code...
Contact
Details:---
Contact me: tal.gal...@gmail.com |
Hi Jim,
Thanks so much for your help. I have read several of your responses on this
mailing list and they have helped me out quite a bit as I have gotten more
and more used to R.
I am still a little confused here by your response.
I think you understood my requirements correctly.
In your
Matrix (which is secretly a vector) can only have one mode
(numeric/factor/character/etc.) for all its elements. If you need
multiple types, go to a data frame
Michael
On Mon, Dec 12, 2011 at 10:39 AM, threshold r.kozar...@gmail.com wrote:
Indeed in txt it looks fine. Anyway, I must stay
Am Montag, den 12.12.2011, 09:41 -0500 schrieb R. Michael Weylandt :
Why can't you just replace temp/data with filename (no quotes)? I'm not
sure I get the question...
… no I am feeling silly and I do not know where I screwed up before.
Here is the full working example.
#!/usr/bin/env
Dear Paul and Sarah
Thanks for the suggestion. I have provided my data here in to make the results
reproducable. I am actually trying to do interpoliation of climate data where
x1 and x2 are my latitude and longitude and sum64-sum 368 are my rainfall
observations which I need to regress
That's not a reproducible example. Paul suggested a list of formulas,
but I recommended creating a list of y variables.
In your attempt, you didn't include the y in the name of the dependent
variable; that's probably why it doesn't work.
Look at this:
y.list - list(y1=runif(10), y2 -
Ah. I did a Silly Thing, and thank you for all the responses helping me track
it down.
The issue lay not in the for (i in 1:C) but in the line VectorName[i] - ...
That's the one that didn't like 0.
Thomas
From: Enrico Schumann [enricoschum...@yahoo.de]
how can I keep the date info after doing unlist to a list of dates?
I have a list of dates were observations were made for each station in
each month
list.obs[[station]][month]
[1] 1979-01-01 1979-01-10 1979-01-25
[1] 0 1 2 3
when i try to unlist i loose the date info.
what
There are methods of the c function for things of class POSIXlt and Date
so do.call(c, dataList) works instead of unlist:
dateList - list(LastWeekend=as.POSIXlt(sprintf(2011-12-%d, 10:11)),
+ Today=as.POSIXlt(2011-12-12))
z - do.call(c, dateList)
z
LastWeekend1
At 15:09 11/12/2011, you wrote:
Dear Michael,
Thanks for the email. This is the kind of forest plot, I want to replicate
http://www.biomedcentral.com/1471-2334/11/122/figure/F2
It would be helpful to cc to r-help in case someone else there knows better.
If you follow my suggestion you get a
Thanks!
Another question also related to the same type of files:
How can I transform it into a time series
as.ts()
I also endup having the same problem, the dates are transformed into a
list of numbers
My file looks like this
datesval
12001-01-121.2
22001-02-12
I don't know why you had c in a(1:2000))
c is a function see ?c ...
and you want a (the row number in SQL_Code) to change with each iteration
in the loop.
Perhaps this might work (I'm not saying this is the best option, just a
potential fix for what you have):
for (a in 1:2000) {
Dataset -
Hi Ryan,
My name is Frederic Andrieu and I work for Mango Solutions in the UK. It looks
like you have swapped the a and c in this for loop statement:
This for (c in a(1:2000)) { should be this for (a in c(1:2000)) {
Which is why the for loop is not running, I would suggest something like this:
Hello,
I'm currently trying to convert a slow and ugly script I made, so that it's
faster and can be computed on a computer grid with the multicore package.
My problem is that I don't see how to turn some loops into an apply-able
function.
Here's an example of my loops :
I got a list of
I am doing a logistic regression, and by accident I included a field
which has the 2digit abbreviation for all 50 states labeled st. I was
surprised to see that the glm did not come up with an error message but
instead appears to have automatically broken down this field into
individual fields
Dear Paul and Sarah
Thanks for the suggestion. I have provided my data here in to make the results
reproducable. I am actually trying to do interpoliation of climate data where
x1 and x2 are my latitude and longitude and sum64-sum 368 are my rainfall
observations which I need to regress
right,
Table - data.frame(matrix(0,8,3))
day = Monday
Table[1,1]=day
Table[1,2]=3
works, thanks a lot. robert
--
View this message in context:
http://r.789695.n4.nabble.com/not-complete-character-in-csv-file-tp4185785p4186639.html
Sent from the R help mailing list archive at Nabble.com.
Kevin,
I'm sorry, but my code has a bug
The correction is:
1) delete the two lines with 'ix'
2) replace them with
for(i in 2:n)
if(abs(x[i] - x[i-1]) delta) x[i] - x[i-1] + delta
That's it. The problem is that it's back to slowness.
Sorry, once again,
Rui Barradas
Hello,
I am dealing with data stored in a database as a 'time' object. I
export the data from the database to a text file and utilize the
'time_to_sec()' function of the database to convert the human readable
time (HH:MM:SS) to seconds so that I can use R to do analysis and
create charts of the
Kevin,
Your problem seems to have three restrictions: (1) abs(x[i] - x[i-1]) =
delta (2) y[i] = x[i]
and (3) minimize sum(y-x). If this is the case I believe I have a better
solution, with a smaller sum
and in much less time.
The problem is restriction (2). If the diffs are negative you can't
...
On Mon, Dec 12, 2011 at 7:00 AM, Jinsong Zhao jsz...@yeah.net wrote:
Apart from the data set here, is there a way to do multiple comparison on
the interaction of one way analysis of covariance?
Interaction in ancova usually means the slopes are different.
...
BTW, when loading HH
--begin inclusion --
I have a matched-case control dataset that I'm using conditional
logistic regression (clogit in survival) to analyze. I'm trying to
conduct k-folds cross validation on my top models but all of the
packages I can find (CVbinary in DAAG, KVX) won't work with clogit
models. Is
Hello everyone,
I want make a variable selection from a dataframe, but when I build this new
object (using cbind) the new table (which is a matrix) I lost all the original
factor names in the variables. I get 1,2,3
Someone would be so kind and tell me if there's is a ny way to get
Without a reproducible example it's impossible to say for certain,
but I'd try
cbind.data.frame() instead of cbind().
You need to have a data frame, not a matrix, for the result.
Sarah
On Mon, Dec 12, 2011 at 2:44 PM, Jose Bustos Melo jbustosm...@yahoo.es wrote:
Hello everyone,
I want make a
Dear Jose,
Am Montag, den 12.12.2011, 19:44 + schrieb Jose Bustos Melo:
Hello everyone,
I want make a variable selection from a dataframe, but when I build this new
object (using cbind) the new table (which is a matrix) I lost all the
original factor names in the variables. I get
Hello,
I am dealing with data stored in a database as a 'time' object. I
export the data from the database to a text file and utilize the
'time_to_sec()' function of the database to convert the human readable
time (HH:MM:SS) to seconds so that I can use R to do analysis and
create charts of the
I really like your varclus function Frank. I have a question about the
quartimax rotation of the 1st two eigen vectors as the SAS VARCLUS algorithm
dictates. According to the SAS user guide example of VARCLUS the only thing
needed to perform the cluster analysis is the variable correlation
I'm glad that worked for you. It's a good idea to reply to the
list, as well as me, so that there's public acknowledgment of
a working solution, and so that the answer appears in the
archives.
Sarah
On Mon, Dec 12, 2011 at 2:56 PM, Jose Bustos Melo jbustosm...@yahoo.es wrote:
Thank you
I would recommend using data.frame(var1, var2, ...)
and not cbind.data.frame(var1, var2, ...). I consider
it bad form to directly call a method of a generic function.
Sometimes it leads to errors, as a method should be
free to assume that its inputs are of the class it
was declared to accept.
Please post only once!
Uwe Ligges
On 12.12.2011 20:32, Tony Stocker wrote:
Hello,
I am dealing with data stored in a database as a 'time' object. I
export the data from the database to a text file and utilize the
'time_to_sec()' function of the database to convert the human readable
time
On 12.12.2011 17:44, Tony Stocker wrote:
Hello,
I am dealing with data stored in a database as a 'time' object. I
export the data from the database to a text file and utilize the
'time_to_sec()'
I get
time_to_sec
Error: object 'time_to_sec' not found
If it is in a package, please tell
On 12.12.2011 19:36, Brian Jensvold wrote:
I am doing a logistic regression, and by accident I included a field
which has the 2digit abbreviation for all 50 states labeled st. I was
surprised to see that the glm did not come up with an error message but
instead appears to have automatically
1. The formula you used is not for a logistic but an ordinal regression
(since you are using the default gaussian family rather than
family=binomial or whatever.
2. R (nor any other software) can deal with perfect separation (nor
quasi-separation) of classes, since the problem is not well
On Dec 12, 2011, at 3:38 PM, Uwe Ligges wrote:
On 12.12.2011 19:36, Brian Jensvold wrote:
I am doing a logistic regression, and by accident I included a field
which has the 2digit abbreviation for all 50 states labeled st.
I was
surprised to see that the glm did not come up with an error
On Dec 12, 2011, at 3:51 PM, Uwe Ligges wrote:
1. The formula you used is not for a logistic but an ordinal
regression (since you are using the default gaussian family rather
than family=binomial or whatever.
this this then produce one version of the Armitage linear test of
trend?
Greetings!
I'm having an issue with the windrose produced by the windrose
function from the circular package. For our weather stations in North
Carolina I'm helping with a script which takes hourly wind speed and
direction data to create windroses for our end users. One of the
stations in the
Hello,
I'm attempting to write a code that automatically imports data from
different files (with different lengths-just one variable) and makes tidy
box plots for comparison. I can successfully import the data and create a
list of the vectors I want to compare. But I cannot, for the life of me,
Hello R folks,
Iâm looking to produce a confusion matrix using the command:
cmx(PData,threshold=0.49,na.rm=FALSE) for a logistic regression model (1 is
habitat, 0 is nonhabitat).
To compile my PData data file (csv file) (as in the example above), Iâm using
the output of the predict(glm)
Hi,
On Mon, Dec 12, 2011 at 4:24 PM, Ryan Utz utz.r...@gmail.com wrote:
Hello,
I'm attempting to write a code that automatically imports data from
different files (with different lengths-just one variable) and makes tidy
box plots for comparison. I can successfully import the data and create
Ryan,
I think you could do what you want by having the vector data written to
separate files; then create a file containing the individual file names. In
R, read the file containing the list of file names and loop through this
reading in the individual vector files. Maybe this is an inelegant,
On 13/12/11 04:39, threshold wrote:
Indeed in txt it looks fine. Anyway, I must stay without 0 because csv is THE
format.
It may well be the case that csv is THE format, but csv is not the problem.
Excel is the problem.
The solution is: Don't use Excel!!! (This is excellent advice in any
Sorry -- previous versiuon prematurely sent. Full version is:
Yikes! You should never have to do this sort of thing (writing stuff
out to files, etc.)
What is wanted, I believe, is ?do.call as in
do.call(boxplot, z)
where z is list(a,b,c) as Sarah described.
However, I think you
Bert,
I may be mistaken, but I thought Ryan wrote write code that
automatically *imports
data* from different files (with different lengths-just one variable), so,
I was referring to doing something with the data before it gets into R. I
understand that one should not need to write out data and
Dear R users,
I know that it was considered many times but having searched through dozens
of posts I have only added to my confusion. My question is - is it possible
to correctly analyse simple mixed (fully)crossed (factorial) 2-way ANOVA?
Assume we have factor A (fixed) and B (random). In the
Ryan,
Do you necessarily have to use list? Have you tried the usecols=TRUE option
in boxplot?
Cheers,
--Kelly V.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Ryan Utz
Sent: Monday, December 12, 2011 1:24 PM
To:
I'm looking for a palette that would maximize the discrepancy with
neighbouring elements.
I have a raster image of objects tagged with integers. I'd like to
highlight the different objects by using a palette that would maximize this
difference with neighbours as close numbers tend to be
Thanks for the help, everyone. However, no solution yet...
My vectors names are neatly stored in a character vector of their own. This
happens based on how they were imported. So, say I have vectors a, b and c,
all with numeric-only contents. Then I have another vector, say z, that
looks like
On 2011-12-12 15:14, Ryan Utz wrote:
Thanks for the help, everyone. However, no solution yet...
My vectors names are neatly stored in a character vector of their own. This
happens based on how they were imported. So, say I have vectors a, b and c,
all with numeric-only contents. Then I have
1 - 100 of 122 matches
Mail list logo
