. Any suggestions?
Thanks ahead of time for any help you can offer,
Austin H. Jones
Department of Mathematics
Wake Forest University
f<-function(v)
{
prodsub<-function(v,q){prod(v[q])}
t<-length(v)
C<-vector("list",t)
for (i in 1:t)
{C[[i]]<-combn(1:t,i)}
e<-rep(0,t)
Dear useRs,
I am trying to build my package (nonpareff) which deals with some models of
data envelopment analysis.
The building worked well, but checking complains when it tests examples.
Zipped nonparaeff.Rcheck is attached.
Following is the log.
-
a
Dear useRs,
I am trying to build my package (nonpareff) which deals with some models of
data envelopment analysis.
The building worked well, but checking complains when it tests examples.
Zipped nonparaeff.Rcheck is attached.
Following is the log.
-
a
ease contact us for more information. Comments
and suggestions will be very welcome.
Cheers,
Paul
---
Paul H. Geissler, Ph.D.
USGS Status & Trends of Biological Resources Program
Coordinator, National Park Monitoring Projec
Hi Scott,
Have you tried rgdal? http://cran.r-project.org/web/packages/rgdal/index.html
Rob
--
Rob Braswell
University of New Hampshire
On Jun 19, 2009, at 10:48 AM, Scott Saleska wrote:
I am looking for an easy way to import data in a GEOTIFF file into
Splus or R. Any suggestions?
Tha
$ sudo R
> install.packages(rgdal)
worked for me on my Ubuntu machine
but, I already had GDAL library installed with
$ sudo apt-get install libgdal1-1.5.0
--
B. H. Braswell (Rob)
University of New Hampshire
On Jun 24, 2009, at 9:24 AM, Luis Ridao Cruz wrote:
R-help,
I'm trying to
Hi expeRts,
I would like to calculate weighted mean by two factors.
My code is as follows:
R> tmp <- by(re$meta.sales.lkm[, c("pc", "sales")],
re$meta.sales.lkm[, c("size", "yr")], function(x)
weighted.mean(x[,1], x[,2]))
The result is as follows:
R>
Hello,
I've run 7 candidate models using mixed-effects logistic regression with
the lmer function from the lme4 package, and I'm getting the following
error for 5 of those models: Warning message: In mer_finalize(ans : false
convergence (8). The candidate models all run with the same data, just
I'm trying to compile R 2.10.1 on AIX, and getting the following errors when it
is trying to build the tools package:
Error in read.dcf(file = descfile) :
Line starting 'Package: tools ...' is malformed!
Calls: makeLazyLoading ... code2LazyLoadDB -> loadNamespace ->
parseNamespaceFile -> read
I need to run some conditional logistic regression models on a multilevel
matched dataset (propensity score matched data from multiple research sites).
I can pretty easily use the clogistic() function in the Epi package to run
conditional logistic regression models on each separate/nested subse
).
Is this understanding correct?
Thanks,
Heather
Heather Hensman Kettrey, PhD
Assistant Professor of Sociology
Clemson University
From: Wolfgang Viechtbauer
Sent: Wednesday, July 3, 2019 9:15:50 AM
To: Marc Schwartz
Cc: Heather H Kettrey; R-help
Subject: Re
3, 2019 1:05:29 PM
To: Marc Schwartz; Heather H Kettrey
Cc: R-help
Subject: Re: [R] Conditional Logistic Regression with Multilevel Data [using
clogistic() in Epi and lmer() in lme4???]
No need to install metafor for lme4.
Best,
Wolfgang
On July 3, 2019 5:19:13 PM GMT+02:00, Marc Schwartz
Hello everyone, I'm trying to upgrade R from 3.6.3 to 4.0.0 using the linux
terminal commands (sudo apt upgrade r-base r-base-dev) but I get a message
that says 3.6.3 is still the latest version. Please see the output below.
How could I fix this?
Thanks much for your help!
Carlos
===
At Thu, 27 Dec 2012 15:38:08 -0500,
Sam Steingold wrote:
> so,
> 1. is there a way for a function to modify a global variable?
Use <<- instead of <-.
> 2. how would you vectorize this loop?
This is hard. Your function has a feedback loop: an iteration depends
on the previous iteration's result.
Hi, Ivo,
At Fri, 28 Dec 2012 16:34:03 -0800,
ivo welch wrote:
> so, it isn't obvious to me how I get to the try-error object
> list, or how to find out what triggered the abort.
A try object is an object that captures the execution context at the
time of an error. Consider the following code:
At Sun, 30 Dec 2012 18:26:45 +0800 (CST),
meng wrote:
>
> hi all:
> Here's a dataframe(dat) and a vector(z):
>
> dat:
> x1 x2x3
> 0.2 1.2 2.5
> 0.5 2 5
> 0.8 3 6.2
>
> > z
> [1] 10 100 100
>
> I wanna do the following:
> 10*x1,100*x2,1000*x3
>
> My solution is us
ou mean multiplied by element.
Here's a better solution using t to transpose the matrix:
dat=data.frame(x1=1:3, x2=11:13)
as.matrix(dat)
x1 x2
[1,] 1 11
[2,] 2 12
[3,] 3 13
t(as.matrix(dat)) * c(2, 3)
[,1] [,2] [,3]
x1246
x2 33 36 39
>
> Andrius
>
At Mon, 31 Dec 2012 12:13:43 +0200,
catalin roibu wrote:
>
> Dear R users,
> I want to group numerical values in classes with different size and count
> the values for each classes.
>
> My data is in this forma:
> d 15 12,5 30,4 20,5 80,4 100,5 8,2 40,5 33 21 11
> And I want the group
At Mon, 31 Dec 2012 15:38:10 -0500,
Stephen P. Molnar wrote:
> The current version of R that is available in Wheezy is 2.15.1. However,
> version 2.15.2 is available at CRAN sites.
> ...
> My question is what should be the format of the line in the sources.lists?
Here's what I have:
deb http:/
At Mon, 31 Dec 2012 22:25:25 +,
Muhuri, Pradip (SAMHSA/CBHSQ) wrote:
> The issue is that, for Utah, I am getting an instead of (42,48.7] in the
> ob_mrj_cat column.
The problem is likely due to comparisons of floating point numbers.
Try moving your lower and upper bounds out a tiny bit. Whe
At Tue, 1 Jan 2013 02:00:14 +,
Muhuri, Pradip (SAMHSA/CBHSQ) wrote:
> Although David's solution (putting the right parenthesis, which I had missed)
> has resolved the issue, I would like to try yours as well.
>
> Could you please clarify the six elements: c(-1e-8, 0, 0, 0, 0, 1e8)?
It's a v
At Fri, 4 Jan 2013 17:17:40 +0530,
Christofer Bogaso wrote:
>
> Hello again,
>
> Let say I have 2 matrices which equal number of columns but different
> number of rows like:
>
> Mat1 <- matrix(1:20, 4, 5)
> Mat2 <- matrix(1:25, 5, 5)
>
> Now for each column 1-to-5 I need to fetch the correspond
At Fri, 15 Mar 2013 09:22:15 -0400,
Amir wrote:
> I have two sets T1={c1,c2,..,cn} and T2={k1,k2,...,kn}.
> How can I find the sets as follow:
>
> (c1,k1), (c1,k2) ...(c1,kn) (c2,k1) (c2,k2) (c2,kn) ... (cn,kn)
I think you are looking for expand.grid:
expand.grid(1:3, 10:13)
Var1 Var2
1
vif.lm(reg) : model contains fewer than 2 terms
Is there a solution or a way to work around?
Thank you very much in advanced.
--
Kind Regards,
Martin H. Schmidt
Humboldt University Berlin
__
R-he
Thanks for your suggestions. I did some further digging and found the problem
was more subtle. I use Dashlane, a password management program, which requires
IE to run in protected mode under Win 8. Disabling protected mode solved the
problem with R help.
Jack
> On May 26, 2014, at 12:02 AM, J
Thanks Gabor. I'll try that.
Jack
Sent from my iPad
> On May 26, 2014, at 2:40 PM, Gabor Grothendieck
> wrote:
>
> On Mon, May 26, 2014 at 8:53 AM, Jack H. Pincus
> wrote:
>> Thanks for your suggestions. I did some further digging and found the
>> problem
Hello,
I was wondering how to evaluate a string variable in R.
eg.
> avar <- "getwd()"
> eval(avar)
gives:
[1] "getwd()"
but I would like to see:
[1] "/home/myhomedir"
Thanks,
Jeremiah
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
Hello,
I want to perform a nonlinear curve fit in order to obtain parameter
estimates from experimentally determined data (y in dependence of x),
but with an implicit function, thus, a function of which I cannot
isolate y on the left-hand side of the equation. As far as I understand,
the function
Hello (and Happy New Year),
When I create a factor with labels in the order I want, write the data as a
text file, and then retrieve them, the factor levels are no longer in the
proper order.
Here is what I do (I tried many variations):
# educ is a numeric vector with 1,001 observations.
# The
Thanks to one and all. I now have a better understanding of the situation.
Original message
>Date: Sat, 1 Jan 2011 11:03:59 -0500
>From: David Winsemius
>Subject: Re: [R] Retrieving Factors with Levels Ordered
>To: h...@udel.edu
>Cc: r-help@r-project.org
>
>
&
How can I get a substring based on the index into the string?
strM <- c("abcde", "cdefg")
ind <- c(1,3,5)
I want to use ind to index into the strings so the result is:
strMind <- c("ace", "ceg")
__
R-help@r-project.org mailing list
https://stat.ethz.c
How can I get a substring based on the index into the string?
strM <- c("abcde", "cdefg")
ind <- c(1,3,5)
I want to use ind to index into the strings so the result is:
strMind <- c("ace", "ceg")
__
R-help@r-project.org mailing list
https://stat.ethz.c
It is a simple problem in that I simply want to convert the For loop to a more
efficient method. It simply loops through a large vector checking if the
numeric element is not equal to the index of that element. The following
example demonstrates a simplified example:
> rows <- 10
> collusionM
I need to increment cells of a matrix (collusionM). The indexes to increment
are in an index (matchIndex). This is some of the code
for (j in 1:(rows-1)) matchIndex[[j]] <- which(mx[j]==mx)
for (j in 1:(rows-1)) collisionM[j,matchIndex[[j]]] <-
collisionM[j,matchIndex[[j]]] + 1
I need to increment cells of a matrix (collusionM). The indexes to increment
are in an index (matchIndex). This is sample code
library(seqinr)
library(Matrix)
x <- "abcabcabc"
mx <- s2c(x)
collisionM <- Matrix(0,nrow=10, ncol=10, sparse=TRUE)
matchIndex <- list()
rows <-
I am working with a large sparse matrix trying replace all 1 values with 0.
The normal method doesn't work. Here is a small example:
> x <- Matrix(0,nrow=10,ncol=10,sparse=TRUE)
> x[,1] <- 1:2
> x
10 x 10 sparse Matrix of class "dgCMatrix"
[1,] 1 . . . . . . . . .
[2,] 2 . . . . . . . . .
[3,]
Anyone have experience specifying the "penalty.factor" option in the
"glmnet" command?
I have 3 variables (out of a million genotype variables) that I want to
force into the model (i.e., set penalty factor to 0), but I can't figure out
how to do that.
[[alternative HTML version deleted]]
rted
to "AB". I tried with car but it very slow as I need to very large
dataframe.
Thanks;
--
Ram H
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
t;non-obvious" use of '[' as
R> ct1b <- sapply(ct1.char, '[', 2) ## a function is intentional :-)
Can anybody help me to get solution out of it?
Ram H
On Sun, May 1, 2011 at 5:03 PM, Steve Lianoglou <
mailinglist.honey...@gmail.com> wrote:
> Hi,
>
&g
, y = list(rot= 90, font = 1, cex =
1) ))
The output is not what I want. I want the orientation of graph like the
following in base R but axis label are in Y axis line and other parameters
as in lattice:
barplot(Y, names.arg = X)
I know this is simple question, but I could not find a true solution.
Dear list,
I'm working with a MySQL-Table and would like to extract selected ids that are
stored in a vector:
Ids <- c(1,2,3)
How do I insert this variable in a for loop like this:
For (a in values) {
Temp <- dbGetQuery(con,
"select * from db
Where
Dear R users
I am new R user, execuse me I bother you, but I worked hard to find a
solution:
# data
ID <- c(1:100)
set.seed(21)
y <- rnorm(100, 10,2)
x1 <- rnorm(100, 10,2)
x2 <- rnorm(100, 10,2)
x3 <- rnorm(100, 10,2)
x4 <- rnorm(100, 10,2)
x5 <- rnorm(100, 10,2)
x6 <- rnorm(100, 10,2)
mydf <- d
Dear R-HELP,
We are compiling R-2.12.1 for 64 bits onto a Sun Sparc machine below:
SunOS 5.10 Generic_142900-13 sun4u sparc SUNW,Netra-T12
The source was compiled successfully with "gcc version 3.4.6" by using the
default configuration and has produced R-2.12.1 in 32 bits.
Compiling R-2.12.1 fo
Dear R-HELP,
We are compiling R-2.12.1 for 64 bits onto a Sun Sparc machine below:
SunOS 5.10 Generic_142900-13 sun4u sparc SUNW,Netra-T12
The source was compiled successfully with "gcc version 3.4.6" by using the
default configuration and has produced R-2.12.1 in 32 bits.
Compiling R-2.12.1 fo
Dear R-Experts
I am sure this might look simple question for experts, at least is problem
for me. I have a large data frame with over 1000 variables and each have
different distribution( i.e. have different quantile). I want to create a
new grouped data frame, where the new variables where the val
Dear R community members
I have been struggling on this simple question, but never get appropriate
solution. So please help.
# my data, though I have a large number of variables
var1 <- rnorm(500, 10,4)
var2 <- rnorm(500, 20, 8)
var3 <- rnorm(500, 30, 18)
var4 <- rnorm(500, 40, 20)
datafr1 <- da
-3.736066 -3.471752
> [2,] 17.37176 34.929800 62.969733 80.224799
> > apply(datafr1, 2, range)
> var1 var2 var3 var4
> [1,] -2.668841 -10.23848 -22.13500 -22.93678
> [2,] 21.803714 40.48186 87.02953 91.92733
>
> Assuming you wanted to do this columnwise
n each element and it is also in a form that you can easily
> process it to create whatever other output you might need.
>
> On Fri, Mar 18, 2011 at 7:24 AM, Ram H. Sharma
> wrote:
> > Hi Dennis and R-users
> >
> > Thank you for more help. I am pretty close, but challeng
in each element and it is also in a form that you can easily
> process it to create whatever other output you might need.
>
> On Fri, Mar 18, 2011 at 7:24 AM, Ram H. Sharma
> wrote:
> > Hi Dennis and R-users
> >
> > Thank you for more help. I am pretty close, but chal
Dear all
My dataframe has > 80,000 variables which I can not everytime load into R
using *.txt files (read.table option), cost me time and sometime computer
decomes not responsive. So I need a way to save my dataframe in my
workdirectory as such.
Execuse me if the problem is too simple. I tried t
Thank you David and Dennis
On Mon, Mar 28, 2011 at 3:43 PM, David Winsemius wrote:
>
> On Mar 28, 2011, at 11:20 AM, Ram H. Sharma wrote:
>
> Dear all
>>
>> My dataframe has > 80,000 variables which I can not everytime load into R
>> using *.txt files (rea
at my data.frame would
look like the follows:
Var1 Var2 Var3 Var4
0 0.05 0.0112
1 0.04 0.0618
2 0.05 0.0814
3 0.01 0.0615
4 0.05 0.07 14
--
Thank you in advance for t
etermined by PR1[1] indicator
determines whether to pick lcd1[1] or lcd2[1] or lcd3[1] element
similarly for PR1[2] indicator whether to pick lcd1[2] or lcd2[2] or lcd3[2]
element
The same process need to be continued other PR2 to PR10 variables. That's
why I attempted to use sapply.
Thank you f
13
5 22 11 12 12 11 11 22
Thank you;
Ram H
On Sat, Apr 9, 2011 at 4:13 PM, Phil Spector wrote:
> Ram -
> I think you'll have to explain what you're trying to do.
> First, you're creating a 10x10 matrix newd, which could
> be done a bit more efficiently by using
to mixed model forum but
I did not get any response. Sorry to post all of you, but my hope is my
question is simple enough and bigger R community can help me !
Ram H
--
Ram H
[[alternative HTML version deleted]]
__
R-help@r-project.org
;")
How can I read this into table where columns first five column are
maintained as such (e.g. 1_1 ch1 0.0 H ) where as the long single column
(filled with A, H, B) split into individual columns. Or alternatively first
read using the following command and then able to split the variable
Sorry NULL in V8 should be corrected as NA
On Wed, Apr 13, 2011 at 7:33 AM, Ram H. Sharma wrote:
> It looks like my question is not clear, I have not get any suggestion, yet
> let me reiterate my problem:
>
> My data looks like this to be read from a text file. As I provided earlie
It looks like my question is not clear, I have not get any suggestion, yet
let me reiterate my problem:
My data looks like this to be read from a text file. As I provided earlier
the A, B, and H in V4 column has much longer chain.
V1 V2 V3V4
1_1ch1 0.0 AHAH
1_2
ct number of dimensions
I can extract fixed effects as well as I do in lme4
b <- fixef(fit1)
Error in UseMethod("fixef") :
no applicable method for 'fixef' applied to an object of class "lmekin"
--
Ram H
__
R-help@r-proj
gt;
> P <- vector (mode="numeric", length = 1000)
>
> P[i] <- lmekin(yvar~ mydata[,i] , data= mydata, random = ~1|id,
> varlist=list(kmat)) $coefficients[2,4]
>
> }
>
> Same errors: I tried lme4 conventioned but did not work !
> I can extract fixed
me to go further in the function and out me the error.
You help is highly appreciated.
Thanks;
--
Ram H
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Thank you Ted and David for prompt reply.
I can accept Inf but can not use for plotting, which I intend to do. May be
I can add 1/(10^308), so that if something comes to 0 will be in -log10
scale be 308.
Ram H
On Sun, Aug 28, 2011 at 12:12 PM, Ted Harding wrote:
> On 28-Aug-11 15:37:06, Ra
i.e. X1 and X2 has correlation of
0.5, X2 and X3 has also correlation of 0.5 and X3 and X4 has correlation of
0.5 and so on.
I have already spend considerable time on it, that is how I got to the
group, sorry for that.
--
Ram H
[[alternative HT
The problem remain unsolved. If you have any idea please do suggest
thank you;
Ram H
On Wed, Aug 31, 2011 at 12:06 PM, Ram H. Sharma wrote:
> Dear R experts:
>
> I have following problem:
>
> # myfunction
> mfun1 <- function(x) {
>if ( x == 2)
save individual Rdata sets and load it using
load ("my.Rdata") function. It is more fast this way. But could not figure
out how to do it !!!
In above function I have put # some of my attempt to solve this problem
Any suggestions
further questions on the Bioconductor mailing list
http://bioconductor.org/docs/mailList.html
Error: package/namespace load failed for 'Rgraphviz'
Help me please.
--
Ram H
[[alternative HTML version deleted]]
__
R-help@r-pr
ted once you've
installed graphviz and before installing Rgraphviz"
-- Forwarded message --
From: Martin Morgan
Date: Mon, Jul 25, 2011 at 9:25 AM
Subject: Re: [R] Rgraphviz installation problem
To: "Ram H. Sharma"
Cc: r-help@r-project.org
On 07/25/2011 06:05 AM, R
Hi Martin
Thank you for help, that helped me to come out the problem !
On Mon, Jul 25, 2011 at 11:46 AM, Martin Morgan wrote:
> On 07/25/2011 07:58 AM, Ram H. Sharma wrote:
>
>> Hi All
>>
>> How can I adjust this path. My current path of Graphviz installation is:
>
;red")
}, as.table=T, subscripts=T)
Thank you;
--
Ram H
[[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/
;red", ylim= c(0,10)
)
panel.lmline(iris$Petal.Length [subscripts], iris$Petal.Width
[subscripts], col = "red")
}, as.table=T, subscripts=T)
Help please:
From: Dieter Menne
> Date: Wed, Jul 27, 2011 at 8:44 AM
> Subject: Re: [R] lattice overlay
> To: r-help@r-project.o
To whomever it may concern,
I'm a young Industrial Engineer working on Senior Design at Georgia Tech
and have found the StructTS method to be excellent for the training set for my
forecasting project. There's only one problem: I don't actually understand what
a Structural Time Series IS
Dear R,
Thanks for helping me locate the source for the StructTS method from
stats, but I've run in to a roadblock in reverse engineering it to locate a
formula for its forecasting because it calls some compiled C code, a function
called KalmanLike. I've looked through that R library that
Hi, All:
I am new to R and tm package. I'm trying to do the stemming using tm_map()
and it doesn't seem to work:
*I used:*
> stemDocument(t_cmts[[100]])
*Where t_cmts is the corpus object, the results is:*
bottle loose box abt airpak sections top plastic bottle squashed nearly
flush neck previo
You might look at John Kruschke's book, Doing Bayesian Data Analysis (AP),
which starts with basics and goes from there. It also relies on R and Bugs.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
Matthew,
My Fault I should have sent you to the current release cycle page. The link was
for the old 2.6 Bioconductor release, sorry about that. :(
Personally the easies way to install any bioconductor package is to open R copy
and paste the following code:
source("http://bioconductor.org/bioc
Lars Fischer sec.informatik.tu-darmstadt.de> writes:
> Hello,
>
> well, I searched list-archive, cran and the references, but found
> nothing. Thus:
>
> Does anybody around here know anything about Dempster-Shafer Theory,
> Evidence Theory or Hints in R? Has anybody stumbled about a package tha
Hi all,
i'm trying to install JGR on linux (see bottom), but it runs into errors.
The bottom of the sequence below shows "JavaJD could not be found"
Any suggestion, please?
$ sudo apt-get install sun-java6-jdk
$ sudo update-alternatives --config java
There are 5 alternatives which provide `j
Thanks.
I basically followed the steps in http://lib.stat.cmu.edu/R/CRAN/, and those
below for JGR and it worked.
...but there are problems with help files e.g. ?matrix generates the
following message:
/tmp/RtmpvxsnQa/.R/doc/html/packages.html (No such file or directory)
g-h-d wrote
Dirk Eddelbuettel wrote:
>
>
> On 3 June 2008 at 08:42, g-h-d wrote:
> | I basically followed the steps in http://lib.stat.cmu.edu/R/CRAN/, and
> those
> | below for JGR and it worked.
> |
> | ...but there are problems with help files e.g. ?matrix generates t
Dirk Eddelbuettel wrote:
>
>
> On 3 June 2008 at 13:44, g-h-d wrote:
> | > | > $ sudo apt-get install sun-java6-jdk
> | > | >
> | > | > $ sudo update-alternatives --config java
> | > | >
> | > | > There are 5 alternatives which prov
Hello,
I recently upgraded from R 2.6.1 to 2.6.2. I uninstalled 2.6.1, installed
2.6.2, and installed the packages I regularly use. Everyting seems to be
running properly, including most of the packages, but when I try to call
the Matrix package, I get the following error:
> require(Matrix)
Load
An undergraduate here, so do not hesitate to let me know if you feel that
I'm heading in a wrong direction.
I have a data frame containing panel data across 10 years (hence 120
months). I want to be able to run regression separately for each month (or
year). The below shows how I ran the regr
(mydf, mydf$TIME), function(x){coefficients(summary(glm(X ~ A +
B, data = x, family=binomial)))})
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Mike H. Ryu
Sent: Tuesday, April 29, 2008 7:27 AM
To: r-help@r-project.org
Subject: [R] Running regression (lm
-0i
[3,] 100680+0i 32520-0i
[4,] 97560-0i 28800-0i
[5,] 88800-0i 33000+0i
[6,] 69840+0i 44400+0i
>?real
> real(pb)
Error in vector("double", length) :
unimplemented type 'complex' in 'asVecSize'
> real(pb[,1], length(dim(pb)[1]))
Error in real(pb[, 1],
I'm trying to drop all rows except for the ones with the most recent year.
So I split the data frame by NPERMNO and keep just the last record of all
groups.
datg=t(sapply(split(datgic, datgic$NPERMNO, drop=TRUE), function(x){return(
x[nrow(x),] )}))
I get something like this...
GVKEY NPE
What is the command for 2D fft? Or 2D convolution and or auto-correlation?
Thanks,
Paul
__
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
Yea the anova object seems to be odd. It's not S4 so that's why I tried
originally the attr() funtion but
summary(pb)$Pr(>F)
Error: unexpected '>' in "summary(pb)$Pr(>"
> summary(pb)$Pr
NULL
> summary(pb)@Pr(>F)
Error: unexpected '>' in "summary(pb)@Pr(>"
> summary(pb)@Pr(F)
Error: no slo
yep That worked Jorge. Thanks!
summary(pb)$"Pr(>F)1"
NULL
> summ<-summary(pb)
> class(unlist(summ))
[1] "numeric"
> unlist(summ)["Pr(>F)1"]
Pr(>F)1
0.02533637
> names(summ)
NULL
Jorge Ivan Velez wrote:
>
> Hi Paul,
>
> Perhaps
>
> # Data set
> set.seed(123)
> x=rnorm(100)
> y=x+2*rnorm(10
e library, but get bogged down in the equation. Also, it doesn't
seem like polymars will generate a tree even if I do get it working. Can rpart
be modified in some way to accomodate multiple response parameters? If
anyone's ever come across this situation before, pointers would be much
app
Thank you for the reply. With the improved formula, mvpart worked like a charm.
Sincerely,
Jeff
Jeff H Bardwell, M.S.
Biology Department
ENV 1101 Lab Coordinator
Goebel 115, OH: Thu 1pm-4pm
710-6596 (e-mail preferred)
From: Gavin Simpson [mailto:[EMAIL
Hi all,
I have a quick question, I want to send a select statement to the
mysql server. With the code below I connect and ask my select statement.
As you can see from the code I want to have a mode = "+". The '+' needs
to be free text and hence the "\"" around it. But it doesn't work. If i
Yea that worked !!
Did
paste(SELECT * FROM MET WHERE molid=", molid, "AND mode=", "\'" mode,
"\'" , sep="")
I guess all I needed was the sep="" which I totally forgot about.
Cheers,
Paul
Eric wrote:
> Does this work for you?
>
> "SELECT * FROM data WHERE molNUM = 165 AND mode = '+'"
>
> Possibl
I am running a series of candidate mixed models using lmer (package lme4)
and I'd like to be able to compile a list of the AIC scores for those
models so that I can quickly summarize and rank the models by AIC. When I
do logistic regression, I can easily generate this kind of list by creating
the
On my current home system, I am getting undesirable output from
graphical devices such as png() and pdf(). The graphical output is
blurry. I haven't experienced the problem on other systems. As you
will see from the attached text file (more information on this file
below), the problem does not occu
On Wed, Dec 17, 2008 at 12:54 AM, Prof Brian Ripley
wrote:
> Your PDF problems indicate a broken viewer. How were you viewing PDF?
You are absolutely correct about the PDF files; I've since checked the
PDF files in other viewers and have not been able to reproduce the
problem. There was definite
On Wed, Dec 17, 2008 at 3:06 AM, Martyn Plummer wrote:
> The artefacts that you see are a normal result of using bitmap graphics
> devices. I have tried to explain these below:
Thanks very much for your explanations, MP; they were quite informative!!
I recognize that others may feel differently
Hello, I have two matrices as shown below:
Matrix 1
IDAB1BC1CD1...
name1 1,1 2,1 0,2 ...
name2 2,0 1,2 1,2 ...
name3 0,2 1,1 2,0 ...
name4 2,0 0,2 0,2
TY OF
ANY KIND. Nonetheless, the authors will endeavor to fix bugs promptly
and to add requested features. Send bug reports, suggestions, and
extensions to either author.
Michael H. Prager - [EMAIL PROTECTED]
Southeast Fisheries Science Center
National Marine Fisheries Service, NOAA
101 Pivers Is
I am trying to generate a graphic with a matrix of 9 line graphs (3 rows, 3
columns), all with the same y-axis range, and only showing the y-axis
labels and title at the left edge of each row of the matrix. I have been
trying to use the widths argument in layout to specify a larger column
width in
Hello all,
I apologize if this is simple or has already been answered somewhere, but
I'm not sure what to search for although I have tried and didn't come up
with anything so.. Here's my question.
How can I interpolate list names or do I have to do it post list creation.
Since that's not v
101 - 200 of 306 matches
Mail list logo
