search for: cdu

...r message. Any help would be much appreciated. > germany2009 <- read.dta("StataforR.dta") > #attach(germany2009) > viewData(germany2009) > # > hcor <- function(data) hetcor(data, std.err=F, pd=T)$correlations > hetcor=hcor(germany2009) > hetcor cdupid MERKELLIKE CDULIKE MERKELPROX CDUPROX cdupid 1.0000000 0.5073577 0.6070514 -0.4064925 -0.4146810 MERKELLIKE 0.5073577 1.0000000 0.7832340 -0.4405571 -0.4200360 CDULIKE 0.6070514 0.7832340 1.0000000 -0.4868895 -0.5042006 MERKELPROX -0.4064925 -0.4405571 -0.4868895 1.0000000...

...b = sub, xlab = xlab, ylab = ylab, ...) : "data" is not a graphical parameter 8: In axis(1, ...) : "data" is not a graphical parameter 9: In axis(2, ...) : "data" is not a graphical parameter > searchpaths() [1] ".GlobalEnv" "/home/cdu/library/lattice" [3] "/home/cdu/library/stats" "/home/cdu/library/graphics" [5] "/home/cdu/library/grDevices" "/home/cdu/library/utils" [7] "/home/cdu/library/datasets" "/home/cdu/library/methods" [9] "Autoloads"...

HI, Dear R community, My original file has 1932 lines, but when I read into R, it changed to 1068 lines, how comes? cdu@nuuk:~/operon$ wc -l id_name_gh5.txt 1932 id_name_gh5.txt > gene_name<-read.table("/home/cdu/operon/id_name_gh5.txt", sep="\t", skip=0, header=F, fill=T) > dim(gene_name) [1] 1068 3 -- Sincerely, Changbin -- Changbin Du DOE Joint Genome Institute Bldg 400 Rm 45...

Hi I just upgraded my system to xen-3.0-testing, hoping it would solve my problems with the virtual network interfaceses (I cannot connect between two domains) and memory usage (I sometime "miss" a few hundret megabytes). The update didn''t help on the networking issues, but it did on memory usage. However, I created a new problem: using htop, the cpu is used at 50% in both

Dear community, I have one file named ca_boost_feature.txt, Feature selection (Boosting:0.0025,5)! H.2.C C.1.D C.3.R E.0.N C.2.S C.0.G H.3.G log file: ep If I want to use the second line of this file, how to read it into R? varr<-read.table("/home/cdu/operon/carbonic/ca_boost_feature.txt", sep=" ", skip=1, header=F, strip.white=TRUE, nrows=1) Warning message: In read.table("/home/cdu/operon/carbonic/ca_boost_feature.txt", : incomplete final line found by readTableHeader on '/home/cdu/operon/carbonic/ca_boost_featu...

...tand much of statistics. Therefore I'd like to give the reader an easy understandable example of Fleiss-Kappa in the Methods part. To achieve this, I obtained a table with the results of the German election from 2005: party number of votes percent SPD 16194665 34,2 CDU 13136740 27,8 CSU 3494309 7,4 Gruene 3838326 8,1 FDP 4648144 9,8 PDS 4118194 8,7 I want to show the agreement of voters measured by Fleiss-Kappa. To calculate this with the kappam.fleiss-function of irr, I need a...

...a111, b111, ..z111.. IF I dont know the relative position of the columns, but I know I need the following variables: var<-c(a1, c1,a11,b11,f111) Can I use read.columns to read the data into R? I have tried the following codes, but it does not work hh<-read.columns("/house/homedirs/c/cdu/operon/gh/dd.dimer", required.col=NULL, text.to.search=var, sep="\t", skip=0, quote="", fill=T) dim(hh) 468, 2402 hh<-read.columns("/house/homedirs/c/cdu/operon/gh/dd.dimer", required.col=var, text.to.search="", sep="\t", skip=0, quote=&q...

somfunc<- function (file) { aa_som<-scale(file) final.som<-som(data=aa_som, rlen=10000, grid=somgrid(5,4, "hexagonal")) pdf(file="/home/cdu/changbin/file.pdf") #output graphic file. plot(final.som, main="Unsupervised SOM") dev.off() } I have many different files, if I want output pdf file with the same name as for each dataset I feed to the function somfunc. HOw should I DO? THANKS! -- Sincerely, Changbin -- Ch...

...created a perspective plot using 'persp' in the graphics package. I'd like to add a second plane of z values to the existing plot, but I cannot seem to do this using 'persp'. Is there an analogue to 'lines' or 'points' for perspectives? Corey. corey.bradshaw at cdu.edu.au

Hey Guys, I just joined this list a couple days ago...what''s up with all the spam? I maintain about 8 or so googlegroups lists and have never received any <cross fingers>. Is this something to (unfortunately) look forward to? -- ______________________________________________________________ Clayton Dukes ______________________________________________________________

...lt;'!Y"E-E<FEE;#TX-S8U-#,R,0H*6T1R:79E($-="E!A=&@]+VUN="]W M:6Y?8PI4>7!E/6AD"DQA8F5L/4U3+41/4PI&:6QE<WES=&5M93UM<V1O<PH* M6T1R:79E($1="E!A=&@]+VUN="]C9')O;0I4>7!E/6-D<F]M"DQA8F5L/4-$ M+5)O;0I&:6QE<WES=&5M/7=I;CDU"@I;1')I=F4@15T*4&%T:#TO;6YT+W=I M;E]C+W=I;F1O=W,O=&5M<`I4>7!E/6AD"DQA8F5L/51M<"!$<FEV90I&:6QE M<WES=&5M/7=I;CDU"@I;1')I=F4@1ET*4&%T:#TD>TA/345]"E1Y<&4];F5T M=V]R:PI,86)E;#U(;VUE"D9I;&5S>7-T96T]=VEN.34*"EM...

