similar to: for help

Dear R-helpers: thank your for your attention. i am a newer to R and i am doing some protein category classification based on the amino acid sequence.while i have some questions urgently. 1. any packages for analysis amino acid sequence 2. given two sequences "AAA" and "BBB",how can i combine them into "AAABBB" 3. based on "AAABBB",how can i get some

I have a data set p1982, its structure is the following Then I take 20 observations from this dataset, and assign to pr. in p1982, p has 1982 levels, in dataset pr, p should have 1 levels. But I do str(pr), it shows that p still has 1982 levels. also for these > pr$aa [1] ARG THR ASP CYS TYR ASN VAL ASN ARG ILE ASP THR THR ALA SER CYS LYS THR ALA LYS Levels: ALA ARG ASN ASP CYS GLN

I have a data set like this I want to assign "outward" to Y if sc <90 and assign "inward" to Y if sc>=90. then cbind(p1982,Y) to get like these p aa as ms cur sc Y 1 154l_aa ARG 152.04 108.83 -0.1020140 92.10410 inward 2 154l_aa THR 15.86 28.32 0.2563560 103.67100 inward 3 154l_aa ASP 65.13 59.16 0.0312137 7.27311 outward 4 154l_aa CYS 57.20 49.85

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

I have a full data set like this: aa bas aas bms ams bcu acu omega y 1 ALA 0 127.71 0 69.99 0 -0.2498560 79.91470 outward 2 PRO 0 68.55 0 55.44 0 -0.0949008 76.60380 outward 3 ALA 0 52.72 0 47.82 0 -0.0396550 52.19970 outward 4 PHE 0 22.62 0 31.21 0 0.1270330 169.52500 inward 5 SER 0 71.32 0 52.84 0 -0.1312380 7.47528 outward 6

Hi, PDC works fine, but Samba BDC doesn't make its job. In srvmgr.exe PDC, BDC appear, but when I kill smb PDC's process, normaly BDC may give a response to smb request. My problem... BDC do not respond, no PDC :: no authentification. any idea. my smb.conf : [global] # Main Config. netbios name = LYS workgroup = TNN server string = Lys (TNN's PDC) security = user domain

Anybody interested in the S maps library for R? I have ported a limited version of the software for a Unix platform, and would be happy to make it available for anybody interested. If there is a lot of interest then I will submit it to CRAN. I know that Ross Ihaka is working on a new version for R, which does not have some of the limitations of the original S version. My version does not do

Dear expeRts, first of all I'd like to thank you for the quick help on my last which() problem. Here is another one I could not tackle: I have data on an absorption measurement which I want to fit with an voigt profile: fn.1 <- function(p){ for (i1 in ilong){ ff <- f[i1] ex[i1] <- exp(S*n*L*voigt(u,v,ff,p[1],p[2],p[3])[[1]]) } sum((t-ex)^2) } out <-

Hi, I am fairly new to 'R' and would like advice on the following. I want to calculate a weekly average number of reports (e.g. of flu, norovirus) based on the same weeks for the last five years. I will then use this to plot a chart with 52 points for the average based on the last five years; another line will then plot the current year, enabling a comparison of current weekly counts

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

I try to type this in my R-winEdt. but I got these. Do you know? Aimin > install.packages('http://rh-mirror.linux.iastate.edu/CRAN/bin/windows/contrib/2.4/plotrix_2.1-6.zip') Warning in download.packages(pkgs, destdir = tmpd, available = available, : no package 'http://rh-mirror.linux.iastate.edu/CRAN/bin/windows/contrib/2.4/plotrix_2.1-6.zip' at the

library(nlme) data(Phenobarb) na.include <- function(x)x phe1 <- nlme(conc~phenoModel(Subject, time, dose, lCl, lV), data = Phenobarb, fixed = lCl+lV~1, random= pdDiag(lCl+lV~1), start = c(-5,0), na.action = na.include, naPattern = ~!is.na(conc)) phe.ranef <- ranef(phe1,augFrame=TRUE) plot(phe.ranef, form=lCl~Wt+ApgarInd) [Error in max(length(x0),

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

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. Is

Hi all! I am desperately trying to figure out the solution to this error, but nothing as of yet is working. As noted in an earlier post I am using GenABEL. In an attempt to read in the phenotype file, in the format .dat, R keeps giving me the error "more columns than column names" I have tried to read in the data without the headers; I have also tried to trim the data to remove any

This is debugs seen on server, whose keys are not accepted by the client: debug1: Seeding random number generator debug1: sshd version OpenSSH_2.5.2p2 debug1: load_private_key_autodetect: type 0 RSA1 debug1: read SSH2 private key done: name rsa w/o comment success 1 debug1: load_private_key_autodetect: type 1 RSA debug1: read SSH2 private key done: name dsa w/o comment success 1 debug1:

Hello, I am new to R, How to delete a matrix column. Thanks,

I use knncat to make a predictive model and get misclass rate > knncat.m<-knncat(training.new,k=c(10,20),classcol=5) > knncat.m Training set misclass rate: 36.88% then I try to calculate prediction accuracy by the following: > pr.knncat.train <- predict (knncat.m,training.new,training.new,train.classcol=5,newdata.classcol=5) > tb.knncat.train <-table (pr.knncat.train,

I am trying to use scatterplot function in CAR like the following: scatterplot(X~Y) I want to label X points and Y ponits using the different color. Any idea for this? Aimin

Please ignore my ignorance, I'm new to R, but starting a longish project that will use R a lot. (a computer scientist learning about stats at the same time). Did anybody alreay write a XML parser for R? XML, as we will have tons of data-interchange with all sorts of other programs and XML is good for giving meaning to raw data. I checked the FAQs and the documentation, but didn't find