search for: pauljohn

...39;m using Ubuntu-9.04. I'm trying to profile an package, and got errors, so I decided to follow BDR's (http://tolstoy.newcastle.edu.au/R/devel/06/02/4254.html) example exactly, and got same problem. Can anyone replicate? $ sudo mkdir -p /var/tmp//usr/lib/R/library/stats/libs $ chown -R pauljohn.pauljohn /var/tmp/usr/lib/R $ export LD_PROFILE=/usr/lib/R/library/stats/libs/stats.so $ R ### in R session, type >example(smooth.spline) > q() Observe that does create the profile file: $ ls -la /var/tmp/usr/lib/R/library/stats/libs total 20 drwxr-xr-x 2 pauljohn pauljohn 4096 2009-0...

I'm working up a set of small working examples in R to show what various distributions are and the beauty of the central limit theorem. Those example programs are in this directory: http://lark.cc.ukans.edu/~pauljohn/R/ExampleCode/ You can feel free to use those if you want, or you can send me other small working example code for R. The key here is small, self contained things that show things that R can do and people can run without too much trouble. Anyway, I have hit these problems (questions) 1. When...

...ograms, the results will typically include N and the number dropped because of missings. Without going back to R and fiddling about to find the total number of rows in the dataframe, there is no way to tell. Somewhat inconvenient. Do you agree? -- Paul E. Johnson email: pauljohn at ku.edu Dept. of Political Science http://lark.cc.ku.edu/~pauljohn 1541 Lilac Lane, Rm 504 University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66044-3177 FAX: (785) 864-5700

I have a linux system with Fedora Core 2 and R-2.0. I was comparing plots made with plot() and dotplot() and discovered a problem. Although the dots are positioned correctly, the numerical labels in the dotplot y axis are not correct. I put copies here: http://lark.cc.ku.edu/~pauljohn/R/plotTrouble1.jpg That is the "correct" one from plot, with the higest value on y showing at 18. http://lark.cc.ku.edu/~pauljohn/R/plotTrouble2.jpg That is the dotplot one. The picture is basically the same, except the numbers on the y axis only go up to 8. But the dots are in the...

...'s something about the environment that I need to tell the function, I just don't know what? pj ps. Now that Frank Harrell is using R, I'm so enthusiastic! I'll do some R in my advanced methods class in Fall, 2001 and continue my updating of the tipsheet: http://lark.cc.ukans.edu/~pauljohn/R/statsRus.html -- Paul E. Johnson email: pauljohn at ukans.edu Dept. of Political Science http://lark.cc.ukans.edu/~pauljohn University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66045 FAX: (785) 864-5700 -.-.-.-.-.-.-.-.-.-...

https://bugzilla.samba.org/show_bug.cgi?id=3844 Summary: --delete with --dry-run no longer reports files to be deleted Product: rsync Version: 2.6.8 Platform: x86 OS/Version: Linux Status: NEW Severity: normal Priority: P3 Component: core AssignedTo: wayned@samba.org

...ith > postscript("../output.ps") and then run the commands again, and then > dev.off() null device 1 In the ps file that is saved, all of the lines look the same. I used the Gimp to take a "screenshot" of the figure as it looks in R: http://lark.cc.ukans.edu/~pauljohn/exampleRGraph.gif And the postscript output looks like this: http://lark.cc.ukans.edu/~pauljohn/exampleRGraph.ps Any help in saving graphs would be appreciated. I feel totally clueless, and I've been looking pretty hard in the V&R and R-intro. 3. In light of difficulty with line styles...

I have made a plot with a legend and on the screen it looks fine, but when I save as jpg or pdf, or print, the legend box is too small, it cuts through the words on the right side. I put an example here: http://lark.cc.ukans.edu/~pauljohn/R/apdftest.pdf Is there a work around? -- Paul E. Johnson email: pauljohn at ukans.edu Dept. of Political Science http://lark.cc.ukans.edu/~pauljohn University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66045 FAX: (785) 864...

...es],sd) dir<-system("pwd",TRUE) newOutput<-rbind(mean,sd) transNewOutput<-t(newOutput) transNewOutput<-rbind(dir,transNewOutput) write.table(transNewOutput,file="Atestoutput",quote=FALSE, sep="\t") The output file looks like this: mean sd dir /home/pauljohn/swarm/src /home/pauljohn/swarm/src T 732.1 148.618520536522 Run 49.5 29.011491975882 Seed 499366074.13 296739559.579976 IntAct 43.85 5.592879672432 Change 0 0 Now the questions: 1. I want 3 columns of output, but here the first line only has two items. I can&...

...a predictor in a logistic regression, then it is tough to plot the S-shaped curve that describes the probabilities. Here is some code with comments describing the difficulties that we have found. Are there simpler, less frustrating approaches? ----------------------------- # Paul Johnson <pauljohn at ku.edu> 2005-11-17 # factorFutzing-1.R myfac <- factor(c(1.1, 4.5, 1.1, 1.1, 4.5, 1.1, 4.5, 1.1)) y <- c(0,1,0,1,0,0,1,0) mymod1 <-glm (y~myfac, family=binomial) p1 <- predict(mymod1, type="response") plot(myfac,p1) # Wait, that's not the picture I wanted. I w...

