Displaying 20 results from an estimated 4000 matches similar to: "generate combination set"
2007 Jan 09
5
a question of substitute
Hi all,
I want to write a wrapper for an analysis of variance and I face a curious
problem. Here are two different wrappers:
fun.1 <- function(formula) {
summary(aov(formula))
}
fun.2 <- function(formula) {
oneway.test(formula)
}
values <- c(15, 8, 17, 7, 26, 12, 8, 11, 16, 9, 16,
24, 20, 19, 9, 17, 11, 8, 15, 6, 14)
group <- rep(1:3, each=7)
# While the first
2006 Nov 20
4
for help about logistic regression model
I have a dataset like this:
p aa
index x y z sdx sdy sdz delta as
ms cur sc
1 821p MET 1 -5.09688 32.8830 -5.857620 1.478200 1.73998 0.825778
13.7883 126.91 92.37 -0.1320180 111.0990
2 821p THR 2 -4.07357 28.6881 -4.838430 0.597674 1.37860 1.165780
13.7207 64.09 50.72 -0.0977129 98.5319
3 821p GLU 3 -5.86733 30.4759
2007 Jan 30
6
jump in sequence
Dear list,
This should be a simple one, I just cannot see it.
I need to generate a sequence of the form:
4 5 6 13 14 15 22 23 24
That is: starting with 4, make a 3 numbers sequence, jump 6, then another 3
and so on.
I can create a whole vector with:
myvec <- rep(rep(c(F, T, F), rep(3, 3)), 3)
Then see which are TRUE:
which(myvec)
[1] 4 5 6 13 14 15 22 23 24
I'd like to avoid
2006 Nov 26
1
plot p(Y=1) vs as
I am trying to fit a logistic regression model for this data set.
Firstly, I want to plot P(Y=1) vs As and P(Y=1) vs Aa.
Does any body know how to do these in R.
Thanks,
Aimin
> p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv")
> str(p5)
'data.frame': 1030 obs. of 6 variables:
$ P : Factor w/ 5 levels "821p","8ABP",..: 1
2006 Dec 07
1
svm plot question
I run the following code, all other is ok,
but plot(m.svm,p5.new,As~Cur) is not ok
Anyone know why?
install.packages("e1071")
library(e1071)
library(MASS)
p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv")
p5.new<-subset(p5,select=-Ms)
p5.new$Y<-factor(p5.new$Y)
levels(p5.new$Y) <- list(Out=c(1), In=c(0))
attach(p5.new)
2006 Nov 30
4
R_WinEdt question
if I want to put fig1plot to the left, figYPplot to the right
figYAaplot on the bottom.
How to modify the following cod to do these?
\begin{figure}[H]
\centering
\begin{minipage}[t]{0.5\textwidth}
\centering
<<label=fig1plot,fig=TRUE,echo=FALSE>>=
<<fig1plot>>
@
%\caption{Caption 1}
\end{minipage}
\begin{minipage}[t]{0.5\textwidth}
\centering
2006 Aug 09
3
objects and environments
Dear list,
I have two functions created in the same environment, fun1 and fun2.
fun2 is called by fun1, but fun2 should use an object which is created in fun1
fun1 <- function(x) {
ifelse(somecondition, bb <- "o", bb <- "*")
## mymatrix is created, then
myresult <- apply(mymatrix, 1, fun2)
}
fun2 <- function(idx) {
if (bb == "o) {
#
2007 Jan 28
2
nnet question
Hello,
I use nnet to do prediction for a continuous variable.
after that, I calculate correlation coefficient between predicted value and
real observation.
I run my code(see following) several time, but I get different correlation
coefficient each time.
Anyone know why?
In addition, How to calculate prediction accuracy for prediction of
continuous variable?
Aimin
thanks,
> m.nn.omega
2007 Jan 21
2
multiple bases to decimal (was: comparing two matrices)
Hi again,
I was contemplating the solution using base 3:
set.seed(3)
mat2 <- matrix(sample(0:2, 15, replace=T), 5, 3)
Extracting the line numbers is simple:
bases <- c(3, 3, 3)^(2:0) # or just 3^(2:0)
colSums(apply(mat2, 1, function(x) x*bases)) + 1
[1] 7 23 25 8 1
The problem is sometimes the columns have different number of levels, as in:
mat1 <- expand.grid(0:2, 0:2,
2007 May 29
2
pie initial angle
Dear all,
I'd like to produce a simple pie chart for a customer (I know it's bad but
they insist), and I have some difficulties setting the initial angle.
For example:
pie(c(60, 40), init.angle=14)
and
pie(c(80, 20), init.angle=338)
both present the slices in the same direction, where:
pie(c(60, 40))
pie(c(80, 20))
present the slices in different directions.
I read everything I
2007 Jan 19
2
split data set
I have a data(ABC) like this:
x y
A 3 4
A 1 3
B 2 6
B 4 8
C 5 4
C 6 7
I want to split this data into
A:
x y
A 3 4
A 1 3
B
B 2 6
B 4 8
C
C 5 4
C 6 7
anyone knows how to do that?
thanks,
Aimin Yan
2007 Jul 27
4
Looping through all possible combinations of cases
Hello!
I have a regular data frame (DATA) with 10 people and 1 column
('variable'). Its cases are people with names ('a', 'b', 'c', 'd',
'e', 'f', etc.). I would like to write a function that would sum up
the values on 'variable' of all possible combinations of people, i.e.
1. I would like to write a loop - in such a way that it
2005 Nov 21
4
attributes of a data.frame
Dear all,
I noticed that a data.frame has four attributes:
- names
- row.names
- class
- variable.labels
While one can use the first three (i.e. names(foo) or class(foo)), the fourth
one can only be used via:
attributes(foo)$variable.labels
(which is kind of a tedious thing to type)
Is it or would be possible to simply use:
variable.labels(foo)
like the first three attributes?
I tried:
varlab
2006 Oct 24
2
for help
I have a question in R.
In directory H:/Delta_angle
I have 19 files like this:
ALA.delta
ASN.delta
ASP.delta
CYS.delta
GLN.delta
GLU.delta
HIS.delta
ILE.delta
LEU.delta
LYS.delta
MET.delta
PHE.delta
PRO.delta
SER.delta
THR.delta
TRP.delta
TYR.delta
VAL.delta
I want to read these files to 19 data sets in R.
All these data sets have "P","AA",index","delta"
2007 Nov 16
2
expand.grid overflows?
>cbn<-as.matrix(expand.grid( rep( list(0:1), 50)))
Error in rep.int(rep.int(seq_len(nx), rep.int(rep.fac, nx)), orep) :
invalid 'times' value
In addition: Warning message:
In rep.int(rep.int(seq_len(nx), rep.int(rep.fac, nx)), orep) :
NAs introduced by coercion
But I'm only interested in cbn matrix rows where:
cbn<- cbn[rowSums(cbn)==5,]
Is there a way to evaluate it
2004 Mar 02
4
how to delete a matrix column
Hello,
I am new to R, How to delete a matrix column.
Thanks,
2006 Nov 23
2
Sweave question
I try Sweave
and get Sweave-test-1.tex
but hot to run LaTeX on 'Sweave-test-1.tex'?
I am using WinEdt.
thanks,
Aimin
> Sweave(testfile)
Writing to file Sweave-test-1.tex
Processing code chunks ...
1 : print term verbatim
2 : term hide
3 : echo print term verbatim
4 : term verbatim
5 : echo term verbatim
6 : echo term verbatim eps pdf
7 : echo term verbatim eps pdf
2004 Oct 13
3
one more Rcmdr problem
Hello,
I'm using R 2.0.0 with the latest Rcmdr package installed from CRAN, on
Windows XP Professional.
When trying to copy some commands or results, either from the upper or lower
text window, this causes Rcmdr to crash:
"R for Windows GUI front-end has encountered a problem and needs to close"
Did anyone have the same problem? I don't think it's my system, as it
2004 Oct 13
3
one more Rcmdr problem
Hello,
I'm using R 2.0.0 with the latest Rcmdr package installed from CRAN, on
Windows XP Professional.
When trying to copy some commands or results, either from the upper or lower
text window, this causes Rcmdr to crash:
"R for Windows GUI front-end has encountered a problem and needs to close"
Did anyone have the same problem? I don't think it's my system, as it
2007 Feb 22
4
Sorting rows of a binary matrix
Hallo,
The command:
x <- 3
mat <- as.matrix(expand.grid(rep(list(0:1), x)))
generates a matrix with 2^x columns containing the binary representations
of the decimals from 0 to (2^x-1), here from 0 to 7. But the rows are not
sorted in this order.
How can sort the rows the ascending order of the decimals they represent,
preferably without a function which converts binaries to decimals