thanks all.your help much appreciated. and its usefull
--
View this message in context:
http://r.789695.n4.nabble.com/HELP-how-to-remove-10-of-data-randomly-in-R-tp4647879p4648098.html
Sent from the R help mailing list archive at Nabble.com.
__
R-hel
HI,
I need to parse an xml where key name are same but values are different.
Some dummy text
one
two
three
When i use xmlGetAttr() function i always get one as value. How can i get
one , two and three? Pls help me out of this pblm .
--
View th
Indeed this seems to be a bug, if the graph is directed. The
workaround is to convert it to an undirected graph, the clique
computation ignores direction anyway:
library(igraph)
adj <- matrix(1, nrow=11, ncol=11) - diag(11)
g <- graph.adjacency(adj)
largest.cliques(g)
# [[1]]
# [1] 10 8 1
#
# [
Hi,
I have a panel data set that I am trying to subset. I am trying to keep
values for years >=1960.
The full set is from 1940 to 2000.
I tried a few things, but none worked. Here are a couple that I am trying to
use.
TableAPS1 <- subset(TableAP, year => 1959)
TableAPS1 <- TableAP[ which(year
There is a problem with the largest clique computation in the recent version
of igraph.
> library(igraph)
> adj <- matrix(1, nrow=11, ncol=11) - diag(11)
> g <- graph.adjacency(adj)
> largest.cliques(g)
[[1]]
[1] 10 8 1
[[2]]
[1] 9 7 1
[[3]]
[1] 8 7 1
Warning message:
In largest.cliques(g) :
Hi
I'd like to fit an asymmetrical curve function to some physiological data. I've
been told a weibull curve is a good place to start, but I'm having trouble
specifying and fitting the function with nls and was wondering if someone could
help.
After some reading, I think the function specificat
WOW Richard. Nice. This really did the trick. I will continue to work
at it to tweak things.
This community is great. Thanks so much all of you. I am making it a
priority to be rid of my noob status.
Kenneth Stephen Dyson
Chercheur postdoctoral
Groupe de recherche sur le système nerveux central
Hi Heiberger,
Thanks for your response. Tolfer had told me the answer, just use print ()
in the xyplot, and it works. Anyway thanks again.
Best regards
Yuanzhen Lin
my data is shown as follow:
> dump("tmp", "")
tmp <-
structure(list(x = c(1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L,
1L, 2L, 3L, 4L,
I made up some data since you didn't send your example.
tmp <- data.frame(x=rep(1:10, 3),
y=rnorm(30),
Fam=rep(c("a","b","c"), each=10),
Spacing=rep(c("A","B"), each=5, times=3))
xyplot(y ~ x | Fam, groups=Spacing, data=tmp, layout=c(3,1), aut
See FAQ 7.22.
cheers,
Rolf Turner
On 01/11/12 12:09, yzh lin wrote:
Hi everybody,
I recently used xyplot, but I found that for control flow can't work in the
xyplot. The program code and some of my data file are shown as follow.
Would someone please tell me where my code is wrong
I would use bwplot in lattice for this example.
library(lattice)
library(latticeExtra)
mydata$dtime <- factor(mydata$dtime)
useOuterStrips(bwplot(BRCLNET ~ dtime | side * group, data=mydata,
xlab="time"))
On Wed, Oct 31, 2012 at 3:45 PM, Rui Barradas wrote:
> Hello,
>
> Like this?
>
>
> library
Hi everybody,
I recently used xyplot, but I found that for control flow can't work in the
xyplot. The program code and some of my data file are shown as follow.
Would someone please tell me where my code is wrong and how can make the
code right?
Thanks.
Best regards
Yuanzhen Lin
Code:
library
I'm having trouble getting ffdfindexget to work right in Windows. Even the
most trivial of examples gives me problems.
> myVec = ff(1:5)
> another = ff(10:14)
> littleFrame = ffdf(myVec, another)
> posVec = ff(c(2, 4), vmode = 'integer')
> ffdfindexget(littleFrame, posVec)
Error in if (any(B < 1
See the CRAN "Robust" Task View (for better ways of doing what you want).
-- Bert
On Wed, Oct 31, 2012 at 2:26 PM, Dieter Best wrote:
> I have distribution that are gaussian to a good approximation. I fit a
> gaussian to these distributons. Once in a while there is an outlier. Could
> someone su
Dear R help,
I am trying to cluster my data according to "group" in a data frame such as
the following:
df=data.frame(group=rep(c("a","b","c","d"),10),(replicate(100,rnorm(40
I'm not sure how to tell hclust() that I want to cluster according to the
group variable. For example:
dfclust=hc
Dear R help,
I am trying to cluster my data according to "group" in a data frame such as
the following:
df=data.frame(group=rep(c("a","b","c","d"),10),(replicate(100,rnorm(40
I'm not sure how to tell hclust() that I want to cluster according to the
group variable. For example:
dfclust=hc
I have distribution that are gaussian to a good approximation. I fit a
gaussian to these distributons. Once in a while there is an outlier. Could
someone suggest a robust method (R package already?) that removes those
outliers and redoes the gaussian fit to get a better fit? Thanks.
[[alte
thank you both! that worked!
--
View this message in context:
http://r.789695.n4.nabble.com/expand-grip-for-permutations-tp4648067p4648076.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat
Okay; I think I get it. One (rather kludgy) way of accomplishing your goal
is as follows: Let the two matrices which you want to check for being
"isomorphic" be A and B. Do:
rA <- apply(A,1,paste,collapse=" ")
rB <- apply(B,1,paste,collapse=" ")
cA <- apply(A,2,paste,collapse=" ")
cB <- apply(
On Wed, 31 Oct 2012, aelmore wrote:
I'm hoping that folks out there with expertise in working with the party
package can help me out here. My team is trying to convert party tree
output into a text file format that can be read by our image processing
software. We are running into difficulties
Sounds like one of your data columns is constant. The variance of a
constant is 0, and scaling would then divide by 0, which is impossible.
Kevin Wright
On Wed, Oct 31, 2012 at 7:47 AM, fillay89 wrote:
> I am trying to run the R Script below, I have actually simplified it to
> just
> this par
Instead of
rep(drugs, block.size)
use
rep(list(drugs), block.size)
as the argument to expand.grid.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf
> Of trekvana
> S
On Wed, Oct 31, 2012 at 8:25 PM, David Winsemius wrote:
>
> On Oct 31, 2012, at 10:46 AM, Hard Core wrote:
>
>> I run R as administrator but when i try to install packages i get this
>> message:
>>
>>> install.packages()
>> Warning in install.packages :
>> apertura non riuscita: stato HTTP '403 F
Compiling this little function gets me some strange behaviour
.initDataDir <- function(){
if(file.exists(LOCATION)) {
if(as.logical(file.info(LOCATION)["isdir"]))return
stop(LOCATION, " exists but is not a directory")
}
Z <- dir.create(LOCATION)
if(!Z){
stop(geterrmessage())
##
Webinar signup:
Tree Ensembles and Classification Trees: A Winning Combination
November 9, 10-11 a.m., PST.
Webinar Registration:
http://2.salford-systems.com/gradientboosting/
Understand major shortcomings of using only decision trees and how tree
ensembles can help overcome these challenges
if i were to have a block size of 4 people and i want to assign a treatment
combination to the entire block, there would be 16 different treatment
combinations (, TTTP, TTPP, PTTP, etc.)
i am trying to get all 16 permutations and i am able to use this code below.
drugs=c('P','T');
comb=expand
On 10/31/2012 02:47 PM, sylvain.giroud wrote:
Dear R-users,
Does somebody know what does the "Estimate" reported by the Lavaan model
tell us?
I assume this tells the relative strength of the dyadic relations.
The 'Estimate' column contains the estimated model parameters. There are
many differe
On Oct 31, 2012, at 10:46 AM, Hard Core wrote:
> I run R as administrator but when i try to install packages i get this
> message:
>
>> install.packages()
> Warning in install.packages :
> apertura non riuscita: stato HTTP '403 Forbidden'
> Warning in install.packages :
> apertura non riuscita
Hello,
Like this?
library(ggplot2)
set.seed(3101)
mydata <- data.frame(
BRCLNET = rnorm(1000),
group = rep(c("1", "2"), each = 500),
side = sample(c("Left", "Middle", "Right"), 1000, replace = TRUE) ,
dtime = rep(1:10, 100))
p <- ggplot(mydata , aes(factor(dtime), BRCL
HI,
May be this helps.
dat1<-read.table(text="
TDate TTime O3 No2 Temp Sun Wspeed Wdirect Hum Indicator
1 19980101 2400 0.065 0.036 31.4 765 9.9 351 NA 1
2 19980102 2400 0.053 0.025 31.8 624 7.7 351 NA 1
3 19980103 2400 0.027 0.033 31.5 852 8.8
This worked for the example you provided. Assumes the header count is the only
numeric value on the 5th line.
epa_extract <- function(address){
doc <- readLines(address, n = 5)[5]
head_count <- as.numeric(gsub("\\D", "", doc))
Hi everyone,
Longtime reader first time questioner. A colleague was showing me the
EViews statistical package and demonstrating some of its regression
capabilities--like DOLS--and commenting how easy it was to do certain kinds
of regression. Being more of an R type, I went to try and replicate his
Dear R friends.
I have a question about running a glm( family= 'binomial', *offset=T*), (I
know offset is a vector of values)
My doubt is about predicting the values on a new data. Does the predict()
function considers the offset? o should I especified something?
Here is the model I´m using:
*mod
I'm hoping that folks out there with expertise in working with the party
package can help me out here. My team is trying to convert party tree
output into a text file format that can be read by our image processing
software. We are running into difficulties because the way the two
different progr
I run R as administrator but when i try to install packages i get this
message:
> install.packages()
Warning in install.packages :
apertura non riuscita: stato HTTP '403 Forbidden'
Warning in install.packages :
apertura non riuscita: stato HTTP '403 Forbidden'
Warning in install.packages :
u
ok this is close, thanks for your effort
i will be more specific about my data set now, maybe that will help.
there are 2 cats, labeled rs29 and rs30
there are a number of experimental results, i will make a graph set for
each one. for example we can look at the results of BRCLNET.
the tasks are
On Oct 31, 2012, at 11:11 AM, David Winsemius wrote:
>
> On Oct 31, 2012, at 9:46 AM, chuck.01 wrote:
>
>> Hello,
>>
>> Here is a link to some data:
>> http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt
>>
>> I am trying to read this in, and want to use:
>> chmva
On Oct 31, 2012, at 9:46 AM, chuck.01 wrote:
> Hello,
>
> Here is a link to some data:
> http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt
>
> I am trying to read this in, and want to use:
> chmval <-
> read.table("http://www.epa.gov/emap/html/data/surfwatr/data/m
Using na.string works better:
> x <-
> read.csv("http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt";,
> skip=84, as.is = TRUE, na.string = '.')
> str(x)
'data.frame': 711 obs. of 75 variables:
$ ALDI: int NA NA NA NA NA NA NA NA NA NA ...
$ ALDS: chr
This worked fine for me:
> x <-
> read.csv("http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt";,
> skip=84, as.is = TRUE)
> str(x)
'data.frame': 711 obs. of 75 variables:
$ ALDI: chr "." "." "." "." ...
$ ALDS: chr "." "S" "S" "S" ...
$ ALDSF : chr
[see in-line below]
On 31-Oct-2012 10:26:14 PIKAL Petr wrote:
> Hi Ted
>
>> -Original Message-
>> From: ted@deb [mailto:ted@deb] On Behalf Of Ted Harding
>> Sent: Tuesday, October 30, 2012 6:41 PM
>> To: r-help@r-project.org
>
>
>
>>
>> The general asymptotic result for the pth quanti
Hello,
Use readLines instead.
?readLines # see argument 'n'
readLines("http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt";,
n = 5)[5]
Hope this helps,
Rui Barradas
Em 31-10-2012 16:46, chuck.01 escreveu:
Hello,
Here is a link to some data:
http://www.epa.gov
Found the problem on this. This is a new computer and does not have Java
installed. Loaded JDK and everything is working. Thanks.
-Original Message-
From: Pascal Oettli [mailto:kri...@ymail.com]
Sent: Wednesday, October 31, 2012 10:15 AM
To: Ashton, Triss
Cc: r-help@r-project.org
Su
Hello,
Here is a link to some data:
http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt
I am trying to read this in, and want to use:
chmval <-
read.table("http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/wchem/chmval.txt";,
sep=",", skip= 84, header=T)
Sorry Sarah.
I want to store them as a vector for use later.
so, similar to this:
links <-
c("http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/benthic/benmet.txt";,
"http://www.epa.gov/emap/html/data/surfwatr/data/mastreams/9396/location/watchr.txt";,
"http://www.epa.gov/emap/htm
Sure.
set.seed(1)
dat1<-data.frame(Height=sample(150:180,12,replace=TRUE),EyeColor=rep(c("Green","Blue","Brown"),each=4))
dat1[,2]=="Blue"
# [1] FALSE FALSE FALSE FALSE TRUE TRUE TRUE TRUE FALSE FALSE FALSE FALSE #
Returns a logical vector which checks whether each row of 2nd column of dat1
Dear Michael ... You're a genius thak you very much
really thank you !!! It worked!!!
--
View this message in context:
http://r.789695.n4.nabble.com/Mean-Help-tp4648000p4648038.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@
Hello,
I'm not sure if this is what you want.
aggregate(delayValue ~ id + Freq, data = sampleData, length)
Hope this helps,
Rui Barradas
Em 31-10-2012 16:38, Edward Patzelt escreveu:
> R-help -
>
> I have this set of aggregated tables (sample data below via dput()). And I
> would like to hav
> I use kinship package of R and it doesn't create a kinship matrix with
repeated id.
The kinship package is out of date, use the kinship2 and coxme packages (it was split into
two parts). Then the problem you describe no longer exists -- the kinship function no
longer requires each subject i
On Oct 31, 2012, at 11:38 AM, Edward Patzelt wrote:
> R-help -
>
> I have this set of aggregated tables (sample data below via dput()). And I
> would like to have delayValue as the column variables with the "temp"
> (temp1, temp2, temp3) values as the row variables. However I would like to
>
On Oct 31, 2012, at 9:56 AM, chuck.01 wrote:
> Sorry, I know I should read a little 1st about this, but I am actually just
> helping somebody really quick and need help too.
>
> I want to grep all of the names of the .txt files mentioned on this html web
> page:
>
> http://www.epa.gov/emap/rem
I am working with the splm package. I use the spgm function: general
estimation of a panel data model. Based on this approach, I know it is
possible to compute a R2, eg the ratio of variation explained by a given
model.
My model is :
bivmod<-spgm(logIKA~NBLITRE0+NBLITRE1,data=mydatap,listw=coms
they're all of the form
http.*txt
but the best way to "grep" them (by which I assume you mean extract the
file names from the page source) depends on what you plan to do with them,
and what sort of output you expect.
It isn't even clear whether you plan to do this in R.
Sarah
On Wed, Oct 31,
R-help -
I have this set of aggregated tables (sample data below via dput()). And I
would like to have delayValue as the column variables with the "temp"
(temp1, temp2, temp3) values as the row variables. However I would like to
have the temp variables *aggregated into single rows* so that I hav
Sorry, I know I should read a little 1st about this, but I am actually just
helping somebody really quick and need help too.
I want to grep all of the names of the .txt files mentioned on this html web
page:
http://www.epa.gov/emap/remap/html/three/data/index.html
Thanks ahead of time.
--
Vi
On Wed, Oct 31, 2012 at 4:31 PM, Hard Core wrote:
> Thank you but with that formula i don't understand how R can select only the
> one with blue eyes
>
Assuming you're referring to my proposal of
aggregate(height ~ eyes, DATA, mean)
It's not blue eyes only: that will summarize height for each e
A possiblie solution might be to use the survey package. You could specify
that the data is clustered using the svydesign function, and then speciy the
replicate weights using the as.svrepdesign function. And then, it would be
possible to use the withReplicates function to bootstrap the clusters
Thank you but with that formula i don't understand how R can select only the
one with blue eyes
--
View this message in context:
http://r.789695.n4.nabble.com/Mean-Help-tp4648000p4648030.html
Sent from the R help mailing list archive at Nabble.com.
_
Thank you ... can you explain what you've done?
--
View this message in context:
http://r.789695.n4.nabble.com/Mean-Help-tp4648000p4648013.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat
Hi Joshua,
The code you put together is very helpful. I have run into a small issue,
however, and I am surprised you weren't getting similar error message when
you tried running the nested for loops? As an example, I've pared down the
for loop a bit to highlight the error message I'm getting (no
Sorry about that, this is my first time posting. Below is the complete post.
I am trying to run the R Script below, I have actually simplified it to just
this part that is causing issues. When I run this script I continue to get an
error that says "cannot rescale a con
> Adam Clark wrote:
>>> I'll go ahead and dig through my C code.
>> My problem was that I had a pointer (*Aest) that had less memory allocated
>> to it than I ended up storing in it ...
On Oct 31, 2012, at 9:04 AM, William Dunlap wrote:
> As long as you use C (or C++ or Fortran ...), using me
R-site search with 'forward backward' returns references to packages for
fitting hidden Markov models.
hth, Ingmar
https://www.google.com/search?q=forward+backward&domains=r-project.org&sitesearch=r-project.org&btnG=Google+Search
On Wed, Oct 31, 2012 at 3:50 PM, quantum wrote:
> How can I wrti
?predict.xxx
where xxx is either lm or glm, as I don't know what you mean by
"general linear model." I also don't know what you mean by "summary
function" but have assumed it's "fitted model."
You need to read ?lm (or ?glm) carefully, and also follow some of the
links therein.
-- Bert
On Wed, O
No difference.
R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-w64-mingw32/x64 (64-bit)
> chooseCRANmirror()
> utils:::menuInstallPkgs()
trying URL 'http://cran.stat.ucla.edu/bin/windows/contrib/2
As long as you use C (or C++ or Fortran ...), using memory that
you don't own is possible. This is one reason people use languages
like R.
(If you program microprocessors or write operating system code
then C's willingness to let you read or write any at address is essential.)
Bill Dunlap
Spotfi
On Oct 31, 2012, at 5:47 AM, fillay89 wrote:
> I am trying to run the R Script below,
There is no script (or data description) appearing in this posting to the
Rhelp Mailing List. Nabble is not Rhelp, despite Nabble's effort to get you to
think it is an archive. Please review the Posting Guid
On Oct 31, 2012, at 5:42 AM, Eugenie wrote:
> tDate tTimeO3 No2 Temp Sun Wspeed Wdirect Hum Indicator
> 119980101 2400 0.065 0.036 31.4 7659.9 351 NA 1
> 219980102 2400 0.053 0.025 31.8 6247.7 351 NA 1
> 319980103 2400 0.027 0.033 31.5 852
Thierry, To answer your first question
summary(casedata)
station year monthday obs
mpe
Length:236 Min. :2011 Min. :8 Min. :28 Min.
:0. Min. :0.0
Class :character 1st Qu.:2011 1st Qu.:8 1st Qu.:28 1st
Qu
Aha - got it.
My problem was that I had a pointer (*Aest) that had less memory allocated
to it than I ended up storing in it (e.g. I had *Aest = 1:100, but stored
into it values at positions 5:105). With that fixed, all works flawlessly.
Thanks a lot for the help. What I hadn't realized was that
please provide a bit of the dataframe in the future using the dput() function
> x<-data.frame(height=1:10,
> color=sample(c("blue","green","brown"),10,replace=T))
> x
height color
1 1 blue
2 2 green
3 3 blue
4 4 brown
5 5 blue
6 6 brown
7 7 green
8
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of Eugenie
> Sent: Wednesday, October 31, 2012 5:42 AM
> To: r-help@r-project.org
> Subject: Re: [R] HELP!! how to remove 10% of data randomly in R
>
> tDate tTimeO3 No2 Tem
You forgot to install a JRE (Java runtime).
And please use R 2.15.2, which is current, if you are installing from
scratch. Some packages built with 2.15.2 (e.g. Matrix) do not work with
2.15.[01].
rJava is on CRAN/rforge.net, not 'R-force' (sic).
> PLEASE do read the posting guide
http://w
Many forms, but I'll recommend this one this time around:
aggregate(height ~ eyes, DATA, mean)
See ?aggregate for an explanation of what aggregate() does and
?formula for an explanation of the tilde syntax.
Note that this assumes your column names are "height" and "eyes."
Adjust as needed.
Mich
On Wed, Oct 31, 2012 at 2:50 PM, quantum wrote:
> How can I wrtie and calculate alpha and beta in the forward backward
> algorithm in R ?
You might need to be more specific Here's a smattering of links
you should probably look into as well:
http://r-manuals.flakery.org/R-intro.html
http://s
Hello,
It is said
> In addition: Warning messages:
> 1: package ‘lsa’ was built under R version 2.15.2
> 2: package ‘Snowball’ was built under R version 2.15.2
Install R 2.15.2 first.
Regards,
Pascal
Le 31/10/12 23:29, Triss.Ashton a écrit :
I just bought a new machine and installed the late
Dear All,
I'd like to set up a loop whereby successive parameter values are used as
start values in gamlss (yes I know this isn't usually necessary ! -
unfortunately for my truncated data it is), to return the estimated
parameters etc. giving the lowest AIC value.
I notice that mu.start can take
How can I wrtie and calculate alpha and beta in the forward backward
algorithm in R ?
--
View this message in context:
http://r.789695.n4.nabble.com/Forward-and-backward-algorithm-in-R-tp4648012.html
Sent from the R help mailing list archive at Nabble.com.
Hi there!
I need to generate a final general linear model from a set of coefficients
of the variables within the summary function.
When I have the summary-function, HOW do I create a final model (response
variable = bla bla bla)?
LordSword
--
View this message in context:
http://r.789695.n4.n
I have a dataframe.
Let's suppose that i have two columns. The first one contains height, the
second one contains eye color that can be Green, Blue or Brown.
I want to calculate the aritmetic mean of the height only for those people
who have Blue eyes. How can I do it?
Thank you for your availabili
Dear R-users,
Does somebody know what does the "Estimate" reported by the Lavaan model
tell us?
I assume this tells the relative strength of the dyadic relations.
Thank you for your help!
Regards,
Sylvain
--
View this message in context:
http://r.789695.n4.nabble.com/Lavaan-model-tp4648004.htm
I am trying to run the R Script below, I have actually simplified it to just
this part that is causing issues. When I run this script I continue to get
an error that says "cannot rescale a constant/zero column to a unit
variance". I cannot figure out what is going on here. I have stripped down
m
I just bought a new machine and installed the latest release of R 2.15.1 two
days ago. Loaded libraries yesterday (all reported successful loads).
While running scripts, rJava and Snowball would not load. Here is Snowball
successful install:
> utils:::menuInstallPkgs()
trying URL
'http://cran.
Thanks for the advice.
I'll go ahead and dig through my C code. It's helpful to know that my C
code can cause R to crash AFTER successfully implementing the code.
I have made sure to account for C's vector indexing, and I think I'm
allocating my C memory, and passing information back and forth be
Dear Adrienne,
What is the output of summary(casestudy) and summary(gridmeta)?
What happens if you set nmax to 10?
krige(formula=bias~1,locations=~lon+lat,data=casedata,newdata=gridmeta
,model=v.fit, nmax = 10)
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for
We would like to announce the following statistics course:
Introduction to linear mixed effects modelling course with R.
11 - 15 February 2013. Lisbon, Portugal.
For details, see: http://www.highstat.com/statscourse.htm
Course flyer: http://www.highstat.com/Courses/Flyer2013FebSIM_Lisbon.pdf
A piece of this is solved by the rms package's Rq and bootcov functions.
-Frank
Roger Koenker-3 wrote
> There is no automatic "clustering" option for QR bootstrapping.
> You will have to roll your own.
>
>
> url:www.econ.uiuc.edu/~rogerRoger Koenker
> email
> rkoenker@
>
My apologies to all, the package is gstat that I am using not geoR.
thanks!
A
On Tue, Oct 30, 2012 at 1:35 PM, Adrienne Wootten wrote:
> Greetings all,
>
> Ran into a strange problem with the krige function from geoR. The
> problem that I am having is that while the krige function seems to
> wor
The rank test inversion option that you are trying to use won't
work with only one coefficient, and therefore with univariate
quantiles, if you use summary(rq(rnorm(50) ~ 1, tau = .9), cov = TRUE)
you will have better luck.
url:www.econ.uiuc.edu/~rogerRoger Koenker
emailrkoen.
There is no automatic "clustering" option for QR bootstrapping.
You will have to roll your own.
url:www.econ.uiuc.edu/~rogerRoger Koenker
emailrkoen...@uiuc.eduDepartment of Economics
vox: 217-333-4558University of Illinois
fax: 217-244-66
Hi Thorn,
May be it is a bug in aggregate.formula(). Not sure about it.
The solutions that I gave and yours were returning the colnames as NA.
I guess this should be much better in terms of getting the result in one step.
do.call(aggregate,list(mF,d,mean))
# a b y z
#1
Dear Arun,
Thanks for your suggestion, that does the trick. Just because I'm curious,
where does the problem come from? I figured that wrapping the formula object in
brackets would work as well:
aggregate((mF), d, mean)
So I guess it has something to do with the scope of mF, or what is the roo
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
John Kane
Kingston ON Canada
> -Original Message-
> From: cat...@bas.ac.uk
> Sent: Tue, 30 Oct 2012 10:27:29 -0700 (PDT)
> To: r-help@r-project.org
> Subject: [R] box() doesnt work
>
> Hi,
>
> when
Hi,
Try this:
res<- aggregate(eval(mF),d,mean)
res
# a b NA NA
#1 1 A -1.48354978 -0.37141485
#2 2 A -0.08862713 0.35359250
#3 3 A 1.17519518 -0.47595290
#4 1 B 0.10214686 -0.70005131
#5 2 B 0.41185154 0.03707291
#6 3 B 0.20507062 -0.67946389
res1<-aggregate(cbind(y, z)
I am not sure I understand exactly what you want but does this do anything like
what you want?
library(ggplot2)
mydata <- data.frame(result = rnorm(100), group = rep(c("1", "2"), each =
50),
side = sample(c("L", "R"), 100, replace = TRUE) , dtime
= rep(1:10, each=10))
On Wed, 31 Oct 2012, Duncan Murdoch wrote:
On 12-10-30 11:13 PM, Adam Clark wrote:
I'm running R 2.15.1x64, though the same problem persists for 2.13.0x32 and
2.13.0x64.
I am trying to run compiled C code using the .C convention. The
code compiles without problems, dynamically loads within the
On 12-10-30 11:13 PM, Adam Clark wrote:
I'm running R 2.15.1x64, though the same problem persists for 2.13.0x32 and
2.13.0x64.
I am trying to run compiled C code using the .C convention. The
code compiles without problems, dynamically loads within the R
workspace with no problems, and even runs
On Wed, Oct 31, 2012 at 10:54 AM, vincent guyader
wrote:
>
> yes, it could be a solution, but i prefer a "simple" operator.. so i will
> build one.
See this thread for both advice and pitfalls:
https://stat.ethz.ch/pipermail/r-help/2012-June/316240.html
Michael Weylandt
_
Hi,
I am working with the RCurl package and I am using the curlPerform
function for an soap-query.
The problem is that the code is usually working well, but sometimes the
connection gets lost.
So I wrote a while-loop to repeat the query if anything might happened
so that the same query runs a
I had the same problem, but found this note in the MICE package documentation
(at http://cran.r-project.org/web/packages/mice/mice.pdf):
Added June 25, 2012: The currently implemented algorithm does not handle
predictors that are
specified as fixed effects (type=1). When using mice.impute.2l.norm(),
1 - 100 of 114 matches
Mail list logo