...If this is a bad idea, you can tell me, my feelings won't be hurt :) pj ps. I just found that SAS has at least one procedure for this, called tpsplines (thin-plate splines), so I know I wasn't misunderstanding this fellow's lecture. -- Paul E. Johnson email: pauljohn at ukans.edu Dept. of Political Science http://lark.cc.ukans.edu/~pauljohn University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66045 FAX: (785) 864-5700 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing...

...e, or what difficulties might be caused with non-English character sets and file systems. Not to mention the problem that the documentation/manuals are in English. If you have any advice that I can collect up for my friends, I would appreciate it. -- Paul E. Johnson email: pauljohn at ku.edu Dept. of Political Science http://lark.cc.ku.edu/~pauljohn 1541 Lilac Lane, Rm 504 University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66044-3177 FAX: (785) 864-5700

...: I have Fedora Core 2 and R-1.9.0 does not build "out of the box". After applying the daily patch file R-release.diff, I find it does build and I've made RPMS and SRPM. In case you want to save yourself a recompile, you can find Fedora Core RPMs in here: http://lark.cc.ku.edu/~pauljohn/software/R These are based on the standard R distribution SPEC file, only the patch was added. Also, I have created an RPM for ESS for Emacs on Fedora Core 2 (not Xemacs) and you can find that here: http://lark.cc.ku.edu/~pauljohn/software/favoriteEmacsFiles -------------------------- While...

...f5 which allows loading and saving of data in hdf5 format. The source for his library is here: ftp://ftp.swarm.org/pub/swarm/src/testing/hdf5_1.4.7.tar.gz I prefer to use RPMs to maintain software across many linux systems, so I have the RPM for this R library here: http://lark.cc.ukans.edu/~pauljohn/software/R-hdf5-1.4.7.R1.4.0-1.i386.rpm The SRPM is in the same directory. Please note in order to use this you have to have hdf5 (duh!) and I have the rpm for the required version 1.4.2 in that same software directory. Because of slight changes in RH7.2 in the ssl libraries, it is necessar...

...hat seems bad to me. I see the lag function in ts, but when I use it, it doesn't change x, so obviously I don't see the point of that. > z <- lag (x) > z [1] 1 3 1 3 1 3 1 3 4 4 So long, thanks in advance, greetings, etc... -- Paul E. Johnson email: pauljohn at ukans.edu Dept. of Political Science http://lark.cc.ku.edu/~pauljohn University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66045 FAX: (785) 864-5700 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing l...

...t;,height=4,width=3.5) I was surprised that the output in the eps file included the colored lines from the plot, even though the screen device was set to "gray" and on screen I did not see colors. I put some of the eps files up here, so you can see what I mean: http://lark.cc.ku.edu/~pauljohn/R/epsfiles And if you don't have an eps viewer handy, I made png versions, since all browsers seem to show them. http://lark.cc.ku.edu/~pauljohn/R/pngfiles/ The docs indicate that the x11 device is being copied, but I rather think that a plot object that exists "behind" the x11 d...

...ry users can read. Here's the shell script I will schedule with cron ------------R_installAll.sh ---------- #!/bin/bash R CMD BATCH /usr/local/bin/R_installAll.R -------------------------------------- And here's the R program -------------R_installAll.R------------- # Paul Johnson <pauljohn _AT_ ku.edu> 2005-08-31 # This should update and then install all packages, except for # ones I exclude because they don't work or we don't want them. options(repos = "http://lib.stat.cmu.edu/R/CRAN/") update.packages(ask=F) theNew <- new.packages() failPackages <- c(...

...e general social survey with Stata8 for linux. The file can be re-opened with Stata, but when I bring it into R, it says all the values are missing for most of the variables. This dataset is called "morgen.dta" and I dropped a copy online in case you are interested http://www.ku.edu/~pauljohn/R/morgen.dta looks like this to R (I tried various options on the read.dta command): > myDat <- read.dta("morgen.dta") > summary(myDat) CASEID year id hrs1 hrs2 Min. : 19721 Min. :1972 Min. : 1 NAP : 0 NAP :...

...lues are replaced by their lables. In R, does one keep a list of variable lables and value lables? If so, where? It seems to me that an R collection type like Map or Hash would be needed so one could put in a variable name and get back a lable. -- Paul E. Johnson email: pauljohn at ukans.edu Dept. of Political Science http://lark.cc.ukans.edu/~pauljohn University of Kansas Office: (785) 864-9086 Lawrence, Kansas 66045 FAX: (785) 864-5700 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailin...

...ry to do something a little different, I cannot get output saved just right. This is on RedHat 7.2 with R 1.4.1. The png output looks fine, but the eps output has the problem that the bounding box on the legend cuts the legend in half. I put a copy of a bad one here: http://lark.cc.ukans.edu/~pauljohn/ResearchPapers/meanProtest-box.eps When I asked about these chores about 6 months ago, I got some tips which have worked other times. I had the impression that if I set the size of the x11 device to match the eventual ps output size, then all would be OK. But its not this time. Here's my c...