Displaying 20 results from an estimated 1000 matches similar to: "lattice panel fine control"
2013 Feb 27
1
lattice xyplot point labelling
This is my reproducible example
tv.ms<-structure(list(inq = structure(4:17, .Label = c("D4", "D5", "D6a",
"D6b", "D6c", "D7", "D8", "F4", "F5a", "F5b", "F6a", "F6b", "F6c",
"F6d", "F7a", "F7b", "F8"), class =
2005 Jul 01
2
loop over large dataset
Hi All,
I'd like to ask for a few clarifications. I am doing some calculations
over some biggish datasets. One has ~ 23000 rows, and 6 columns, the
other has ~620000 rows and 6 columns.
I am using these datasets to perform a simulation of of haplotype
coalescence over a pedigree (the datestes themselves are pedigree
information). I created a new dataset (same number of rows as the
pedigree
2005 Dec 20
2
Extracting data from .zip file in WINDOWS version of package
Hello,
I am building a R-package for Genetics analysis. The accepted data is in pedigree (.ped) file format.
To load the data (say CAMP.ped) from "data" directory, I have a function "CAMP.R", which does the job.
The package builds successfully in Linux (.tar.gz) and the data loads successfully by "data(CAMP)".
However, when I build the package in WINDOWS, the data
2011 Apr 15
1
no solution yet, please help: extract p-value from mixed model in kinship package
I am making the question clear. Please help.
> Dear R experts
>
> I was using kinship package to fit mixed model with kinship matrix.
> The package looks like lme4, but I could find a way to extract p-value
> out of it. I need to extract is as I need to analyse large number of
> variables (> 10000).
>
> Please help me:
>
> require(kinship)
>
> #Generating
2011 Mar 22
1
help need on working in subset within a dataframe
Dear R-experts
Execuse me for an easy question, but I need help, sorry for that.
>From days I have been working with a large dataset, where operations are
needed within a component of dataset. Here is my question:
I have big dataset where x1:.....x1000 or so. What I need to do is to work
on 4 consequite variables to calculate a statistics and output. So far so
good. There are more vector
2005 Jun 02
3
How to change all name of variables
Dear R-helpers,
First I apologize if my question is quite simple
I have a large datasets which more 100 variables.
For a research I need to change all name of variables with add one or
more letters on each variables.
For example,
> data(Pima.tr)
> Pima.tr[1:5,]
npreg glu bp skin bmi ped age type
1 5 86 68 28 30.2 0.364 24 No
2 7 195 70 33 25.1 0.163 55 Yes
3 5
2010 Mar 18
2
Pedigree / Identifying Immediate Family of Index Animal
I have a data frame containing the Id, Mother, Father and Sex from about
10,000 animals in our colony. I am interested in graphing simple family
trees for a given subject or small number of subjects. The basic idea is:
start with data frame from entire colony and list of index animals. I need
to identify all immediate relatives of these index animals and plot the
pedigree for them. We're
2018 Jul 10
4
Construcción de archivo de texto
Hola a todos,
A partir de los siguientes datos:
d <- list(`1` = structure(list(ped = c(1L, 1L, 1L, 1L, 1L, 1L, 1L),
id = 1:7, father = c(2L, 0L, 0L, 2L, 2L, 2L, 2L), mother = c(3L,
0L, 0L, 3L, 3L, 3L, 3L), sex = c(2L, 1L, 2L, 2L, 2L, 1L,
2L), affected = c(1L, 2L, 1L, 1L, 2L, 2L, 2L)), row.names = c("1",
"2", "3", "4", "5",
2011 Sep 03
2
problem in applying function in data subset (with a level) - using plyr or other alternative are also welcome
Dear R experts.
I might be missing something obvious. I have been trying to fix this problem
for some weeks. Please help.
#data
ped <- c(rep(1, 4), rep(2, 3), rep(3, 3))
y <- rnorm(10, 8, 2)
# variable set 1
M1a <- sample (c(1, 2,3), 10, replace= T)
M1b <- sample (c(1, 2,3), 10, replace= T)
M1aP1 <- sample (c(1, 2,3), 10, replace= T)
M1bP2 <- sample (c(1, 2,3), 10, replace= T)
2005 Mar 18
1
How to show which variables include in plot of classification tree
Dear all
For my research, I am learning classification now.
I was trying some example about classification tree pakages, such as
tree and rpart, for instance,
in Pima.te dataset have 8 variables (include class=type):
library(rpart)
library(datasets)
pima.rpart <- rpart(type ~ npreg+glu+bp+skin+bmi+ped+age,data=Pima.te,
method='class')
plot(pima.rpart, uniform=TRUE)
text(pima.rpart)
2010 Feb 23
1
GenABEL - problems with load.gwaa.data
Hi all! I am using GenABEL on R for GWAS analysis. I am having a couple of
issues:
First, I am having a problem reading files (.map, & .ped, size 900Mb, using
windows 32-bit) onto R in the "convert.snp.ped" statement. I am thinking
this problem is likely due to the large size of the files & my version of R
is not able to handle them, since I can read in smaller files.
2009 Jan 22
1
infer haplotypes phasing trios tdthap
Dear R mailing list,
I have a dataset with genotypes from trios and I would like to infer
haplotypes for each mother, father and child. The package that I could
find that can do this is tdthap.
But when the mother is homozygous (e.g., 2/2) the haplotype is called as
not possible to infer (0); I would prefer for it to call the genotype
(2). From what I understand it is doing what I would like
2008 Feb 06
2
kinship package: drawing pedigree error
Hi
Im using the kinship package to draw a pedigree. On my data set this works fine but when i add indivudals to the pedigree i keep getting an error i hope someone can help me!
This is the code im using:
Data<-read.table("Tree.txt", header=T, sep=",")
attach(Data)
ped<-pedigree(id, dadid, momid, sex, aff)
par(xpd=T)
plot.pedigree(ped)
This is my data looks like
2011 Jun 21
4
Re; Getting SNPS from PLINK to R
I a using plink on a large SNP dataset with a .map and .ped file.
I want to get some sort of file say a list of all the SNPs that plink is
saying that I have. ANyideas on how to do this?
--
Thanks,
Jim.
[[alternative HTML version deleted]]
2010 Jun 30
6
Multiline and grouping in R
Hi All,
this is my first mail here.
I'm trying to plot a multiline chart grouping values with no success. I have
read a lot in the official Wiki and also searched via Google, but I did not
find anything.
I'm importing some data from a cvs file. Here is a sample:
YEAR,AREA,CASES
1988,CONTRACTS,286
1988,INTERNATIONAL,189
1988,FAMILY,385
1988,TAXATION,177
1989,CONTRACTS,233
2013 Feb 06
4
[Bug 60369] New: src/nouveau_exa.c:142:31: error: 'CREATE_PIXMAP_USAGE_SHARED' undeclared (first use in this function)
https://bugs.freedesktop.org/show_bug.cgi?id=60369
Priority: medium
Bug ID: 60369
CC: airlied at freedesktop.org
Assignee: nouveau at lists.freedesktop.org
Summary: src/nouveau_exa.c:142:31: error:
'CREATE_PIXMAP_USAGE_SHARED' undeclared (first use in
this function)
QA Contact:
2009 Jun 21
1
help with installation of local gzip-ped packages
I was suggested to install two tar-red and gzip-ped packages that are not part of CRAN or BioConductors yet.
I read the R manual about Administration and could only find a good description of how to install packages
not canonically included in CRAN repository, on UNIX systems.....
I work on Linux/SuSE and Windows ... so I am stuck with such an installation.
Any suggestion in very welcome.
Thank
2013 Oct 01
1
sonewconn: pcb 0xfffffe00c7223498: Listen queue overflow
Hello,
I updated our main server to 9.2-STABLE today and afterwards I noticed a
bunch of these messages, does anyone know what they mean? I was unable
to find anything on this error message. Things appear to be working OK
so far.
Sep 30 22:08:56 illidan kernel: sonewconn: pcb 0xfffffe00c7223498:
Listen queue overflow: 193 already in queue awaiting acceptance
Sep 30 22:12:27 illidan kernel:
2005 Sep 26
1
reading SAS data files
I am attempting to read in a SAS 9.1 data file. After starting R I
change to the directory containing the sas data file and use the "dir"
command to confirm that it is there. Then I run the following R-code:
library(foreign)
sashome <- "/Program Files/SAS/SAS 9.1"
test<-read.ssd(file.path(sashome), "pcb",
sascmd = file.path(sashome,
2012 Feb 14
2
how to test the random factor effect in lme
Hi
I am working on a Nested one-way ANOVA. I don't know how to implement
R code to test the significance of the random factor
My R code so far can only test the fixed factor :
anova(lme(PCB~Area,random=~1|Sites, data = PCBdata))
numDF denDF F-value p-value
(Intercept) 1 12 1841.7845 <.0001
Area 1 4 4.9846 0.0894
Here is my data and my hand