Displaying 20 results from an estimated 2601 matches for "p5".
2006 Dec 08
1
please help me for 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))...
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))...
2006 Nov 23
2
random effect question and glm
consider p as random effect with 5 levels, what is difference between these
two models?
> p5.random.p <- lmer(Y
~p+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
> p5.random.p1 <- lmer(Y
~1+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
in addtion, I try these two models, it seems they are same.
what is the difference between these two model. I...
2006 Nov 20
4
for help about logistic regression model
...50 75.3691
5 821p LYS 5 -4.27189 27.7594 6.272780 1.205650 1.20123 1.633780
53.3304 41.98 57.68 0.1305950 91.1431
6 821p LEU 6 0.05675 27.5178 6.309750 1.370120 0.64664 1.656920
27.4681 0.00 0.00 0.0000000 94.0851
here p is random effect, and aa is nested in p
I do like this:
p5 <- read.csv("p_5_angle.csv", header=T, sep=",")
Y<-p5$sc>=90 # probability of pointing inward
library(MASS)
mp5.null <- glmmPQL(Y~1,data=p5,random=~1|p/aa,family=binomial(logit))
summary(mp5.null)
mp5.full<-glmmPQL(Y~as*ms*cur,data=p5,random=~1|p/aa,family=binom...
2012 Feb 13
3
Change dataframe-structure
Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt.
Name: nicht verf?gbar
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120213/d2a5afa6/attachment.pl>
2008 Feb 04
2
make dataframe from table
Dear R-experts,
I have got a dataframe:
data
ID disease
V1 V2
1 p1 1
2 p1 3
3 p3 3
4 p3 5
5 p5 1
From which I extract a usefull table: affect
affect
1 3 5
p1 1 1 0
p3 0 1 1
p5 1 0 0
I want to merge this with anotherdataframe:
age
p1 23
p2 24
p3 23
p4 11
p5 45
If have tried as.data.frame(affect) and other solutions to get the
following comment going:
merge(age,affec...
2016 Apr 19
2
Creating two new variables conditional on retaining values from previous rows
...1 and groupdate = date or else if not first ID and date - previous date <= 10 or date - previous group date<=10 then group=previous group # and groupdate = previous date.
Input:
ID DATE ITEM
1 1/1/2014 P1
1 1/15/2014 P2
1 1/20/2014 P3
1 1/22/2014 P4
1 3/10/2015 P5
2 1/13/2015 P1
2 1/20/2015 P2
2 1/28/2015 P3
2 2/28/2015 P4
2 3/20/2015 P5
Desired Output
ID DATE ITEM GROUP GROUPDATE
1 1/1/2014 P1 1 1/1/2014
1 1/15/2014 P2 2 1/15/2014
1 1/20/2014 P3 2 1/15/2014
1 1/22/2014 P4 2 1/15/2014
1 3/10/20...
2011 Aug 13
1
Own R function doubt
...Pi*Hembras)
P1<-Pi-B
R1<-P1*Hembras*R
M1<-P1*M
T2<-(R1-M1+B)/(P1-M1+R1)
P2<-P1-B
R2<-P2*Hembras*R
M2<-P2*M
T3<-(R2-M2+B)/(P2-M2+R2)
P3<-P2-B
R3<-P3*Hembras*R
M3<-P3*M
T4<-(R3-M3+B)/(P3-M3+R3)
P4<-P3-B
R4<-P4*Hembras*R
M4<-P4*M
T5<-(R4-M4+B)/(P4-M4+R4)
P5<-P4-B
R5<-P5*Hembras*R
M5<-P5*M
T6<-(R5-M5+B)/(P5-M5+R5)
P6<-P5-B
R6<-P6*Hembras*R
M6<-P6*M
T7<-(R6-M6+B)/(P6-M6+R6)
P7<-P6-B
R7<-P7*Hembras*R
M7<-P7*M
T8<-(R7-M7+B)/(P7-M7+R7)
P8<-P7-B
R8<-P8*Hembras*R
M8<-P8*M
T9<-(R8-M8+B)/(P8-M8+R8)
P9<-P8-B
R9...
2006 Nov 22
2
help
consider p as random effect with 5 levels, what is difference between these
two models?
> p5.random.p <- lmer(Y
~p+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
> p5.random.p1 <- lmer(Y
~1+(1|p),data=p5,family=binomial,control=list(usePQL=FALSE,msV=1))
thanks,
Aimin Yan
2006 Dec 08
0
svm code, what is wrong here?
> install.packages("e1071")
Warning: package 'e1071' is in use and will not be installed
> library(e1071)
> library(MASS)
> p5 <- read.csv("http://www.public.iastate.edu/~aiminy/data/p_5_2.csv")
> attach(p5)
The following object(s) are masked from p5 ( position 3 ) :
Aa As Cur Ms P Y
The following object(s) are masked from p5 ( position 4 ) :
Aa As Cur Ms P Y
&g...
2011 Dec 23
1
FreeBSD Security Advisory FreeBSD-SA-11:09.pam_ssh
...Category: contrib
Module: pam
Announced: 2011-12-23
Credits: Guy Helmer, Dag-Erling Smorgrav
Affects: All supported versions of FreeBSD.
Corrected: 2011-12-11 20:40:23 UTC (RELENG_7, 7.4-STABLE)
2011-12-23 15:00:37 UTC (RELENG_7_4, 7.4-RELEASE-p5)
2011-12-23 15:00:37 UTC (RELENG_7_3, 7.3-RELEASE-p9)
2011-12-11 20:38:36 UTC (RELENG_8, 8.2-STABLE)
2011-12-23 15:00:37 UTC (RELENG_8_2, 8.2-RELEASE-p5)
2011-12-23 15:00:37 UTC (RELENG_8_1, 8.1-RELEASE-p7)
2011-12-11 1...
1997 Dec 30
0
Weird problem with smb 1.9.17p4 and p5
Overall environment is SunOS and Solaris 2.X machines running a mixture of
smb.1.9.17.p4 and p5.
This machine is a SunOS 4.1.3_U1 running smb.1.9.17p5. (debug level 3)
AAAA is Primary Domain Controller (NT4.0) 206.63.206.24
zotz is Win95 machine
>From smb.conf:
security = server
password server = aaaa
>From log.zotz:
[snip]
Connecting to 206.63.206.24 at port 139
connected to pas...
2003 Jun 19
2
Subseting by more than one factor...
Is it possible in R to subset a dataframe by more than one factor, all at
once?
For instance, I have the dataframe:
>data
p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 pred
1 0 1 0 0 0 0 0 0 0 0 0.5862069
4 0 0 0 0 0 0 0 0 0 1 0.5862069
5 0 0 0 0 0 0 1 0 0 0 0.5862069
6 0 0 0 0 0 0 0 1 0 0 0.5862069
7 0 0 1 0 0 0 0 0 0 0 0.5862069
9 0 0 0 0 1 0 0 0 0 0 0...
2011 Dec 19
1
calculating correlation coefficients on repeated measures
...9 repeated measures (measurement variable == 'Delta13C') for
individuals (ID variable == 'Individual_ID'. Each repeated measure is
"indexed" (right term?) by the variable 'FeatherPosition' and given as
c('P1', 'P2', 'P3', 'P4', 'P5', 'P6', 'P7', 'P8', 'P9'). I would like
to calculate a correlation coefficient (r) and p.value for all
measures of 'Delta13C' by individual. the function 'cor' only seems to
work when comparing two individual measures (e.g. P1 and P2, P2 and
P3, e...
2004 Jun 10
2
odesolve: lsoda vs rk4
...ningful?
Any help much appreciated,
Thanks in advance,
Chris
func <- function(t, y, p)
{
Ad <- p["p2"]*(p["p1"]*y["A"]*y["D"])/(p["p2"]+p["p3"]) +
p["p6"]*(p["p4"]*y["B"]*p["p10"])/(p["p5"]+p["p6"]) -
p["p1"]*y["A"]*y["C"]
Bd <- p["p3"]*(p["p1"]*y["A"]*y["D"])/(p["p2"]+p["p3"]) +
p["p5"]*(p["p4"]*y["B"]*p["p10"])/(p["p5"]+p[&qu...
2004 May 06
5
Orthogonal Polynomial Regression Parameter Estimation
...The following values are coefficients of 0-6th order
(for n=8) polynomial collected from Pearson, Hartley
(1958) table, page 212:
> p0<-rep(1,8)
> p1<-c(-7,-5,-3,-1,1,3,5,7)
> p2<-c(7,1,-3,-5,-5,-3,1,7)
> p3<-c(-7,5,7,3,-3,-7,-5,7)
> p4<-c(7,-13,-3,9,9,-3,-13,7)
> p5<-c(-7,23,-17,-15,15,17,-23,7)
> p6<-c(1,-5,9,-5,-5,9,-5,1)
Now, the estimated parameters of the orthogonal
polynomial is calculated by the following formula:
> alpha0<-sum(Y*p0)/sum(p0^2);
alpha1<-sum(Y*p1)/sum(p1^2);
alpha2<-sum(Y*p2)/sum(p2^2);
alpha3<-sum(Y*p3)/sum(p3^2...
2011 Feb 02
2
Help with one of "those" apply functions
Hello there,
I'm still struggling with the *apply commands. I have 5 people with id's
from 10 to 14. I have varying amounts (nrep) of repeated outcome (value)
measured on them.
nrep <- 1:5
id <- rep(c("p1", "p2", "p3", "p4", "p5"), nrep)
value <- rnorm(length(id))
I want to create a new vector that contains the sum of the values per
person.
subject.value[1] <- value[1] # 1 measurement
subject.value[2] <- sum(value[2:3]) # the next 2 measurements
...
subject.value[5] <- sum(value[11:15]) # the next 5...
2011 Feb 09
3
Loop in variable names
...I should also point out that my data set has several missing observations and sometimes the observations that are missing are not the same for all my variables.
What I would like to do:
table(StoreData$CompanyID,
!is.na(StoreData$P2))
table(StoreData$CompanyID,
!is.na(StoreData$P5))
If I run the above code, I get:
> table(StoreData$CompanyID,
+ !is.na(StoreData$P2))
FALSE TRUE
2 940 0
3 0 323
4 288 0
5 306 0
> table(StoreData$CompanyID,
+ !is.na(StoreData$P5))
FALSE TRUE
2 940 0
3 0 323
4...
2012 Mar 16
1
R merge two dataframes with different row?
...101 name1
2 200 201 name2
2 300 301 name3
3 400 401 name4
3 500 501 name5
4 600 601 name6
4 700 701 name7
File2 is:
V1 V2 V3 V4
1 50 55 p1
3 402 449 p2
4 550 650 p3
4 651 660 p4
2 150 250 p5
2 250 350 p6
3 450 499 p7
2 100 250 p8
I hope to have the merged file3 meet the following three criteria:
(1) File1$V1==File2$V1, and
(2) File1V2>=File2$V2, and
(3) File1V3<=File2$V3.
In this case, we can see that there should be four records meet these three
crit...
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 1 1 1 1 1 1 1 1 1 ...
$ Aa : Factor w/ 19 levels "ALA","ARG",&q...