similar to: How to use pakcage SEM

Hi All. I am quite a newbie to R. This is a basic question. I like to modify elements of a vector. For Example: a1 <- c(1,2,3,4,3,5) TThe following program sentence does not work but the intention is; if (a1==3) a1*3 . 3 in the vector should be changed to 9, and the resulted vector is (1,2,9,4,9,5). How can I get the result? Thanks in advance for help. -------========-------- mitsu5

Hi All. I am quite a newbie to R. This is a next basic question. I have a matrix; > x <- matrix(1:10.,5) > x [,1] [,2] [1,] 1 6 [2,] 2 7 [3,] 3 8 [4,] 4 9 [5,] 5 10 I like to get a modified matrix as follows; [,1] [,2] [1,] 1 6 [2,] 2 7 [3,] 3 8 * 5 -> 40 [4,] 4 9 [5,] 5 10 The following expression does not work.

Hi. I am wondering whether Package SEM can do with intercepts and means in its structural analysis. If it can not calculate, how can I make a supplemental function in R? Many thanks in advance. --------========---------- Mitsuo Igarashi mitsu5 at ruby.famille.ne.jp

Hello, I'm trying to convert a character column in several dataframes to lower case. ### # # Sample data and 'spp' column summaries: # dput(ban.ovs.1993[sample(row.names(ban.ovs.1993), 20), 1:4]) ban.ovs.93 <- structure(list(oplt = c(43L, 43L, 38L, 26L, 35L, 8L, 39L, 1L, 34L, 50L, 10L, 29L, 31L, 24L, 18L, 12L, 27L, 49L, 28L, 51L), rplt = c(NA_integer_, NA_integer_,

Hi, I am trying to run samba with bind_dlz (bind-9.9.1 - P1) on a multi-homed network. I have configured the setup as per Samba4 Howto. But when I try to do "samba_dnsupdate --all-names" it fails with error: dns_tkey_negotiategss: TKEY is unacceptable The kerberos ticket being used by samba_dnsupdate shows follwoing principals: klist -c /tmp/tmp6cxfgY Ticket cache: FILE:/tmp/tmp6cxfgY

Dear R-users, I have the following problem. I have performed an experiment for which I gathered a lot of data which I now want to test. The problem is that I cannot find an appropriate test in R (I am a starter) and someone might give me a hand. This is what I have done: Across three sites (Site), I have laid out five transects (Trans)...meaning five transects in each sites. In each transect I

Dear Help-list, How can I do a multiple comparison test (mct) on selected contrasts from a linear model while using packages lme4 and multcomp? I am running R 2.13.0 under Windows 7. The following linear model and mct produces a global mct of 15 paired contrasts of the combined (Site, Position) factor SitePos of which only 9 are of interest. Model.G = lmer(log10(SrCa) ~ SitePos + (1 | Eel),

Hi there, I hope this question is not as stupid as the one before ... I tried to shorten my histogram (because the distribution is quite skewed and I simply don't want to see the long tail but still use the histogram plot). How can I do something like this? (The example does not work but I don't know why...) data <- rnorm(100) # as example, of course this is not skewed... h <-

Hi: In R, how do I display some intermediate results calculated in a "for" loop within a function? For example, in the attached code, how do I get it to print the intermediate variable "mh.new" for each simulation, when I call the function "MHsim.ind"? thanks for any help, Ravi. #################################################################### MHsim.ind

Hi, I found some weird behavior of the function multhist in the plotrix package and I would like to ask you about it: so, it works well to put two histograms next to each other, but I don't understand the breaks… so for example, library(plotrix) mh <- list(rnorm(200, mean=100, sd=10), rnorm(200, mean=100, sd=10)) multhist(mh) #produces a plot with

Dear Help List, Thanks in advance for reading...I hope my questions are not too ignorant. I have an experiment looking at evolution of wing size [centroid] in fruitflies and the effect of 6 different experimental treatments [treatment]. I have five replicate populations [replic] in each treatment and have reared the flies in two different temperatures [cond] to assay the wing size, making

Hi! I get the following message trying to install the matrix pakcage, can anyone help me please? trying URL `http://cran.r-project.org/bin/windows/contrib/2.0/Matrix_0.95-5.zip' Error in download.file(url, destfile, method, mode = "wb") : cannot open URL `http://cran.r-project.org/bin/windows/contrib/2.0/Matrix_0.95-5.zip' In addition: Warning message: cannot open: HTTP

Dear list, I have following data set, which I want to plot the "Scale" variable on the x-axis and "Mean"´on the y-axis for each Ageclass and for each sex. The Mean value of each Ageclass for each sex would be connected by a line. Totally, there should be 6 lines, from which three present the Mean values of each Ageclass for respective sex. Are there any easy ways to do

Dear list, Is there a way to reorder the xaxis using lattice. Using the following data, the x-axis is ordered as BP GH MH PF RE RP SF VT but I would like the x-axis to be ordered as PF RP BP GH VT SF RE MH. Kön Skalor Tillfälle Medelvärde 1 Kvinnor BP 1-inskrivning 36.45283 2 Kvinnor GH 1-inskrivning 38.62255 3 Kvinnor MH 1-inskrivning 62.88889 4

Hi! I am trying to print from an NT machine to a simple samba printer. That NT machine is the PDC and the samba box is having logins authenticated by the NT box. When I try to access the printer with my own username, it works fine. As soon as I have the account that is destined to do the printing print, it can't log in. I have included my smb.conf file and two log excerpts - one successful,

Dear all Sorry for simple question: I want to put the following option into look as number of X is large 1000 variables X1 <- sample(c(1,2, 3, 4),10, replace = T, prob = c(0.4, 0.2, 0.2, 0.2)) cv1 <- round(runif(2, 1, 10)) # X2 is copy of X1 X2 <- X1 # now X2 is different in cv1 random positions X2[cv1] <- 5 cv2 <- round(runif(2, 1, 10)) # X3 is copy of X2 X3 <- X2

Dear R users, In metafor 1.6-0, the Mantel-Haenszel method is implemented by the rma.mh() function. I have observed that the sum of the weights computed by weights(x) doesn't add to 100% when x is an object of class rma.mh. The consequences of this fact can be clearly seen when a forest diagram is drawn with forest(x), which calls weights(x) (or more precisely, the method weights.rma.mh()

Hi, I have a function that generates a set of data but I am having problems determining the parameters using the nls fitting procedure. #### "MH"<-function(field,diameter,mu=10e-7,sig=0.1,Ms=100,chi=0){ #variables mu, sig, chi, Ms #input: field and diameter #all in CGS rho <- 5 kb <- 1.38e-16 t <- 300 length.d<-length(diameter) length.H<-length(field)

Hello All, I have a technical difficulty installing RMySQL. I am running openSUSE11.1 and R 2.12 I have installed MySQL from the website. and following installation , as root This is the part where trouble begin, ...... checking mysql.h usability... no checking mysql.h presence... no checking for mysql.h... no checking for mysql_init in -lmysqlclient... no checking for mysql_init in

Hi! This is my share: |[root] | comment = Root-Partition | path = / | guest ok = no I have a file called /etc/exim/config that is required to be owned by mail.mail, mode 664. My user account is called mh and I am member of the group mail. I'd like to edit that file via samba with an editor that is running on my desktop Windows box. Whenever I write to the file, it gets reset to mh.mh 744