Displaying 20 results from an estimated 200 matches similar to: "Linear separation"
2006 Jun 14
1
Estimate region of highest probabilty density
Estimate region of highest probabilty density
Dear R-community
I have data consisting of x and y. To each pair (x,y) a z value (weight) is assigned. With kde2d I can estimate the densities on a regular grid and based on this make a contour plot (not considering the z-values). According to an earlier post in the list I adjusted the kde2d to kde2d.weighted (see code below) to estimate the
2009 Mar 08
1
Summary of data.frame according to colnames and grouping factor
A dataframe holds 3 vars, each checked true or false (1, 0). Another
var holds the grouping, r and s:
### start:example
set.seed(20)
d <- data.frame(sample(c(0, 1), 20, replace=T), sample(c(0, 1), 20,
replace=T), sample(c(0, 1), 20, replace=T))
names(d) <- c("A", "B", "C")
e <- rep(c("r", "s"), 10)
### end:example
How do I get the
2010 Nov 17
2
slicing list with matrices
A list contains several matrices. Over all matrices (list elements) I'd like to access one matrix cell:
m <- matrix(1:9, nrow=3, dimnames=list(LETTERS[1:3], letters[1:3]))
l <- list(m1=m, m2=m*2, m3=m*3)
l[[3]] # works
l[[3]][1:2, ] # works
l[[1:3]][1, 1] # does not work
How can I slice all C-c combinations in the list?
S?ren
--
S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag,
2009 Mar 07
2
Recode factor into binary factor-level vars
How to I "recode" a factor into a binary data frame according to the
factor levels:
### example:start
set.seed(20)
l <- sample(rep.int(c("locA", "locB", "locC", "locD"), 100), 10,
replace=T)
# [1] "locD" "locD" "locD" "locD" "locB" "locA" "locA" "locA"
2009 Nov 13
1
shrink list by mathed entries
Hello
a <- c("Mama", "Papa", "Papa; Mama", "", "Sammy; Mama; Papa")
a <- strsplit(a, "; ")
mama <- rep(F, length(a))
mama[sapply(a, function(x) { sum(x=="Mama") }, simplify=T) > 0] <- T
papa <- rep(F, length(a))
papa[sapply(a, function(x) { sum(x=="Papa") }, simplify=T) > 0] <- T
# ... more
2009 Mar 11
3
chisq.test: decreasing p-value
A Likert scale may have produced counts of answers per category.
According to theory I may expect equality over the categories. A
statistical test shall reveal the actual equality in my sample.
When applying a chi square test with increasing number of repetitions
(simulate.p.value) over a fixed sample, the p-value decreases
dramatically (looks as if converge to zero).
(1) Why?
(2) (If
2009 Nov 22
3
Define return values of a function
I have created a function to do something:
i <- factor(sample(c("A", "B", "C", NA), 793, rep=T, prob=c(8, 7, 5,
1)))
k <- factor(sample(c("X", "Y", "Z", NA), 793, rep=T, prob=c(12, 7, 9,
1)))
mytable <- function(x){
xtb <- x
btx <- x
# do more with x, not relevant here
cat("The table has been created,
2008 Jun 25
1
help_transformation
heya,
i am fitting linear mixed effect model to a response Y. Y shows an s-shaped distribution when using QQ-plots (some zero values and some very high values). hence, which transformation should i apply that Y follows a normal distribution? any r-function/package available to do this?
thanks for any hint,
regards,
lukas
???
Lukas Indermaur, PhD student
eawag / Swiss Federal Institute of
2009 Nov 13
0
Craddock-Flood Test in R?
Hello
The "Craddock-Flood Test" is recommended for large tables with small
degrees of freedom and low-frequency cells. Is there an R procedure
and/or package which does the test?
Thank you for your help!
S?ren Vogel
--
S?ren Vogel, Dipl.-Psych. (Univ.), PhD-Student, Eawag, Dept. SIAM
http://www.eawag.ch, http://sozmod.eawag.ch
2010 Dec 27
0
Nagelkerke R square for Prediction data
Hello
I found some small postings dated to 22 Oct 2008 on the message subject. Recently, I have been working with binary logistic regressions. I didn't use the design package. Yet, I needed the "fit" indices. Therefore, I wrote a small function to output the Nagelkerke's R, and the Cox-&-Snell R from a fitted model. I am no professional programmer by far, yet, I hope, that
2010 Jul 23
0
ROpenOffice (which requires Rcompression)
Hello, for my data preparation and administration (data, labels, etc.) I use OpenOffice.org ODS spreadsheet files with several sheets in one file. However, I find it inconvenient to export every single sheet to a csv file whenever I apply changes to the labelling or so, and I haven't found a plugin or an application which does this batch in OOorg. Is there new development on a direct import
2007 Feb 20
1
testing slopes
Hello
Instead of testing against 0 i would like to test regression slopes against -1. Any idea if there's an R script (package?) available.
Thanks for any hint.
Cheers
Lukas
???
Lukas Indermaur, PhD student
eawag / Swiss Federal Institute of Aquatic Science and Technology
ECO - Department of Aquatic Ecology
?berlandstrasse 133
CH-8600 D?bendorf
Switzerland
Phone: +41 (0) 71 220
2007 Feb 28
1
bootstrap
Hi,
I would like to evaluate the frequency of the variables within the best selected model by AIC among a set of 12 competing models (I fit them with GLM) with a bootstrap procedure to get unbiased results. So I would ike to do the ranking of the 12-model-set 10'000 times separately and calculate the frequency of variables of the 10'000 best ranked models. I wrote a script doing the model
2017 Nov 18
0
Using cforest on a hierarchically structured dataset
Hi,
I am facing a hierarchically structured dataset, and I am not sure of
the right way to analyses it with cforest, if their is one.
- - BACKGROUND & PROBLEM
We are analyzing the behavior of some social birds facing different
temperature conditions.
The behaviors of the birds were recorder during many sessions of 2 hours.
Conditional RF (cforest) are quite useful for this analysis
2006 Dec 01
1
mixed effects model and r-squared
Heya
I am fitting linear mixed effects model in R and want to assess the model fit (with Animal number as random factor; repeated measures for Animals):
ts.model <- lme(LOG_FOC_MW ~ R_DN_SUM + ANIMAL + SEX+ YY, data = t.data, random = ~ 1 | ANIMAL, correlation=corCAR1(0.2, form = ~1 | ANIMAL ), method='ML', na.action=na.omit)).
Is there a possability to easly compute an
2009 Jan 23
2
Write to multiple connections or multiple text files
Hi all,
I want to modify a large number of text files (ca 4000) by replacing a
value found on a particular line in them with a value from an R object.
For a single file I would normally use:
con<-file ("foo.txt", open="r+")
content<-readLines(con)
content[n]<-"test"
writeLines(content,con)
close(con)
For repeating this for several files I can
2007 Feb 27
4
fitting of all possible models
Hi,
Fitting all possible models (GLM) with 10 predictors will result in loads of (2^10 - 1) models. I want to do that in order to get the importance of variables (having an unbalanced variable design) by summing the up the AIC-weights of models including the same variable, for every variable separately. It's time consuming and annoying to define all possible models by hand.
Is there a
2007 Jan 11
3
batch job GLM calculations
Hello
I want to batch job the calculation of many GLM-models, extract some values and store them in a file. Almost everything in the script below works (read file, extract values and write them to file) except I fail in indexing the GLM with the modelstructure it should run. Running GLM's conventionally is no problem.
Conventionally a GLM is calculated as:
2000 Nov 16
2
newbee question
Dear All
Where can I lookup good methods to compute
p from q=bin(m,n)p^n*(1-p)^(m-n) such that
q<=alfa, alfa small. Are there such libs,
code and source in R?
Best Regards
--
Jan Burse SIAM, EAWAG
Scheuchzerstr. 67 ?berlandstr. 133
8006 Z?rich 8600 D?bendorf
tel: +41-1-364 17 66 tel: +41-1-823 55 34
2009 Nov 23
0
R-help Digest, Vol 81, Issue 23
Hi,
keine ahnung. Das liegt jetzt bei Hr. Feld. Frag mal bei ihm nach.
Gr??e
Thushyanthan
r-help-request at r-project.org wrote:
> Send R-help mailing list submissions to
> r-help at r-project.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> https://stat.ethz.ch/mailman/listinfo/r-help
> or, via email, send a message with subject or body 'help'