Hi Marcos,
I'd be glad to help, but your example is not usable for anyone but
yourself - I don't have the files that your code reads and I don't have
any idea what the resulting data frames look like. Please provide
examples of the data and the result.
Andreas
Marcos Amaris Gonzalez schrieb
Dear Renoir,
are you referring to:
http://econ.la.psu.edu/~hbierens/TVCOINT.PDF
?
If so, no, but you could up this framework fairly easily and hereby employ the
functions of urca. But this should already be evident from the package's
manual.
Best,
Bernhard
> -Ursprüngliche Nac
Dear list,
I have the model below which i am using to account for spatial autocorrelation:
exponential <-corExp(form = ~ Longitude + Latitude)
explanation_mod_all <-
gls(Lower_PD~Area+Elevation+Temperature+Preceipitation+Agriculture+Urban+Human.footprint+Population,
correlation = exponential)
Dear R users,
pROC is a package to compare, visualize, and smooth receiver operating
characteristic (ROC) curves.
The package provides the following features:
* Partial or full area under the curve (AUC) computation
* Comparison of two ROC curves (curves and AUC)
* Calculating and plotting confid
On Thu, Mar 31, 2011 at 9:04 PM, Shi, Tao wrote:
> This question has been asked by many people already. The easiest way is:
>
> 1) install the new version
> 2) copy all or the libraries that you installed later from the "library"
> folder
> of older version to the new version
> 3) uninstall the
Dear R Users,
I use for my data crunching a combination of MySQL and GNU R. I have
to handle huge/ middle seized data which is stored in a MySql
database, R executes a SQL command to fetch the data and does the
plotting with the build in R plotting functions.
The (low level) calculations like sum
On Fri, Apr 1, 2011 at 1:46 PM, Henri Mone wrote:
> Dear R Users,
>
> I use for my data crunching a combination of MySQL and GNU R. I have
> to handle huge/ middle seized data which is stored in a MySql
> database, R executes a SQL command to fetch the data and does the
> plotting with the build i
Hi there,
Could someone help me with the following programming problem..?
I have written a function that works for my intended purpose, but it
is quite closely tied to a particular dataframe and the names of the
variables in this dataframe. However, I'd like to use the same
function for different
On Fri, 1 Apr 2011, Henri Mone wrote:
Dear R Users,
I use for my data crunching a combination of MySQL and GNU R. I have
to handle huge/ middle seized data which is stored in a MySql
database, R executes a SQL command to fetch the data and does the
plotting with the build in R plotting function
On Fri, Apr 1, 2011 at 6:46 AM, Henri Mone wrote:
> Dear R Users,
>
> I use for my data crunching a combination of MySQL and GNU R. I have
> to handle huge/ middle seized data which is stored in a MySql
> database, R executes a SQL command to fetch the data and does the
> plotting with the build i
On Fri, Apr 1, 2011 at 11:46 AM, Henri Mone wrote:
> Dear R Users,
>
> I use for my data crunching a combination of MySQL and GNU R. I have
> to handle huge/ middle seized data which is stored in a MySql
> database, R executes a SQL command to fetch the data and does the
> plotting with the build
See the warning in ?subset.
Passing the column name of lvar is not the same as passing the 'contextual
column' (as I coin it in these circumstances).
You can solve it by indeed using [] instead.
For my own comfort, here is the relevant line from your original function:
Data.tmp <- subset(Fulldf, l
> The R results I'm getting are similar to the SAS results but not very
> close. For example, my coefficient for TRT is .2913 in R but SAS
> gives .3450. I'd like to be able to run the R Code with method =
> "exact" to make it as same as the AS example but I can't seem to get
> it to work
You n
On Thu, Mar 31, 2011 at 3:26 AM, Marius Hofert wrote:
> Dear Deepayan,
>
> thanks for answering. It's never too late to be useful.
>
> I see your point in the minimal example. I checked the z-axis limits in my
> original problem for the point to be inside and it wasn't there. I can't
> easily
> r
Thanks Nick and Juan for your replies.
Nick, thanks for pointing out the warning in subset(). I'm not sure
though I understand the example you provided -- because despite using
subset() rather than bracket notation, the original function (myfunct)
does what is expected of it. The problem I have is
This should be a version that does what you want.
Because you named the variable lvarname, I assumed you were already passing
"lvar" instead of trying to pass lvar (without the quotes), which is in no
way a 'name'.
myfunct.better <- function(subgroup, lvarname, xvarname, yvarname,
dataframe)
{
Setting Z=Q-A would be the incorrect dimensions. I could Z=Q/A. Is
fitting a nls model the same as fitting an ols? These data are
hydraulic data from ~47 sites. To access predictive ability I am
removing one site fitting a new model and then accessing the fit with
a myriad of model assessment c
Ruben,
One more thing you might try is to add a 'size' component to
the lines list. The default is size=5, but since you have
enought space in the plot, try size=8 or 9; the line types
will show up more clearly:
lines=list(pch=1:4, lty=1:4, type='b', size=8),
Peter Ehlers
useR! 2011 Conference, University of Warwick, Coventry, UK
16-18 August, 2011
(International R user conference; see
http://www.r-project.org/conferences.html
for the history.)
The abstract submission period for this year's useR! conference, which
was due to end today, is *extended* by 2 days.
Dear all,
I am a happy user of R console, but I would like to see syntax
coloring. I use R 2.12 in Ubuntu Linux.
I have found the packages "xterm256" and "highlight", but I was not
able to figure out how to use it to highlight the syntax in console
output.
Also, I tried several GUI interfaces, b
On Fri, Mar 11, 2011 at 12:28 AM, wrote:
> Good afternoon,
>
> I am trying to create a plot where the bottom and top axes have the same
> scale but different tick marks. I tried user-defined xscale.component
> function but it does not produce desired results. Can anybody suggest
> where my use
2011/4/1 Nick Sabbe :
> This should be a version that does what you want.
Indeed it does, thank you very much!
> Because you named the variable lvarname, I assumed you were already passing
> "lvar" instead of trying to pass lvar (without the quotes), which is in no
> way a 'name'.
Sorry about th
Dear R-Users,
Am I right if I say that comparisons among gls models with differing fixed
effects, but same correlation stuctures, fitted by REML is not meaningful? I
need to use ML fit to get comparable AIC, BIC values? And the same is
applied to likelihood ratio statistics?
I don't really know wh
Hello,
I am writing an R-guide which is addressed especially for students of German
Universities
(http://www.rrzn.uni-hannover.de/buch.html?&no_cache=1&titel=statistik_r).
For the title page, I want to use the R logo. I found a similar question in
the R mailing list about the sam
RStudio is a new interface. I just started using it and find it very good.
You might want to check it out. It does have Ubuntu capabilities, but you
have be sure that your graphic drivers are the same as is specified in the
RStudio system.
Hope this helps
Steve Friedman Ph. D.
Ecologist / Spati
Hello
On Fri, Apr 1, 2011 at 12:32 PM, January Weiner
wrote:
> I have found the packages "xterm256" and "highlight", but I was not
> able to figure out how to use it to highlight the syntax in console
> output.
>
I've been in this position before, without finding a solution for
syntax highlighti
On Fri, Mar 11, 2011 at 12:29 PM, Marius Hofert wrote:
> Dear Deepayan,
>
> many thanks for answering.
>
> Another thing I am wondering is the following: I know you can have (3d-like)
> "crosses" in the wireframe plot. But are there any other 3d-like plot
> symbols? Of course one can use differe
Dear all,
I am trying the barplot command but some of the labels are disappearing as
there is not enough place on the graph to put them all.
Here is an example of code that doesn't show all the labels:
barplot(sort(runif(9,0,0.2),decreasing=TRUE),xlim=c(0,20),width=2,names.arg=c("first
nam
Try the cex.names argument in barplot function.
On Fri, Apr 1, 2011 at 10:34 AM, Samuel Le wrote:
> Dear all,
>
>
>
> I am trying the barplot command but some of the labels are disappearing as
> there is not enough place on the graph to put them all.
>
> Here is an example of code that doesn't s
On Sat, Mar 19, 2011 at 2:53 AM, Guy Jett wrote:
> Due to an error on my part, I have renamed the previously attached file from
> T_5-04b_LTC-SE-SO-Compared.csv to
> T_5-04b_LTC-SE-SO-Compared.txt.
> It remains a comma-delimited file so the extension can be changed and used
> per the script, or
On Sat, Mar 19, 2011 at 6:11 AM, Evans, David G (DFG)
wrote:
> Hi,
>
> >From the following code (tweaked from another user):
>
>
>
> variable<-sample(rep(1:2,100))
>
> individual<-rep(1:3, length(variable))
>
> group<-rep(LETTERS[1:2],length(variable)/2)
>
> mydata<-data.frame(variable,individual,
Hello
I have two issues
1. I have a linear model with several thousands of Y and X variables, so if I
do not want to write them down one by one in the lm(Y ~ X) model how should I
go about it.
2. Also, if I suspect that some of the Y variables could not be independent,
what function should I
> Recently RStudio was introduced [1] and, although beta, it received
> quasi-unanimous acclaim from the community, so you risk finding it
> useful too.
Dear Liviu,
RStudio might be a fine program, but it does not feature syntax
highlighting, which is the only thing I am missing from R Console (i
I have an Affybatch object called "batch" :
>
> batch
AffyBatch object
size of arrays=1050x1050 features (196 kb)
cdf=HuGene-1_0-st-v1 (32321 affyids)
number of samples=384
number of genes=32321
annotation=hugene10stv1
notes=
>
>
Is there a way of converting a portion of this data into a matrix?
Hi all,
I have a huge data set. All I want to do is to change the column names of
the it. So if I use read.table, I can read the data in and change colnames
and then write back, such as,
t <- read.table("a.txt", header = T, sep = "\t")
colnames(t)
colnames(t) <- c(") # new column names
wri
Dear January,
Have a look at Eclipse with the STAT-ET plugin http://www.walware.de/goto/statet
Best regards,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek
team Biometrie & Kwaliteitszorg
Gaverstra
Hi,
I am running a test to call an R script with in excel using VBA. My VBA
code is shown bellow. The middle section of this mail also includes the
content of my Rscript. The bottom part shows the error message form the
R console.
It seems that Excel is opening the R console without any proble
Jorge,
This looks like an R on Windows problem, it is unrelated to VBA.
You have the statement
("source('X:\Trading\Energy\JorgeSpace\TMPholder\cpixe\home
models\toto.R')")
which is invalid R. You need to use either forward slashes or
double-back-slashes in R commands.
See the R for Windows F
On Fri, Apr 1, 2011 at 3:48 PM, January Weiner
wrote:
> RStudio might be a fine program, but it does not feature syntax
> highlighting, which is the only thing I am missing from R Console (it
> only colors the commands typed).
>
The idea is that you shouldn't use the R console for you main
program
On 01-04-2011, at 16:12, Jorge Nieves wrote:
>
> Hi,
>
> I am running a test to call an R script with in excel using VBA. My VBA
> code is shown bellow. The middle section of this mail also includes the
> content of my Rscript. The bottom part shows the error message form the
> R console.
>
>
On Fri, Mar 25, 2011 at 3:59 PM, JP wrote:
> Hi there David,
>
> Many thanks for your time and reply
>
> I created a small test set, and ran your proposed solution... and this is
> what I get http://i.imgur.com/vlsSQ.png
> This is not what I want - I want separate grp_1 and grp_2 panels and in eac
Use "/" instead of "\" in the path.
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Software/Embedded Contro
Jorge:
You can run save your scripts in the same folder where your wokbook is and run
it like this:
Sub Scatter()
Call rinterface.StartRServer
'Put the dataframe into R
Call rinterface.PutDataframe("scatter", DownRightFrom(Range("RData!A1")),
WithRowNames:=False)
Call rin
Interesting variety of solutions! Thanks very much.
Murali
-Original Message-
From: Henrique Dallazuanna [mailto:www...@gmail.com]
Sent: 31 March 2011 18:26
To: Menon Murali
Cc: r-help@r-project.org
Subject: Re: [R] choosing best 'match' for given factor
Try this:
bestMatch <- function
On Sun, Mar 27, 2011 at 8:50 PM, Marius Hofert wrote:
> Dear expeRts,
>
> I partially managed to obtain what I wanted by using latticeExtra. However,
> the
> following questions remain:
> 1) why do not all x-axis labels appear? [compare bw and bw2]
You have unnecessarily asked for relation="free
Thanks to Deepayan for sending me the code for the correct approach
using subscripts.
On Mar 13, 2011, at 9:46 PM, David Winsemius wrote:
On Mar 13, 2011, at 8:51 PM, Mark Linderman wrote:
David, thank you for your quick reply. I spent a few minutes
getting your
command to work with some
On Apr 1, 2011, at 10:00 AM, hongsheng wu wrote:
Hi all,
I have a huge data set. All I want to do is to change the column
names of
the it. So if I use read.table, I can read the data in and change
colnames
and then write back, such as,
t <- read.table("a.txt", header = T, sep = "\t")
co
Dear January,
Have you looked at Emacs + ESS? http://ess.r-project.org/
It highlights in the text editor and the actual R process besides
coming with a rich set of features and a mailing list filled with
helpful Emacs & R users. I've tried several different interfaces and
ended up being happies
Dear R-listers,
I would like to have your helps on make a good strategy of quality
control by several quality control variables. This means I need a good
strategy to choose cutoff for each quality control variables, or make
cutoff for all such quality control variables in one time.
For quality co
Am 01.04.2011 17:02, schrieb Joshua Wiley:
Dear January,
Have you looked at Emacs + ESS? http://ess.r-project.org/
It highlights in the text editor and the actual R process besides
coming with a rich set of features and a mailing list filled with
helpful Emacs& R users. I've tried several di
Dear Prof. Murdoch and Dr. Winsemius,
Thank you very much for your suggestion. They were big help.
One of my friends told me another approach by using "smooth.spline". Here is
the code:
x <- c(1,2,3,4,5,6)
y <- c(1,4,9,16,25,36)
sp <- smooth.spline(x,y)
x <- seq(1, 6, length=100)
pred <- predict(
Dear R users,
I'm trying to add a regression line on my boxplots (something
like:boxplot(c(1:3),c(4:6),c(5:8)))
But I can't see it.
Please help !!!
It's not a April fool's joke!!!
Yan
[[alternative HTML version deleted]]
__
R-help@r-project.o
Hello list,
we are beginners in R and we are trying to fit the following time series:
*> x <- c(1.89, 2.46, 3.23, 3.95, 4.56, 5.07, 5.62, 6.16, 6.26,
6.56, 6.98, 7.36, 7.53, 7.84, 8.09)*
The reason for choosing the present time series is that the we
have previously calculated analytically th
Hi Terry and Kevin,
Â
Thanks for your replies to my post. And thanks particularly to Terry for
pointing out my confusion about the "exact partial likelihood" vs. the "exact
marginal likelihood." I'm learning and so am likely to be confused about such
things from time to time.
Â
I had thought
Dear Liviu,
thanks for the programming tip! However, I do all my editing in vim,
which has had syntax highlighting for quite a while, as well as
auto-completion and a number of other goodies. But while I do most of
my programing in vim, I do most of my scientific studies simply using
the R interfa
On Apr 1, 2011, at 11:24 AM, Yan Jiao wrote:
boxplot(c(1:3),c(4:6),c(5:8))
df <- structure(list(grp = c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3),
val = c(1L, 2L, 3L, 4L, 5L, 6L, 5L, 6L, 7L, 8L)),
.Names = c("grp", "val"),
row.names = c(NA, -10L),
class = "data.frame")
boxplot(val~grp,data=df)
Am 01.04.2011 11:24, schrieb Yan Jiao:
Dear R users,
I'm trying to add a regression line on my boxplots (something
like:boxplot(c(1:3),c(4:6),c(5:8))) But I can't see it. Please help
!!! It's not a April fool's joke!!!
Your sample data does not make any sense (at least to me). I would do it
a
Please move this thread off r-help. It's about Rstudio, not R.
-- Bert
On Fri, Apr 1, 2011 at 8:05 AM, January Weiner
wrote:
> Dear Liviu,
>
> thanks for the programming tip! However, I do all my editing in vim,
> which has had syntax highlighting for quite a while, as well as
> auto-completion
On Fri, Apr 1, 2011 at 5:05 PM, January Weiner
wrote:
> Dear Liviu,
>
> thanks for the programming tip! However, I do all my editing in vim,
> which has had syntax highlighting for quite a while, as well as
> auto-completion and a number of other goodies. But while I do most of
> my programing in
See inline.
On 2011-03-31 22:22, Yusuke Fukuda wrote:
Thanks Bert.
I have read "?formula" again and again, and I'm still struggling;
lm(body_length ~ head_length-1)
This removes intercept from each individual regression (for male, female,
unknown).
When they are taken together,
lm(body_
On 2011-04-01 05:44, stephen sefick wrote:
Setting Z=Q-A would be the incorrect dimensions. I could Z=Q/A. Is
fitting a nls model the same as fitting an ols? These data are
hydraulic data from ~47 sites. To access predictive ability I am
removing one site fitting a new model and then accessin
R community:
We're excited to post a preliminary agenda for the upcoming 3rd
conference on R and Applied Finance, to be held in Chicago on April
29th and 30th.
In addition to keynotes from John Bollinger, Mebane Faber, Stefano
Iacus and Louis Kates, we are excited to have 31 additional talks
cove
Dear expeRts,
I know I can't paste expressions in the normal way, but I just couldn't figure
out
how to get the following (I want to paste a character vector to an expression
vector)
right with bquote() or substitute.
vec1 <- c("a", expression(tilde(b)), "c")
vec2 <- c("1", "2", "3")
main <- a
Great. thank you Bernd! Learned a new thing here.
John
From: Bernd Weiss
Cc: r-help@r-project.org
Sent: Thu, March 31, 2011 6:19:25 PM
Subject: Re: [R] regular expression
Am 31.03.2011 21:06, schrieb array chip:
> Ok then this code didn't do what I wanted. I
Try this also:
> grep("arg", "arg.symptom", value = TRUE, invert = TRUE)
character(0)
> grep("arg", "liver.symptom", value = TRUE, invert = TRUE)
[1] "liver.symptom"
>
On Fri, Apr 1, 2011 at 1:33 PM, array chip wrote:
> Great. thank you Bernd! Learned a new thing here.
>
> John
>
>
>
>
> __
I have a matrix with 2 columns and I want to do fishers exact test for these
with the totals for each row being 100 say.
The data has the form:
23 12
32 21
12 2
and these represents the tables:
23 12
77 88
32 21
78 79
12 2
88 98
How do I use apply to speed up aclculation of the fishe
Dear R users,
How to show mean in the boxplot instead of median ?
Many thanks
Yan
[[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://w
On Apr 1, 2011, at 12:32 PM, Marius Hofert wrote:
Dear expeRts,
I know I can't paste expressions in the normal way, but I just
couldn't figure out
how to get the following (I want to paste a character vector to an
expression vector)
right with bquote() or substitute.
vec1 <- c("a", expre
On Apr 1, 2011, at 1:20 PM, Yan Jiao wrote:
Dear R users,
How to show mean in the boxplot instead of median ?
Read ... help(boxplot)
Use boxplot to create a set of values, substitute th emean for the
groups then pass the mangled set of values to bxp().
--
David Winsemius, MD
West Hartf
Dear David,
thanks for your reply. The paste() meant to paste the vectors vec1 and vec2
together, so main should be a vector of length 3 of the form
"a 1", "b 2", "c 3" # with b being tilde(b)
However, with c() it is a vector of length 6:
expression("a", tilde(b), "c", "1", "2", "3")
Do you know
On Apr 1, 2011, at 1:56 PM, Marius Hofert wrote:
Dear David,
thanks for your reply. The paste() meant to paste the vectors vec1
and vec2 together, so main should be a vector of length 3 of the form
"a 1", "b 2", "c 3" # with b being tilde(b)
However, with c() it is a vector of length 6:
exp
On Apr 1, 2011, at 2:34 PM, David Winsemius wrote:
On Apr 1, 2011, at 1:56 PM, Marius Hofert wrote:
Dear David,
thanks for your reply. The paste() meant to paste the vectors vec1
and vec2 together, so main should be a vector of length 3 of the form
"a 1", "b 2", "c 3" # with b being tild
On Apr 1, 2011, at 2:34 PM, David Winsemius wrote:
On Apr 1, 2011, at 1:56 PM, Marius Hofert wrote:
Dear David,
thanks for your reply. The paste() meant to paste the vectors vec1
and vec2 together, so main should be a vector of length 3 of the form
"a 1", "b 2", "c 3" # with b being tild
I'm stumped, can anyone find my error in this sequence.
for(j in 1:(varsize[4]-1))
temp <- get.var.ncdf(nc=input,
varid="p_foraging",c(1,1,j),c(varsize[1],varsize[2],1))
filled.contour(x, y, temp, color = terrain.colors,
plot.title = title(main = paste("Everglades Wo
> -Original Message-
> From: stephen sefick [mailto:ssef...@gmail.com]
> Sent: April-01-11 5:44 AM
> To: Steven McKinney
> Cc: R help
> Subject: Re: [R] Linear Model with curve fitting parameter?
>
> Setting Z=Q-A would be the incorrect dimensions. I could Z=Q/A.
I suspect this is con
On 2011-04-01 10:19, Jim Silverton wrote:
I have a matrix with 2 columns and I want to do fishers exact test for these
with the totals for each row being 100 say.
The data has the form:
23 12
32 21
12 2
and these represents the tables:
23 12
77 88
32 21
78 79
12 2
88 98
How do I us
> On Fri, Mar 11, 2011 at 12:28 AM,
> wrote:
> > Good afternoon,
> >
> > I am trying to create a plot where the bottom and top axes have the
> > same scale but different tick marks. I tried user-defined
> > xscale.component function but it does not produce desired results.
> > Can anybody
On 2011-04-01 10:28, David Winsemius wrote:
On Apr 1, 2011, at 1:20 PM, Yan Jiao wrote:
Dear R users,
How to show mean in the boxplot instead of median ?
Read ... help(boxplot)
Use boxplot to create a set of values, substitute th emean for the
groups then pass the mangled set of values to
Thank you very much Deepayan. This does exactly what I needed!
Cheers,
Guy Jett
gj...@itsi.com
-Original Message-
From: Deepayan Sarkar [mailto:deepayan.sar...@gmail.com]
Sent: Friday, April 01, 2011 6:51 AM
To: Guy Jett
Cc: r-help@R-project.org
Subject: Re: [R] XYPlot Conditioning Varia
Aren't you missing a set of parentheses?
I can't run your code since it's not reproducible, but to
my aging eyes it seems that you need a set of '{}' around
the contents of your loop:
for(j in 1:(varsize[4]-1)) { loop stuff }
Peter Ehlers
On 2011-04-01 11:48, steve_fried...@nps.gov wrote:
I'
Hi Peter,
Thanks for taking the time to consider the problem. I realize the procedure
is not reproducible. I use a 4 dimensional netCDF file (4.2 GB) in size to
pull data into this process. Nobody in their right mind should work with
such things. It's a spatial temporal database with 10 years of d
Hi, I tried to pass the function dist() as an argument, but got an error
message. However, almost the same code with mean() as the function to be
passed,
it works ok.
foo<-function (x,
xfun = dist)
{
xfun(x)
}
foo(matrix(1:100,nrow=5))
Error in foo(matrix(1:100, nrow = 5)) : could not fi
Hi R helpers... I am having troubles because of the discrepancy
between the dendrogram plotted from hclust and what is wrote in the
hc2Newick file. I've got a matrix C:
hc <- hclust(dist(C))
plot(hc)
with the:
write(hc2Newick(hc),file='test.newick')
both things draw completely differen
Hi Peter,
Thanks for taking the time to consider the problem. I realize the procedure
is not reproducible. I use a 4 dimensional netCDF file (4.2 GB) in size to
pull data into this process. Nobody in their right mind should work with
such things. It's a spatial temporal database with 10 years of d
Hi all,
I have made an overdispersion test for a data set and get the following result
Overdispersion test Obs.Var/Theor.Var Statistic p-value
poisson data 16.24267 47444.85 0
after deleting the outliers from the data set I get the following result
Overdispersion test
On Apr 1, 2011, at 5:38 PM, steve_fried...@nps.gov wrote:
Hi Peter,
Thanks for taking the time to consider the problem. I realize the
procedure
is not reproducible. I use a 4 dimensional netCDF file (4.2 GB) in
size to
pull data into this process. Nobody in their right mind should work
w
Steve,
I use filled.contour but have a semicolon, ";", between the two
axis calls. Do you need one after
"axis(1, seq(45 , 58, by = 1))"
?
Clint
--
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.e
On Apr 1, 2011, at 5:56 PM, array chip wrote:
Hi, I tried to pass the function dist() as an argument, but got an
error
message. However, almost the same code with mean() as the function
to be passed,
it works ok.
foo<-function (x,
xfun = dist)
{
xfun(x)
}
foo(matrix(1:100,nrow=5))
Erro
OK, I figured it out, need to add stats::: before dist
foo<-function (x,
xfun = stats:::dist)
{
xfun(x)
}
John
To: r-help@r-project.org
Sent: Fri, April 1, 2011 2:56:06 PM
Subject: [R] function in argument
Hi, I tried to pass the function dist() as an ar
Or:
foo<-function (x, xfun = dist) {
match.fun(xfun)(x)
}
On Fri, Apr 1, 2011 at 7:06 PM, array chip wrote:
> OK, I figured it out, need to add stats::: before dist
>
> foo<-function (x,
> xfun = stats:::dist)
> {
> xfun(x)
> }
>
>
> John
>
>
>
>
> To: r-h
I have a question about wireframe 3-D plots and how to apply colors. I have
a large dataset of river flow (m^3/s) over time, and I have coded these
flows based on their height. I would like to produce a wireframe plot that
colors the graph based on the flow code, i.e. I would like high flows to b
Steve,
I was just about to send off another suggestion which
began with ...
Try the following or wait for David W. to chime in.
But I see that he's already done that. Anyway, in case
it's still of some use:
First I would check that you don't have a variable 'T' in
your workspace. (Never use
Hello everyone!I have a data set with missing observations that I am trying to
impute. I am using MICE and I would like the imputed values to all be positive.
I have two types of variables: prices (P1 to P136) and quantities (Q1 to Q136)
and I also want the range of these two types to be differ
Hi Peter and everyone else who has offered suggestions. I am now at home
and will try your suggestions this evening.
The initial effort with the stripped down verion I have used earlier did
access the file correctly so I have no reason to think it is not formatted
correctly. That however
On Apr 1, 2011, at 6:21 PM, Pam Allen wrote:
I have a question about wireframe 3-D plots and how to apply
colors. I have
a large dataset of river flow (m^3/s) over time, and I have coded
these
flows based on their height. I would like to produce a wireframe
plot that
colors the graph bas
I'm stuck here. The following code works great if I just use it on the
command line in a 'for' loop. Once I try to define it as a function as seen
below, I get bad results. I'm trying to return a vector of column indices
based on a short string that is contained in a selection (length of about
70)
Hi Ben,
I am having some trouble figuring out what it is exactly you want. A
workable example would be nice and then (even if just manually) typed
out what you would like it to return based on a given input. I almost
think you just want grep().
grep("test", c("test", "not", "test2", "not"), val
Hi all!
I have a vector, let's say for example int <- sample(1:20,10);
for now:
now I have a matrix...
M = m x n
where the first column is a "feature" column and most likely shares at least
one of the int (interesting) numbers.
I want to extract the rows where int[] = M[,1]
I thought:
rownames(
Hi:
Here' s one approach:
> int <- sample(1:20,10)
> m <- matrix(sample(1:40, 20), nrow = 10)
> int
[1] 7 12 4 6 1 19 17 20 15 5
> m
[,1] [,2]
[1,]9 15
[2,] 23 32
[3,] 40 14
[4,] 19 38
[5,] 286
[6,] 26 18
[7,] 34 22
[8,]7 35
[9,] 21
HI,
here is another solution:
int <- sample(1:20,10)
int
[1] 10 4 5 2 14 17 9 11 16 13
mat<-matrix(11:30,ncol=4)
mat
[,1] [,2] [,3] [,4]
[1,] 11 16 21 26
[2,] 12 17 22 27
[3,] 13 18 23 28
[4,] 14 19 24 29
[5,] 15 20 25 30
mat[apply(mat,1, funct
100 matches
Mail list logo