...attached my output, the left one dont have any labels in the note, the right one has more than one label in each note. I need to have only one label for each note either by majority vote or prediction. Can anyone give some suggestions or advice? Thanks so much! alex<-read.table("/home/cdu/operon/alex2.txt", , sep="\t", skip=0, header=T, fill=T) alex1<-alex[,c(1:257)] levels(alex1$Label) alex1$outcome<-as.numeric(alex1$Label) alex1$outcome[1:20] #self-organizing maps(unsupervised learning) library("kohonen") #SOM, the supervised learning, train the...

...of you. Thanks in advance for your kind help. Have a good weekend. Cheers Xavier Xavier Hoenner PhD Student Nesting Ecology, Migrations and Diving Behaviour of Hawksbill Turtles School of Environmental Research, Charles Darwin University Darwin, NT 0909 Ph: (08)8946.7721 email:xavier.hoenner@cdu.edu.au [[alternative HTML version deleted]]

Hello colleagues, I am attempting to determine the nonlinear least-squares estimates of the nonlinear model parameters using nls. I have come across a common problem that R users have reported when I attempt to fit a particular 3-parameter nonlinear function to my dataset: Error in nls(r ~ tlm(a, N.fix, k, theta), data = tlm.data, start = list(a = a.st, : step factor 0.000488281

...lot using 'persp' in the graphics package. >I'd like to add a second plane of z values to the existing plot, but I >cannot seem to do this using 'persp'. Is there an analogue to 'lines' or >'points' for perspectives? > >Corey. >corey.bradshaw at cdu.edu.au > -- Romain FRANCOIS : francoisromain at free.fr page web : http://addictedtor.free.fr/ (en construction) 06 18 39 14 69 / 01 46 80 65 60 _______________________________________________________ Etudiant en 3eme ann??e Institut de Statistique de l'Universit?? de Paris (ISUP) Fili??re...

HI, Dear R-community, I AM using the following codes to grow tree and plot tree: # Classification Tree with rpart library(rpart) pdf(file="/home/cdu/changbin/dimer_tree.pdf") # grow tree fit.dimer <- rpart(outcome ~ ., method="class", data=p.dimer[,2:402]) plotcp(fit.dimer) # visualize cross-validation results # plot tree plot(fit.dimer, uniform=TRUE, main="Classification Tree for AA.dimer") text(fit.dimer, use.n...

the old one back in out of desperation and, BAM, all the problems went away. There were no other symtoms that pointed at that switch so we didn't suspect it all. Re: [Samba] Win98 troubles - Resolved Date: Fri, 08 Aug 2003 14:59:21 +0000 Message-ID: <20030808.cdU.52098400@ocicat.crandell.net> MIME-Version: 1.0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 8bit Content-Disposition: inline X-Mailer: AngleMail for phpGroupWare (http://www.phpgroupware.org) v 0.9.14.006 It turns out this all happened because of a brand...

Hello, After updating to CentOS-5.7, I have a (small) problem : The context of /dev/megadev0 is now defined (in /etc/selinux/targeted/contexts/files/file_contexts) as system_u:object_r:removable_device_t:s0. This cause smartmontools to fail : avc: denied { read write } for pid=2847 comm="smartd" name="megadev0" dev=tmpfs ino=8284

Dear R users, We have encountered a slight problem when using the lmer() function: 1. Data description: 11 locations; Nt: monthly mosquito population density from 1994-2005 in each location. 2. Question: to examine the degree of spatial heterogeneity in the system by testing model support for single versus multiple intercepts and slopes for the location effect. We applied the lmer()

hello all, I am trying to use the function getCovariateFormula(nlme) in conjunction with the library lme4. When I load both packages I get the following message and the getCovariateFormula function no longer works: library(nlme) library(lme4) Attaching package 'lme4': The following object(s) are masked from package:nlme : contr.SAS getCovariateFormula